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eight-page supplements giving details of reproductions of German technical 
publications issued under the authority of the Alien Property Custodian 
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SECOND-HAND SCIENTIFIC BOOKS. An extensive stock of books in 
all branches of Pure and Applied Science may be seen In this department. 
Large and small collections bought. Back volumes of Scientific Journals. 
Old and rare Scientific Books. Mention Interests when writing. 

140 GOWER STREET. 


SCIENTIFIC LENDING LIBRARY 


Annual subscription from One Guinea. _ Details of terms and prospectus free on request. 
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of authors and subjects: to subscribers !2s. 6d. net., to non-subscribers 25s. net., postage 8d. 
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Telegrams: “ Publicavit, 
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Tree ANNOUNCEMENT 
RECTIFIERS Megatron Selenium 
DOUBLE-VOLTAGE Fnoto: Cells 
WESTALITE RECTIFIERS oe 

(Size and weight halved) in PRODUCTION 
| M INTA TU RE Technical Enquiries to 
METAL RECTIFIERS Dp Gr A; VESZI; 
CONSTANT 15a, Sateen N.4 
POTENTIAL A.C. Archway 3739 
(THE STABILISTOR) eens 
ATIC PHASE CONVERSION! | MEGATRON LTD. 


are a few of the latest 


96/9/,.CHANDOS. = HOUSE: 
BUCKINGHAM GATE, S.W.|I 
Abbey 7216 


DEVELOPMENTS 


Miniature ‘‘ Westector ”’ 


[se Z/E 
FOR INDICATING - RECORDING: CONTROLLING 2° 0 | E 
TEMPERATURE “= 


We shall be glad to supply full particulars 


BOWEN INSTRUMENT Co, Ltd. 
9/13, NEWTON ROAD - LEEDS, 7 


Telephones: 41036/7 


“Stabilistor '' 
STINGHOUSE BRAKE & SIGNAL Co. Ltd. 
Pew Hill House, Chippenham, Wilts. 
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VITREOSIL | 
APPARATUS _ 


AVIMO HIGH PRESSURE COUPLINGS 


The only really efficient pressure-resisting coupling 
for pipe connections. Approved by the Air Ministry 
and standarised by leading aircraft manufacturers. For 
water, oil and glycol, at any temperature. Flexible. 
Vibration proof. Available in sizes from 34” o/d to 3}” o/d. 
Avimo facilities are at your disposal. 


This photograph—by the way, it won |] 
International Exhibition in 1939—shows | 
some pre-war apparatus made to speci- | 
fication in VITREOSIL pure fused silica, 
but you need hardly be reminded that |} 
DESIGN : DEVELOPMENT - MANUFACTURE this work can always be done for you. |f} 


OPTICAL-MECHANICAL TRICAL INSTRUMENTS THE THERMAL SYNDICATE LTE 
gps Wirt hihed tema 2 NereOne Head Office: Wallsend, Northumberland. 


London Depot: 12-14, Old Pye Street, 
Westminster, S.W.|1. 


SMALL GEARED 
MOTOR UNITS 


Drayton ““ RQ” motors are available unidirectional 

or reversing with or without self-switching for 100 

110 or 200 250 Volts A.C. 

FIL NASHAFT SPEEDS : 

Reversing type. ... Maximum 600 rev./min. ; 
minimum 27 min./rey. 

Unidirectional type Maximum 280 rev./min. ; 
minimum 60 min./rev. 


TORQUE 
Reversing ce ... 60 1b. in. maximum. 
Unidirectional ... ... 37 1b.in. maximum. 


TYPICAL APPLICATIONS :—Operating valves, dampers, 

or rheostats actuating Geneva movements for drum- 
type switches, rocking baths, continuous turning, feed 
of light strip under process, work movement, eg. | 
soldering and welding fixtures. 


The 


Drayton Regulator & Instrument Co Lid. 
West Drayton : : Middiese 
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LABORATORY FURNITURE 


Functional suitability, 
technical convenience 
and personal  con- 
venience are three 
important factors to 
be observed in design- 
ing a laboratory. 

It is evident that the 
successful planning of 
a laboratory requires 
the exercise of a 
degree of technical 
knowledge which few 
fully appreciate, It 
involves the solution 
of complex design 
problems (ventilating, 
fume removal, drain- 
age, layout, disposal of 
supply services) and a 
thorough knowledge 
of structural materials 
and their suitability 
for use in various parts 
of the laboratory. 


We are experienced 
designers of laboratory 
furniture, manufac- 
tured in our own 
works. 


Grrr INand lATLOCK Ly 


Established as Scientific Instrument Makers in 1826 


LONDON MANCHESTER GLASGOW EDINBURGH 
Kemble Street, W.C.2 19, Cheetham Hill Rd., 4 45, Renfrew Street, C.2 7, Teviot Place, | 


BIRMINGHAM : STANLEY BELCHER & MASON LTD., CHURCH STREET, 3 


ie oe 


| i i itter T h 
A ree © Model H.s. 300 
HIGH SPEED CAMERA 


Available on hire for high speed research 
work requiring accurate observation of 


fast moving machinery, projectiles, etc. 


a 


= a 


si 


logue 
F 


cation 


> 
junicate 


rements 


nber of 
schnical 
1 staff 
isit your 
y 

Bogient. 


COST OF HIRE:- 
5.15.0 PER DAY (including operator) 


W. VINTEN LIMITED 


North Circular Road, Cricklewood, N.W.2. 
Telephone : GLlAdstone 6373 


a2 
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SELENIUM PHOTOCELLS 


for use with direct reading instruments or valve amplifiers. 


TEST TUBE COLORIMETERS 


for Haemoglobin and other determinations. 


REFLECTOMETERS 


for comparing reflectance from all surfaces. 


EVANS ELECTROSELENIU 


LIMITED 
BISHOP'S STORTFORD. HERTS. 


REGD, 


MAGNETIC SCREENS 


MUMETA| 


MUMETAL is a renowned TELC 


Passed for the shield; ON alloy uneur- 


5 of delic : 
one uni-directional ae ate equipment 


—— fluxes. D 


alternating magnetic 


uring Wartime MU 


METAL has been 


* The experience 


your peacetim 


inquiries are invite d © problems and 


THE TELEGRAPH CONSTRUCTION & MAINTENANCE CO, 
Head Office: 22 OLD BROAD ST., LONDON, E.C.2. Tel: LONdon We 


| 
Founded 1864 
| 


Enquiries to: TELCON WORKS, GREENWICH, S.E.10. Tel: Greenwie¢ 
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\/ ) \V FOR NEARLY A CENTURY WRAY LENSES 
i\ HAVE ENJOYED THE CONFIDENCE OF SCIENTISTS 


SCIENTIFIC 
LENSES 


Vv 


THE LENS’ HERE ILLUSTRATED 
HAS BEEN SPECIALLY DESIGNED 
FOR USE ON CONSTANT SPEED 
CAMERAS TO RECORD ON 35 mm. 
FILM FROM A CATHODE RAY TUBE. 
IT IS AN ANASTIGMAT OF 2” FOCAL 
LENGTH, HAS AN APERTURE OF F.I. 
AND IS CORRECTED FOR AN OBJECT 
TO IMAGE RATIO OF 4:1. 


RAY (OPTICAL WORKS) LTD. BROMLEY, KENT 


| SIMPLIFIED 


BA IN 


"INSTRU MENT COMPANY LIMITED 


| MEASUREMENT 


Laboratory methods of determination requiring time 

and skill can often be replaced by direct-reading 

instruments that can be successfully used by the non- 

technical worker. Our development service for industry 

covers the field of Idea, Prototype and Manufacture. 
Ask us about it. 


A AND 


; CONTROL 


BALDWIN INSTRUMENT COMPANY LIMITED 
CUMNOR - OXFORD 


riginators and makers of scientific instruments for measurement and control. 
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TEMPERATURE and HUMIDITY-CONTROL 
with Adjustable Contact Thermometers 


AUTOMATIC VACUUM-CONTROL for 
LABORATORY and ALTITUDE TESTS 


SUPER-SENSITIVE RELAYS and 
MERCURY SWITCHES 


COMPONENTS for HERMETIC SEALING 


LABGEAR 


MANUFACTURERS OF 
Laboratory Apparatus 
scientific Instruments & 
Electronic Equipment, etc. 


INVITE YOUR INQUIRIES FOR SPECIAL REQUIREMENTS OR CATALOGUES 


WRITE OR PHONE LABGEAR PHONE : 2494/5 
WILLOW PLACE CAMBRIDGE 


i | 
} 
ML) | 


ELECTRO METHODS Ltd. 12 Svustaoy4. en 


| 
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A typical G.E.C. Cathode 
Ray Tube 4102, with 237 
screen, widely used in 
industry during the war, 
and with many peace-time 
applications. 


SF 


PHOTOSGELES 


made it possible/ 


On the screen of the Electron Microscope, the scientist 
examines the characteristics of tiny organisms less than a 
millionth of an inch long. Electronics have, in fact, brought 
under scientific examination many germs and bacteria so small 
that they have never previously been observed in detail by 
the human eye. 

G.E.C. Cathode Ray Tubes and OSRAM Valves cover every 
electronic application, and will bring to the pursuits of peace 
many well-tried electronic devices to speed, smooth, and 
make safer our way of life. 


GEC  ogam 


CATHODE RAY TUBES 


of The General Eleciric Co. Ltd., Magnet House, Kingsway, London, W.C.2, 
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Current, A.C. and D.C. (0. to 10 amps.) 

Voltage, A.C. and D.C. (0 to 1000 v.) 

Resistance (up to 40 megohms.) 

Capacity (0 to 20 mfds.) 

Audio-frequency Power Output (0 to 
2 watts.) 

Decibels (—25 Db. to +16 Db.) 


The UNIVERSAL AVOMINOR 


Electrical Measuring Instrument 


An accurate moving coil meter providing 22 ranges 
of readings on a 3-inch scale. Total resistance 
200,000 ohms. Self-contained for resistance 
measurements up to 20,000 ohms, and by using an 
external source of voltage the resistance ranges 
can be extended to 10 megohms. The ohms com- 
pensator for incorrect voltage works on all ranges. 
Suitable for use as an output meter when the A.C. 
voltage ranges are being used. Complete with 
leads, testing prods, crocodile clips, and instruction 
booklet. 


D.C. Voliage : 0—500 volts. 
A.C, Voltage  : 0—500 volts. 
D.C. Current : 0—500 milliamps. 
Resistance : O—10 meghoms. 


Size: 43” x 327 x 12’. Nett Weight: 18 ozs. 


Sole Proprietors and Manufacturers :— 


THE AUTOMATIC COIL WINDER & ELECTRICAL EQUIPMENT CoO., LT 


The MODEL 7 50-Range Universal 


AVOMETER | 


Electrical Measuring Instrument 


The Model 7 Universal AvoMeter is a multi-range A.C./O 
Measuring Instrument of B.S. Ist grade accuracy, provia 
50 ranges of readings on a 5” hand calibrated scale. Ra 
selection is by means of two rotary switches, for A.C. and Oy 
respectively. A press button provides an additional range 
each value of current and voltage shown on the switch kne 
Current consumption is | mA. or 2 mA. at full scale deflec 
A. self-contained 


The MODEL 40 Universal 
AVOMETER 


according to whether the press. button is used or not. 
total resistance of the Meter is 500,000 ohms. 
The Meter is protected by an automatic cut-out aga’ 
damage through severe over-load, and compensation is 
vided for variations in ambient temperature. 
range A.C./D.C. ing 
ment providing 40 ra 
of readings of cur 
voltage and resist} 
without external sh 


Size: 8” x7} x 44”. 

Nett Weight: 6% Ibs. (including leads). 

Various accessories are available for extending the above ran 
transformers or 

pliers. Range selecti 

accomplished by m| 


of A.C. and D.C. sw 
knobs, also a press bull 
which halves the val 
any current or vol 
range. } 
A 3mA. movemeni 
incorporated, full g 
deflection on vo 
ranges being gene 
obtained with the 
sumption of 3 or 6. 
The resistance of 
meter on the 600 
1,200 volt ranges 
200,000 ohms. | 
5” hand calibrated sq 
with parallax mi 
Accuracy within B.S.. 
grade limits guaran) 
on both A.C. and ff 
Provided with auton# 
overload cut-out, 


Size : 8” x 72” x 4h”. 


Weight: 64 Ibs. 
(including leads). 


i 


! 


|) 
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Registered Trade Mark 


The “‘AVO”’ ALL-WAVE 
OSCILLATOR 


An inexpensive, accurate, modulated oscillator covering a 
continuous fundamental frequency band from 95 Ke. to 40 Mc. 
by means of 6 separate coils mounted on a rotary turret 
selector. A harmonic calibration extends the range to 80Mc. 
A large, clearly marked dial is directly calibrated throughout, 
accuracy being within 1%. Externally modulated, internally 
modulated, or R.F. signals obtainable at will. Internal modula- 
tion at 30% at 400 cycles, the L.F. signal being available for 
external use. Calibrated double attenuator enables signal 
to be varied from a few microvolts to 50-millivolts, with a 
force output of Iv. Self-contained, fully shielded. Shielded 
output lead and interchangeable dummy aerials. A.C. Mains 
and Battery Models. 


Size: 10}” x 84” x 44” Nett Weight:134 Ibs 


thew AVO. = VALVE TESTER 


Complete 
with Univer- 
sal valve panel 
fitted with Selector 
Switch enabling all 
English, American and 
Continental valves to be 
easily tested. Compre- 
hensive valve data booklets, 
containin test details o 
A portable self-contained 50-cycles bridge of nearly 9,060 English eat 
exceptional accuracy and utility for direct measure- Krnericda: valvesdics supplied 
ment of all normal values of condensers and ? Se 
resistances. Except at the extreme ends of the 


scale, accuracy is within 5%. ae 
Facilities also provided for 6 F , 
CAPACITY condenser Bowen factor mea- Indicates mutual conductance direct in mA/V. Tests at makers 


i f English 

aaa t d leakage tests specified electrode voltages. One panel for all types o y 

Ne at ek hg By cha flashing Pe had: and American valves. ee cele oeie (aan ae 
: =e { i i ; easily operated multi-roller Universal Selector Switch. Tests 
- a abies a earrontces Cathode/Heater insulation of an indirectly heated valve when 


i dards. the valve is hot; emission of diodes and rectifiers ; indicates 
| SAN scary riage See filament continuity and“electrode “‘ shorts. Coloured scale 
> 


: Ae Age iverascc 1 2lndi ” oR 
500—5 000000hms efficient valve voltmeter indi- indicates state o va' e as ke) ie ndifferent or = C= 
50. 060 ohms to cator for measurement of place.’’ No calculations of any kind involved. Works entirely 

2 


50 Megohms. both audio and radio fre- from A.C. mains. 
quency voltages. 


9ER HOUSE, DOUGLAS STREET, LONDON, S.W.|I Telephone: ViCtoria 3404-8 


casa scccss eee EEE 


A Standardof. Zero Loss Angle 


VARIABLE AIR CONDENSER TYPE D-14-A 


This three-terminal double-screened condenser is provided with a guard 
circuit which ensures that the dielectric of the plate-to-plate capacitance is 
composed entirely of air. This, together with the special surface treatment 
of the plates reduces the plate-to-plate power loss to a quantity which 
can be disregarded even when measuring the smallest power factors. 


BRIEF SPECIFICATION 


CAPACITANCE. 50 upF min. 1,250 »uF max. SCALE READING. To 1 part in 5,000 direct reading ; 
LOSS ANGLE. Approximately 1 micro-radian in a dry To 1 part in 20,000 by interpolation. 
atmosphere: 7 micro-radians in 75% relative humidity, BACKLASH. Not exceeding 1 part in 20,000. 


for the frequency range 50 c.p.s. to 10,000 c.p.s. DIMENSIONS. 12. 7/8” x 10’’ x 13 5/8”. 


DRIVE. Worm reduction gear, 50: 1 ratio. Write for Bulletin B-537-A giving further particulars, 


MUIRHEAD & COMPANY LTD., ELMERS END, BECKENHAM, KENT. TEL: BECKENHAM (041-2. 


] 
FOR OVER 60 YEARS DESIGNERS AND MAKERS OF PRECISION IN STRUMEN 
~HEO ne ae Be Ae ar Pe Te SEER SENET aE E STE ets : - ; i 


O.R.C. 45 
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RESEARCH & 
ENGINEERING 


‘B TH Products are available 
as instruments and _ apparatus 
for research, development, and 
engineering, including all branches 
of electronic control. 


THYRATRONS MAGNETRONS 
CRYSTALS REGHIEIERS 
SPECIAL VACUUM ELECTRONIC DEVICES 
GLASS, AND GLASS METAL SEALS 
RECTIFIERS 
FRACTIONAL HORSEPOWER MOTORS 
ELC: 


LIMITED, RUGBY, ENGLAND. 


A3593N 
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RAY TUBE 


WV henever the Oscillograph trace requires recording ILFORD LIMITEL 
can supply the MOST SUITABLE photographic material. Ilford Limited marke 
a full range of papers, films and plates for photographic recording of all YP 
used in scientific and industrial investigations. The following is a selection af 
materials specially made for OSCILLOGRAPH RECORDING : 


CODE NATURE CHIEF USES . } 
a i] 
| 


5R101 Film—Panchromatic. Specially suitable for recording traces on red fluoresce 
Very high speed. | screens (Phosphate screens) and sometimes used foi 
Medium contrast. | green and blue screens. | 


F.P.3 Plate—Panchromatic. Sometimes used for recording single transients as it if 
Medium contrast. | slightly faster than 5RIOI. 


5G9| Film—Orthochromatic. | The fastest material for recording traces on green anif 
Very high speed. | blue fluorescent screens. Specially suitable for record 
Medium contrast. | ing single fast transients. | 


i sea ae rai mc 
5B52 Film—-Blue sensitive. A new film specially for recording traces on blu 
High speed. fluorescent screens. On account of its special methoff 
Medium contrast. | of manufacture there is a minimum of image sprealf 
On orange base. | enabling traces, in which the spot velocity varies ove 
a wide range, to be successfully recorded. Also suitablif 
for high-speed processing. 


1 


BPI Paper—Glossy. Very widely used for recording traces on blue fluo#f 
Blue sensitive. escent screens and also for many other forms of recoraf 
High speed. ing where a paper base is required. | 
High contrast. 

8B2| Film—Blue sensitive. Suitable for the direct recording of high-speed electro | 
Medium speed. in an evacuated tube. | 


High contrast. i 


For an explanation of the code numbers used and for fuller details of these and other terormmenal 


CF 


materials, write for a copy of the Ilford booklet ‘‘ Photography as an Aid to Scientific Work. | 


ILFORD _LIMITED + ILFORD - LONDON 
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THE KENT 


MULTELEC 


is a strong. and sturdy instrument which, upon activation at two- 
second intervals by any minute change of potential, sets up a 
powerful servo-motor force for controlling any desired form of 
correction, yet without imposing any restraint upon the sensitivity 
of the instrument. Furthermore, it indicates and records. 


VARIOUS USES : 

The Kent Multelec is extensively used in various industries for 
recording, indicating and controlling temperatures, its potentiometric 
basis ensures a fine degree of accuracy under industrial conditions. 


It is also used in the control of pH value in various processes, and 
for the correction of boiler feed water by measurement of 
electrical conductivity. 


It will be recognised that the Multelec can be adapted to many 
forms of controlling, recording or indicating besides the above. 


George Kent, Ltd., also supply meters for accurate measurement 
of the flow of steam, water, gas, air, oil or any other fluid, 
whatever the quantity or pressure. 


The Multelec booklet (No. 897) is available on request. 


GEORGE KENT, LTD. LUTON & LONDON 


XE 
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GLASS TO METAL SEALS 


* 


GLASS TO METAL SEALS 


x 


GLASS TO METAL SEALS 


ah 


OSEL 


(LOW EXPANSION STEEL ) 


SPECIAL ALLOY STEEL 


’ 


‘“« Nicosel ’’ is a valuable British contribution to the Arts of electric 
lamp and vacuum tube manufacture. Its low co-efficient of expansion 
makes it eminently suitable for modern developments of the glass and 


electrical industries. 


© 
Ate 


SPELL VV EIN) 


THOS. FIRTH € JOHN BROWN LID = SHEFFIELD 
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RESISTANCES 


Although present circumstances render 
it difficult for us to give our pre-war 
service to all customers, we are still work- 
ing in their interests. 


New materials and manufacturing pro- 
cesses, which we are now using to increase 
output, also contribute in large measure to 
improved performance and reliability of 
our products. Thus, when normal times 
return, all users of Berco Resistances 
will benefit by our work to-day. 


THE BRITISH ELECTRIC RESISTANCE CO. LTD. 
QUEENSWAY, PONDERS END, MIODLESEX 


Telephone: HOWARD 1492. 
Telegrams: ‘‘ VITROHM, ENFIELD.” 


a nite 


or HEAVY DUTY 


than any other non-ferrous alloy, a higher conductivity than 

any of the bronzes and excellent resistance to corrosion and 

wear—these characteristics of Mallory 73 Beryllium Copper 

have made it first choice for instrument springs, diaphragms 

and bellows, current-carrying springs, snap action switch Properties of 

blades, contact blades and clips. MALLORY 73 BERYLLIUM COPPER 
after heat treatment 


Greater tensile, elastic and fatigue strengths BERYLLI UM COPPER . 


Supplied annealed or lightly cold worked, it has good 
forming properties and is readily fabricated into springs 


and parts of complicated shape. A simple heat treatment dace y 5 Se no 
then develops its remarkable properties. ar Nee: : 
Limit of proportionality 


I 
Available as sheet, strip and wire, in a range of tempers to tons per sq. inch 47-50 | 
suit users’ requirements, and as rod, tube, precision rolled Fatigue limit | 
hair-spring strip, and silver-faced contact bi-metal strip. tons per sq. inch + 19-20 
é A Vickers pyramid hardness 
: 2 350-420 
Electrical Conductivity 
MALLORY METALLURGICAL PRODUCTS LTD. per cent 1.A.C.S. 23.25 
An Associate Company of JOHNSON, MATTHEY & CO. LTD. | 
78 Hatton Garden, London, E.C.I Full details are given in our beoklet, 
Telephone: HOLborn 5027 which will be sent on request. 
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CAMBRIDG | 
POLAROGRAPH 


(HEYROVSKY DROPPING MERCURY ELECTRODE METHOD} 


Some Fields of Application : 


PHYSICAL CHEMISTRY 
METALLURGY 
MINERALOGY 

SUGAR 

EXPLOSIVES 

GASES 

AGRICULTURE 
MEDICINE 

OIL DISTILLING 
DYEING & BLEACHING | 
CERAMICS 

WATER PURIFICATION | 


at “HIS 36-page publication (109-L) will be sent to responsible 
executives and research workers on request. It contains 
much valuable information on the technique of polaro- 
graphic analysis, with formule for base solutions, tables and 
charts of half-wave potentials, and a bibliography of 168 impor- 
tant references. A special feature is the reproduction of actual I 
records obtained in original experiments in our own laboratory. 
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ABSTRACT. Six-figure tables of the trichromatic coefficients of the colours of 
Planckian (total or black-body) radiators expressed on the C.I.E. trichromatic system 
(1931) are given for the temperature range 0-1-660 mireds. The colours refer to energy 
distributions calculated from Planck’s formula with a value 14384-8 micron degrees for 
the constant C,. Additional tables give the differences in the coefficients of the colours 
from the tabulated values for a change in the value of C, from 14384:8 to 14320 and from 
14384-8 to 14350.. These additional tables can be used to calculate the colours if other 
values of C, are used in Planck’s formula. 


radiators expressed on the C.I.E. system (1931) for the temperature 

range 1500-10-9000° k. have been published (Harding, 1944). These 
six-figure tables enable the coefficients of the colours of total radiators at any 
temperature within the given range to be calculated to five-figure accuracy by 
linear interpolation. To extend the range of these tables to cover the tem- 
perature range between 9000° k. and infinite temperature, the colours have been 
recalculated on a mired temperature scale and the range extended to infinite 
temperature. A mired, or micro-reciprocal-degree, is a million divided by the 
temperature of the radiator expressed in °K. (Priest, 1933). 

The colours of the total radiators for the range 0-1—-660 mireds were obtained 
by calculating the coefficients of the colours first for 0-10-660 mireds, then 
interpolating these coefficients to tenths. 

The colours of the total radiators for 0-10-110 mireds had not been published 
and so they were calculated in the following way. ‘The spectral distributions of 
the energies of the radiators were calculated from Planck’s formula, 


=f 
LeydA= Ci | exp a = i] dr, . 


where EF, dA represents the amount of energy radiated at MM mireds between the 


Ele ABLES which give the trichromatic coefficients of the colours of total 


wave-lengths A+ a (microns) and C, and C, are constants. Since the energy 


values were made equal to 100-0000 at 0:56, the value of C, was not required. 
The value 14384-8 micron degrees was given to C, (Birge, 1941). The distribu- 
tions of energy were checked by fifth differences and were estimated to be correct 
to about one part ina million. The trichromatic coefficients of the colours of the 
radiators were then calculated from the spectral distributions of energy, using the 
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Condensed Tables for Colour Computation for the equal-energy stimulus ( Smith, 
1934). 

The colours of the radiators at 120-10-660 mireds were calculated from 
the previously published tables and the trichromatic coefficients were altered in 
some instances by a unit in the sixth decimal place to make the fifth differences 
smooth. 

The trichromatic coefficients obtained by the methods just explained for 
the range 0-10-660 mireds were then interpolated to tenths, giving values at 
0-1-660 mireds. Alterations of a unit in the sixth decimal place were made to 
the Z coefficients where necessary to make the sum of all three coefficients exactly 
unity. These coefficients are given in table 1. 

The alterations to be made to the trichromatic coefficients if C, is changed 
from 14384:8 to 14320 were calculated for the mired range 0-10—660 mireds in 
the following way. Two illuminants have the same spectral distribution of energy, ~ 
consequently the same colour, if the value of C,M in Planck’s formula remains 
unaltered. It follows that the trichromatic coefficients of the colour of a total 
radiator at M,4329 mireds, calculated with C,=14320, willbe the same as those of a 
radiator at 14320 (14384-8)- 1My4sg4-3=0-9954952 M4394. mireds if 14384-8 has 
been used for the colour calculation. The differences between the trichromatic 
coefficients for Mj4354-3 maireds and 0-9954952 M,43,., mireds were therefore 
calculated from the trichromatic coefficients in table 1 and are tabulated in table 2. 
If the value of C, is 14350, the alterations to the coefficients in table 1 are given 
by the differences between the coefficients for Mj4354.3 mireds and 14350 (14384-8)4 
Moyazga-g = 9°9975808 My a5g4- mireds. ‘These differences are tabulated in table 3. 
All the differences in tables 2 and 3 were differenced to the third order and an 
alteration of a unit in the sixth decimal place was made to a few values to make the 
third differences smooth. 

For other likely values of C, different from 14320 or 14350, linear interpolation 
of the values given either in table 2 or table 3 gives the differences with errors 
less than one or two units in the sixth decimal place; for instance, if C, = 14360 
the X coefficient corresponding to 150 mireds is equal to 

14384-8 — 14360 


0-310998 — 0-000237 4384-8 — 14350 = 0-310998 — 0-000169 =0-31083. 


The six figures for the trichromatic coefficients are only quoted here in order that 
the smoothed tables might serve as a basis for the preparation of other tables 
that might be required with a value of C, other than 14384:8. 

The values quoted in tables 1, 2 and 3 might in a few instances be one or two 
units in the sixth decimal place different from those in the previously published 
tables. These differences were considered to be necessary to smooth the values 
in the present tables. 
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Colours of total radiators expressed on the C.L.E. system 


‘Table 1 


Mireds 
(CHa are le 
(Cy=14384°8) 


‘Trichromatic coefficients 


exo lO’ i alOS 2x108 
239906 Pee AUSyI| 525963 
240163 ero} 525324 
240423 234898 524679 
240687 235287 524026 
240953 235680 523367 
241221 236077 522702 
241493 236477 522030 
241768 236882 521350 
242046 237290 520664 
242327 237702 SHA 
242610 238118 519272 
242897 238538 518565 
243187 238961 517852 
DE SMD 239388 BW 6 
Mesos 239819 516406 
244074 240254 515672 
244376 240692 514932 
244680 241134 514186 
244988 241580 513432 
245299 242030 512671 
245613 242483 511904 
245930 242940 511130 
246250 243401 510349 
246574 243865 509561 
246900 244333 508767 
DAT D9 244805 507966 
247562 245280 507158 
247897 245759 506344 
248236 246241 505523 
248578 246727 504695 
248923 247216 503861 
249271 247709 503020 
249622 248206 502172 
249977 248705 501318 
250334 249209 500457 
250695 249715 499590 
251058 250225 498717 
291425 250739 497836 
251795 IWS 496950 
252168 201775 496057 
Dao 45 252299 495156 
252924 252825 494251 
253307 253355 493338 
253692 253888 492420 
254081 254424 491495 


N 


Mireds 
CIP 2 NEN! 
(C,=14384-8) 


a 
at ee SS ee 
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Table 1 (continued) 


Trichromatic coefficients 


| y 108 ax 108 
254473 254963 490564 
254868 255505 489627 
255266 256050 488684 
255668 256598 487734 
256072 257150 486778 
256480 257704 485816 
256890 258261 484849 
257304 258821 483875 
257721 259384 482895 
258141 259950 481909 
258564 260518 480918 
258990 261090 479920 
259419 261664 478917 
259851 262240 477909 
260287 262820 476893 
260725 263402 475873 
261166 263986 474848 
261611 264573 473816 
262058 265162 472780 
262509 265754 471737 
262962 266349 470689 
263419 266945 469636 
263879 267544 468577 
264341 268145 467514 
264807 268749 4660444 
265275 269355 465370 
265747 269962 464291 
266221 270572 463207 
266698 271184 462118 
267178 271798 461024 
267662 272414 459924 
268148 273032 458820 
268636 273652 457712 
269128 274274 456598 
269623 274898 455479 
270120 ZISS28 454357 
270621 276150 453229 
271124 276779 452097 
271630 277409 450961 
272138 278041 449821 
272650 278675 448675 
273164 279310 447526 
273681 279946 446373 
274201 280584 445215 
274723 281224 444053 


Colours of total radiators expressed on the C.I.E. system 
Table 1 (continued) 


Wureds Trichromatic coefficients 
Cire )ia 1108 
8) x x 108 | y x 108 | zx 108 
eee ree eer ae ee eee tle fees he ee ee 

90 275249 281864 442887 
94 275776 282506 441718 
92 276307 283150 440543 
93 276840 283794 439366 
94 277376 284440 438184 
95 277914 285086 437000 
96 278455 285734 435811 
97 278999 286383 434618 
98 279545 287033 433422 
99 280093 287684 432223 
100 280644 288336 431020 
101 281197 288989 429814 
102 281754 289642 428604 
103 282313 290296 427391 
104 282874 290951 426175 
105 283437 291607 424956 
106 284003 292263 423734 
107 284572 292920 422508 
108 285143 293578 421279 
109 285716 294236 420048 
110 286291 294894 418815 
tla 286869 295553 417578 
112 287449 296212 416339 
113 288031 296872 415097 
114 288615 297533 413852 
aS 289202 298193 412605 
116 289791 298853 411356 
117 290382 299514 410104 
118 290976 300175 408849 
119 291571 300836 407593 
120 292169 301497 406334 
121 292769 302158 405073 
2 293371 302819 403810 
123 293975 303480 402545 
124 294581 304141 401278 
125 295190 304802 400008 
126 295800 305463 . 398737 
Wy) 296412 306123 397465 
128 297027 306783 396190 
129 297643 ; 307443 394914 
130 298261 308103 393636 
131 ‘ 298881 308762 392357 
132 299503 309421 391076 
133 300127 310080 389793 


134 300752 310738 388510 
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Table 1 (continued) 


Mireds Trichromatic coefficients | 
(T°K.)+ x 108 : i 
(Cy=14384°8) xx 108 | ven Oe | BS NOe 

135 301379 311396 387225 
136 302009 312053 385938 
137 302640 312709 384651 
138 303273 313365 383362 
139 303908 314020 382072 
140 304544 314675 380781 
141 305182 315329 379489 
142 305822 315982 378196 
143 306463 - 316634 376903 
144 307106 317286 375608 
145 307751 317937 374312 
146 308397 318587 373016 
147 309045 319236 371719 
148 309695 319884 370421 
149 310346 320531 369123 
150 310998 321177 367825 
151 311652 321822 366526 
152 312307 322466 365227 
153 312964 323109 363927 
154 313623 323751 362626 
155 314282 324392 - 361326 
156 314943 325031 360026 
157 315606 325669 358725 
158 316270 326307 357423 
159 316935 326942 356123 
160 317601 397577 354822 
161 318269 328210 353521 
162 318937 328842 352221 
163 319608 329473 350919 
164 320279 330102 349619 
165 320951 330730 348319 
166 321625 331356 347019 
167 322300 331980 345720 
168 322976 332604 344420 
169 323653 333226 343121 
170 324331 333846 341823 
171 325010 334465 340525 
172 325691 335082, 339227 A 
173 326372 335697 337931 
174 327054 336311 336635 
175 327738 336923 335339 
176 328422 337534 * 334044 
177 329108 338142 332750 
178 329794 338749 331457 


ey) 330481 339359 330164 
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Table 1 (continued) 


Mireds Trichromatic coefficients 
CE.) S108 
a, | | 
(Gm 14384 8) x > 106 | y & 108 | 2x 108 
a 
180 331169 339958 328873 
181 331858 340560 327582 
182 332548 SHAS 326293 
183 333238 Solwion 325005 
184 BSS I29 342353 323718 
185 334621 342947 322432 
186 335314 343539 321147 
187 336008 344129 319863 
188 336702 3447 17, 318581 
189 337397 345304 317299 
190 338093 345888 316019 
191 338789 346470 314741 
1 339486 347050 313464 
193 340184 347629 312187 
194 340882 348205 310913 
195 341581 348779 309640 
196 342281 349351 308368 
197 342981 349920 307099 
198 343681 350488 305831 
199 344382 351054 304564 
200 345084 351617 303299 
201 345786 352178 302036 
202 346489 SOS 300774 
203 SOM BS2o4. ZOOS. 
204 347895 353849 298256 
205 348599 354401 297000 
206 349303 354951 295746 
207 350007 355499 294494 
208 350712 356045 DOS Ta 
209 351418 356588 291994 
210 BO 123 357129 290748 
PAL 352829 357668 289503 
212 353535 358204 288261 
PAN 354241 358738 287021 
214 354947 359270 285783 
INS: 355654 359799 284547 
216 356361 360326 283313 
217 357068 369851 282081 
218 357775 361373 280852 
219 358482 361893 279625 
220 359190 362410 278400 
227A 359898 362925 : ZITAT 
PypIpL 360606 363437 275957 
Ds 361314 363947 274739 


224 362022 364454 | 273524 
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Table 1 (continued) 


Trichromatic coefficients 


Mireds 

(FeK) =< 108 

(C,=14384:-8) xX 108 y x 108 ex 108 
225 362730 364959 DIP 
226 363438 365461 DAMON 
227 364146 365961 269893 
228 364854 366458 268688 
229 365563 366953 267484 
230 366271 367445 266284 
Dei 366979 367935 265086 
232 367688 368422 263890 
233 368396 368906 262698 
234 369104 369389 261507 
255 369812 369868 260320 
236 370520 370345 259135 
237 371228 370819 257953 
238 371935 371291 256774 
239 372643 371760 255597 
240 373350 372227 254423 
241 374057 372691 253252 
242 374764 373i 252084 
243 375470 373611 250919 
244 SH Oni 374067 249756 
245 376883 374521 248596 
246 377589 374972 247439 
247 378294 375420 246286 
248 378999 375866 245135 
249 379704 376309 243987 
250 380409 376749 242842 
Dail 381113 377187 241700 
5) 381817 377622 240561 
253 382521 378054 239425 
254 383224 378483 238293 
255 383927 378910 237163 
256 384630 379335 236035 
Di 385332 379756 234912 
258 386034 380175 233791 
259 386735 380591 232674 
260 387436 381005 231559 
261 388136 381416 230448 
262 388836 381824 229340 
263 389535 382230 228235 
264 390234 382632 227134 
265 390933 383033 226034 
266 391630 383430 224930 
267 392328 383825 223847 
268 393024 384217 222759 
269 393720 384606 221674 


Colours of total radiators expressed on the C.L.E. system 


Table 1 (continued) 


Mireds 
Ci) = 1108 
(Cy=14384-8) 


x x 108 


394416 
SEIU 
395805 
396499 
397193 


397885 
398577 
399269 
399959 
400650 


401339 
402028 
402716 
403403 
404090 


404776 
405461 
406146 
406830 
407513 


408195 
408876 
409557 
410237 
410915 


411593 
412271 
412947 
413623 
414297 


Oa 
415644 
416316 
416987 
417657 


418327 
418995 
419663 
420330 
420995 


421660 
422324 
422987 
423648 
424309 


‘Trichromatic coefficients 


| y 108 | 


384993 
385377 
385758 
386137 
386513 


386886 
387257 
387625 
387990 
388352 


388712 
389069 
389423 
389775 
390124 


390470 
390814 
391154 
391493 
391828 


392161 
392491 
392819 
393144 
393466 


393786 
394103 
S944 
\ 394729 
395038 


Sess 
395648 
395949 
396248 
396543 


396837 
SAT 
397415 
397701 
397984 


398264 
398542 
398816 
399089 
399358 


eS Ke 


220591 
219512 
218437 
217364 
216294 


215229 
214166 
213106 
212051 
210998 


209949 
208903 
207861 
206822 
205786 


204754 
203725 
202700 
201677 
200659 


199644 
198633 
197624 
196619 
195619 


194621 
193626 
192636 
191648 
190665 


189685 
188708 
187735 
186765 
185800 


184836 
183878 
182922 
181969 
181021 


180076 
179134 
178197 
177263 
176333 


9 


Io 
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Table 1 (continued) 


Trichromatic coefficients 


Mireds 
(T °K.) 108 
(Cy= 14384°8) xx 108 y x 108 
315 424969 399626 
316 425628 399890 
317 426286 400152 
318 426943 400412 
319 427599 400669 
320 428254 400923 
Bt 428908 401175 
322 429561 401424 
BY8 430213 401671 
324 430863 401915 
B25 431513 402157 
326 432162 402397 
B27. 432809 402633 
328 433456 402868 
329 434102 403100 
330 434746 403329 
331 435389 403556 
332 436032 403780 
333 436673 404002 
334 437313 404222 
SoD 437952 404439 
336 438590 404654 
337 439226 404866 
338 439862 405076 
339 440497 405283 
| 

340 441130 405488 
341 441762 405691 
342 442394 405891 
343 443024 406089 
344 443653 406284 
345 444280 406477 
346 444907 406668 
347 445533 406856 
348 446157 407042 
349 446780 407226 
350 447402 407407 
351 448023 407586 
352 448642 407763 
353 449261 407937 
354 449878 408109 
355 450494 408279 
356 451109 408446 
357 451722 408611 
358 452335 408774 
359 452946 408935 


SOO! 


175405 
174482 
173562 
172645 
GAS 


170823 
169917 
169015 
168116 
167222 


166330 
165441 
164558 
163676 
162798 


161925 
161055 
160188 
159325 
158465 


157609 
156756 
155908 
155062 
154220 


153382 
152547 _ 
151715 
150887 
150063 


149243 
148425 
147611 
146801 


145994 


145191 
301 
143595 
142802 
142013 


141227 
140445 
139667 
138891 
138119 


Colours of total radiators expressed on the C.I.E. system 


Table 1 (continued) 


Muireds 
Ce) 1108 


Trichromatic coefficients 


(Ca— 143848) xX 108 y x 108 | zx 108 
360 453556 409093 137351 
361 454165 409249 136586 
362 454772 409403 135825 
363 455379 409555 135066 
364 455984 409704 134312 
365 456588 409851 133561 
366 457191 409997 132812 
367 457793 410139 132068 
368 458393 410280 131327 
369 458992 410419 130589 
370 459590 410555 129855 
371 460187 410689 129124 
372 460782 410821 128397 
373 461376 410951 127673 
374 461969 411079 126952 
375 462561 411205 126234 
376 463152 411329 125519 
377 463741 411450 124809 
378 464329 411569 124102 
379 464916 411687 123397 
380 465501 411802 122697 
381 466085 411915 122000 
382 466668 412026 121306 
383 467250 412135 120615 
384 467830 412242 119928 
385 468409 412347 119244 
386 468987 412450 118563 
387 469564 412550 117886 
388 470139 412649 117212 
389 470713 412746 116541 
390 471286 412841 115873 
391 471858 412934 115208 
392 472428 413025 114547 
393 472997 413114 113889 
394 473565 413202 113233 
395 474131 413287 112582 
396 474697 413370 111933 
397 475261 413452 111287 
398 475823 413532 110645 
399 476385 413609 110006 
400 476945 413685 109370 
401 477504 413759 108737 
402 478062 413831 108107 
403 478618 413901 107481 
404. 479173 413969 106858 


If 


I2 


Mireds 
(Gi) e108 
(C,=14384°8) 


me ———$—————— ——— — 


410 


414 


415 


419 


420 


424 


AD 


4a) 


430 


434 


435 


439 


440 


445 


449 
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Table 1 (continued) 


x < 108 


479727 
480280 
480831 
481381 
481930 


482477 
483023 
483568 
484112 
484654 


485195 
485735 
486273 
486810 
487346 


487881 
488414 
488946 
489477 
490007 


490535 
491062 
491588 
AODAT2 
492635 


=e Sy7 
493678 
494197 
494715 
495232 


495748 
496262 
496775 
497287 
497798 


498307 
498815 
ANS 22 
499827 
500332 


500835 
501336 
501837 
502336 
502834 


Trichromatic coefficients 


| y x 108 


414035 
414100 
414163 
414224 
414283 


414340 
414396 
414449 
414501 
414552 


414600 
414647 
414692 
414735 
414777 


414817 
414855 
414892 
414927 
414960 


AIA9O2 
415022 
415050 
415077 
415102 


415126 
415148 
415168 
415187 
415204 


415220 
415234 
415247 
415258 
415267 


ALSZ 75 
415281 
415286 
415290 
415292 


15292 
415291 
415289 
415285 
415280 


2 x 10° 


106238 
105620 
105006 
104395 
103787 


103183 
102581 
101983 
101387 
100794 


100205 
99618 
99035 
98455 
97877 


97302 
96731 
96162 
95596 
95033 


94473 
93916 
93362 
92811 
92263 


DAT 
91174 
.90635- 
90098 
89564 


89032 
88504 
87978 
87455 
86935 


86418 
85904 
85392 
84883 
84376 


83873 
83373 
82874 
82379 
81886 


Colours of total radiators expressed on the C.I.E. system 


Table 1 (continued) 


Mireds 
(ii) 1108 
(C.=14384-8) 


‘Trichromatic coefficients 


503331 
503827 
504321 
504814 
505306 


505796 
506286 
506774 
507261 
507747 


508231 
508714 
509196 
509677 
510156 


510635 
511112 
511588 
512062 
512536 


513008 
S349 
513949 
514418 
514885 


S1S351 
515817 
516280 
516743 
517205 


517665 
518124 
518582 
519039 
519494 


519949 
520402 
520854 
521305 
521754 


522203 
522650 
523097 
523542 
523986 


A273 
415265 
ASD S)5. 
415244 
A sy2 337 


415218 
415203 
415186 
415169 
415149 


An) 
415107 
415084 
415060 
415034 


415007 
414979 
414950 
414919 
414887 


414854 
414820 
414784 
414747 
414709 


414070 
414630 
414588 
414545 
414501 


414456 
414410 
414362 
ANAS IA: 
414264 


414213 
414161 
414108 
414054 
413998 


413942 
413885 
413826 
413767 
413706 


81396 
80908 
80424 
79942 
79462 


78986 
78511 
78040 
77570 
77104 


76640 
76179 
75720 
75263 
74810 


74358 
73909 
73462 
73019 
T2a Td 


72138 
71701 
71267 
70835 
70406 


C2979 
69553 
69132 
68712 
68294 


67879 
67466 
67056 
66647 
66242 


65838 
65437 
65038 
64641 
64248 


63855 
63465 
63077 
62691 
62308 


¥ 
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Table 1 (continued) 


Trichromatic coefficients 


Mireds 

(QP Nae NOE 

(C= 14384-8) xX 108 yx 108 2x 108 
495 524428 4136044 61928 
496 524870 413582 61548 
497 525310 413518 61172 
498 525750 413453 60797 
499 526188 412388 60424 
500 526625 413321 60054 
501 527061 413253 59686 
502 527496 413185 59319 
503 a27929 413115 58956 
504 528362 413044 58594 
505 528793 412973 58234 
506 529223 412900 57877 
507 529652 412826 57522 
508 530080 412752 57168 
509 530507 412676 56817 
510 530933 412600 56467 
511 531358 412523 56119 
52 531781 412444 S75) 
513 532204 412365 55431 
514 532625 412285 55090 
515 533045 © 412204 54751 
516 533464 412122 54414 
Sly 533882 412039 54079 
518 534299 411956 53745 
519 534715 411871 53414 
520 535130 411786 53084 
Sil 535544 411700 52756 
522 535956 411613 52431 
523 536368 411525 52107 
524 536778 411437 51785 
525 537187 411347 51466 
526 537596 411257 51147 
Si 538003 411166 50831 
528 538409 411074 50517 
529 538814 410981 50205 
530 539218 410888 49894 
531 539621 410794 49585 
a2 540023 410699 49278 
588) 540424 410603 48973 
534 540824 410507 48669 
535 541223 410409 48368 
536 541620 410311 48069 
587 542017 410213 47770 
538 542413 410113 47474 
539 542807 410013 47180 


Colours of total radiators expressed on the C.1.E. system 


Table 1 (continued) 


Mireds 
ike ee al(08 
(C.=14384-8) 


‘Trichromatic coefficients 


pore | ay NOP | 2 10° 
540 543201 409912 46887 
541 543594 409810 46596 
542 543985 409708 46307 
543 544376 409605 46019 
544 544765 409501 45734 
545 545154 409397 45449 
546 545541 409292 45167 
547 545928 409186 44886 
548 546313 409080 44607 
549 546698 408973 44329 
550 547081 408865 44054 
551 547463 408757 43780 
Dow 547845 408647 43508 
Soe 548225 408538 AVS 
594, 548604 408427 42969 
555 548983 408316 42701 
556 549360 408205 42435 
557 549737 408092 42171 
558 550112 407980 41908 
559 550487 407866 41647 
560 550860 407752 41388 
561 551233 407637 41130 
562 551604 407522 40874 
563 551975 407406 40619 
564 552344 407290 40366 
565 DOING 407173 40114 
566 553081 407056 39863 
567 553447 406937 39616 
568 553813 406819 39368 
569 554178 406700 39122 
570 554542 406580 38878 
5/1 554905 406460 38635 
SW 555267 406339 38394 
573 555628 406218 38154 
BA 555988 406096 37916 
575 556348 405973 37679 
576 556706 405850 37444 
577 557063 405727 rocco 
578 557420 405603 36977 
579 S577 1S: 405479 36746 
580 558130 405354 36516 
581 558484 405229 36287 
582 558836 405103 36061 
583 559188 404977 35835 
584 559539 404850 35611 


LS 
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Table 1 (continued) 


Miureds 
(IPS BSAC 
(C,=14384-8) 


585 
586 
587 
588 
589 


590 
591 
592 
593 
594 


So5 
596 
597 
598 
SoH) 


600 
601 
602 
603 
604 


605 
606 
607 
608 
609 


610 
611 
612 
613 
614 


615 
616 
617 
618 
619 


620 
621 
622 
623 
624 


625 
626 
627 
628 
629 


Trichromatic coefficicnts 


xox 108 y X 108 2x 108 
559889 404723 35388 
560239 404595 35166 
560587 404467 34946 
560934 404338 34728 
561280 404209 34511 
561626 404080 34294 
561971 403950 34079 
562314 403820 33866 
562657 403689 33654 
562999 403558 33443 
563340 403426 33234 
563680 403294 33026 
564019 403162 32819 
564358 403029 32613 
564695 402896 32409 
565032 402762 32206 
565368 402628 32004 
565703 402494 31803 
566037 402359 31604 
566370 402224 31406 
566703 402088 31209 
567035 401952 31013 
567366 401816 30818 
567696 401679 30625 
568025 401542 30433 
568353 401405 30242 
568681 401267 30052 
569007 401129 29864 
569333 400991 29676 
569658 400852 29490 
569982 400713 29305 
570306 400573 29121 
570628 400433 28939 
570950 400293 28757 
571271 400153 28576 
571591 400012 28397 
571910 399871 28219 
572229 399730 28041 
572547 399588 27865 
572863 399446 27691 
573180 399304 27516 
573495 399161 27344 
573809 399018 DATS 
574123 398875 27002 
574436 398732 26832 


a 


Colours of total radiators expressed on the C.I.E. system 


Table 1 (contined) 


Mireds Trichromatic coefficients 
(T°x.)~1 x 108 
(C.=14384-8) cio y 108 2108 
630 574748 398588 26664 
631 575059 398444 26497 
632 575370 398300 26330 
633 575680 398156 26164 
634 575989 398011 26000 
635 576297 397866 25837 
636 576604 397721 25675 
637 576911 397576 25513 
638 Swaia 397430 25353 
639 SH S2D 397284 25194 
640 577826 397138 25036 
641 578130 396992 24878 
642 578432 396845 24723 
643 578734 396698 24568 
644 579036 396551 24413 
645 579336 396404 24260 
646 579636 396256 24108 
647 579935 396109 23956 
648 580234 395961 23805 
649 580531 395812 23657 
650 580828 395664 23508 
651 581124 - 395515 23361 
652 581420 395367 23213 
653 581714 395218 23068 
654 582008 395068 22924 
655 582301 394919 22780 
656 582594 394769 22637 
657 582886 394620 22494 
658 583177 394470 22353 
659 583467 394319 22214 
660 583757 394169 22074 
Table 2 
Change in trichromatic coefficients if C, 
Mireds is altered from 14384°8 to 14320 


(T °x.)-2<108 
(C,=14384-8) 
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x x 108 


y x 108 


2x 108 


17 


18 


nn  —————————— 


Mireds 


AE Ae) el: 
(Cz=14384°8) 


50 
60 
70 
80 
90 


100 
110 
120 
130 
140 


150 
160 
170 
180 
190 


200 
210 
220 
230 
240 


250 
260 
270 
280 
290 


300 
310 
320 
330 
340 


350 
360 
370 
380 
390 


400 
410 
420 
430 
440 


450 
460 
470 
480 
490 
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Table 2 (continued) 


Change in trichromatic coefficients if C, 
is altered from 14384-8 to 14320. 


213 


248 


401 


y x 108 


—125 
157 
191 
PENS, 
260 


zx 10° 


Se 


+217 
276 
339 
404 
473 


542 
611 
680 
748 
814 


877 
938 
994 
1047 
1095 


1140 
1179 
1214 
1244 
1270 


1290 
1306 
1317 
1324 
1327 


1326 
1321 
1312 
1301 
1286 


1268 
1248 
1226 
1202 
EAS 


1148 
1119 
1090 
1059 
1028 


+996 
964 
933 
901 
869 
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Table 2 (continued) 


Change in trichromatic coefficients if C, 
Muireds is altered from 14384-8 to 14320 
. (PEN EESAOE ; 

Rig near) xx 108 yx 108 =x 108 

500 —986 — 149 + 837 

510 980 174 806 
520 O78 198 7754 

530 966 221 745 

540 959 244 TANS 

550 952 266 686 

560 944 286 658 

570 936 306 630 

580 928 325 603 

590 920 343 Sa 

600 ! 1) 360 552 

610 904 377 S27, 

620 896 392 504 

630 888 407 481 

640 880 421 459 

650 871 434 437 

660 863 447 416 

‘Table 3 
Change in trichromatic coefficients if @ 
Mureds is altered from 14384-8 to 14350 
| AGE SR SG 

(C.=14384-8) x x 108 yX 108 zx 108 

0 — 0 pany tad 

10 Zi 10 17 

20 15 22 317) 

30 25 36 61 

40 36 51 87 

50 49 67 116 

60 64 84 148 

70 79 103 182 

80 96 121 217 

90 114 140 254 

100 133 158 291 

110 153 NA 328 

120 173 192 365 

140 215 222 | 437 
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Table 3 (continued) 


Change in trichromatic coefficients if C2 


Mireds is altered from 14384-8 to 14350 

CEekajaca 108 

(C,=14384-8) xx 108 y x 108 | 2x 10° 
150 —237 +234 +471 
160 258 246 504 
170 279 yapys) 534 
180 DOD 263 562 
190 320 268 588 
200 340 Dae 612 
210 358 275 633 
220 Sy YS 652 
230 394 274 668 
240 411 271 682 
250 426 266 692 
260 441 260 701 
270 454 253 707 
280 467 244 711 
290 479 233 fle 
300 489 222 711 
310 499 210 709 
320 507 197 704 
330 515 183 698 
340 521 169 690. 
350 526 154 680 
360 531 L3oun 669 
370 535 122 657 
380 538 106 644 
390 540 90 630 
400 542 73 615 
410 543 57 600: 
420 544 41 585 
430 544 25 569 
440 543 2. 552 
450 542 a 535 
460 540 22 518 
470 538 37 501 
480 535 52 483 
490 532 66 466. 
500 529 81 448 
510 526 94 432 
520 522 107 415 
530 S19 119 400: 
540 515 131 384 
550 511 143 368 
560 507 154 353 
570 503 165 338. 
580 498 75 323 
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Table 3 (continued) 


Change in trichromatic coefficients if C, 
M:reds | is altered from 14384-8 in 14350 
Clie ke 10S es 
(C.=14384-8) | 


| 
nae IO y x 108 SPOS 


| 
600 _—489 +194 +295 
610 485 202 283 
620 481 211 270 
630 476 219 257 
640 472 226 246 
650 467 233 234 
660 463 240 223 


HYSTERESIS AND EDDY LOSSES IN SINGLE 
CRYSTALS OF AN ALLOY OF IRON AND 
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ABSTRACT. ‘The total energy dissipated, in an alternating magnetic field, in single 
crystals of iron containing 2:1% silicon has been measured calorimetrically. Results are 
given for the field applied in the three crystallographic directions [100], [110], [111]. 
The total losses are separated into the conventional ‘‘ hysteresis’ and “‘ eddy’ com- 
ponents by means of measurements at different frequencies. ‘The “‘ eddy ”’ losses do not 
appear to be dependent on direction ; the “ hysteresis ’’ losses for [100] are about one- 
third of those for the other directions measured. A tentative explanation is suggested. 


§1. INTRODUCTION 


continuous heat production, even when the rate is small. A striking 
application was the measurement of the energy of beta-disintegration 
of Radium E (Ellis and Wooster, 1928). Other examples of its use are the 
measurement of losses in dielectrics at radio frequencies (MacGregor Morris 
and Grisedale, 1939) and the study of the breakdown of dielectrics under 
high voltage (Whitehead and Nethercot, 1935). Bates and his co-workers 
have applied a calorimetric method to the somewhat different problem of the 
energy changes taking place in various parts of a single magnetic cycle (Bates 
and Weston, 1941 ; Bates and Healey, 1943). 
In an investigation of electrical sheet steels undertaken at the Cavendish 
Laboratory for the Electrical Research Association, it was necessary to devise 


T= calorimetric method is convenient for measuring the rate of 
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means for measuring energy losses in small samples of material (~5 grams). 
The normal methods require considerably more material, so the possibilities 
of using a calorimetric method were investigated. It has proved somewhat 
slow for routine use, but results have been obtained on single crystals mag- 
netized in three important crystallographic directions. 

The principle of the method is very simple. A thermocouple is enclosed 
in a pack of two or more sheets, and the increase in temperature caused by 
applying an alternating field for a known time is noted. The loss in ergs 
per gram per cycle is then given by 


Cé 


L= aa 


where -C is the heat capacity of the assembly in ergs per degree, 0 is the tem- 
perature increase, ¢ is the time of application of the field in seconds, M is the 
mass of the iron and f is the frequency. ‘T'wo difficulties prevent the use of 
equation (1) in this simple form; first, the observed temperature increase is 
less than it should be because of heat losses during the time of application 
of the field, and secondly, the heat capacity of the assembly is doubtful, because 
of uncertainties in the heat capacity of Fe-Si alloys and lack of knowledge 
as to how much of the thermocouple leads and the surrounding air is effectively 
part of the assembly. It is therefore necessary to correct the observed tem- 
perature rise, and to calibrate the assembly by means of a small heating coil 
incorporated for the purpose. ‘This calibration has the advantage that the 
thermocouple acts only as a transfer instrument, and it is unnecessary to 
know the factor for converting-change in e.m.f. to change in temperature. 


§2, APPARATUS, SPECIMEN AND MEASUREMENTS 


d 


The “‘ calorimeter ’’ was constructed so as to have as small a heat capacity 
as possible; it consisted merely of six junctions of a copper-constantan thermo- 
element, a 1-ohm constantan heating element, and a coil to measure the mag- 
netic induction in the specimen. The junctions and the heating element 
were placed side by side, making a flat plate about 1 cm. square, and held in 
position with cotton thread and varnish. The coil was wound to a length 
of about 0-4 cm. on a rectangular former about 0-25 by 2-8 cm., varnished to 
prevent it from uncoiling, and then fixed to the thermocouple leads so that 


the junctions and heater were supported in the middle of it. ‘There were - 


thus two pockets, one on each side of the junctions and heater, into which 
the specimen to be investigated could be put. 

The specimen consisted of four single-crystal discs of Fe-Si, cut from a 
single crystal about 7 cm. across. After the discs had been annealed at 1000°c. 
for half an hour, their orientation was found by means of back-reflection Laue 
photographs. It was found that [100] made an angle of 26°+2° with the 
rolling direction, and that the plane of the sheet was within 6° of (011). The 
lattice parameter of a portion of the large crystal was 2-8581 kx. (2-8638 a.), 


corresponding to a silicon content of 2:1%. The average diameter of the- 


discs was 2:42) cm., the measured total thickness 0-15 cm., and the total 
mass 5°32; g. ‘The thickness calculated from the mass, density and diameter 
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would be about 0-15, cm. A pair of the discs, separated by a sheet of thin 
paper to prevent electrical contact, was placed in each of the pockets men- 
tioned above, and secured by tying with fine thread. The whole assembly 
was suspended in a Dewar flask and placed in a solenoid, 

The flux-measuring coil was connected directly to a rectifier voltmeter 
whose resistance was high compared with that of the coil. It is easily shown 
that for an ideal rectifier the maximum flux linking the coil is given by ° 


108 
N= Noa acters a(Z) 
where F is the voltmeter reading, and that therefore 
; 10°F 
47] — \/2nfnA pe 8 Oe eae acy oie (3) 


where J’ is the maximum magnetic intensity uncorrected for the field, n is 
the number of turns in the coil and A is the area of cross-section of the iron. 
The value of 47J’ needs correction for the actual magnetic field; the correction 
was zero within the limits of error for 47J less than about 10000, and amounted 
to about 300 at 477=15000. The apparent fields varied from 150 to 1000. 
The largest error in the corrected values of 47I is uncertainty in what ought to 
be taken as the area of cross-section of the discs, since the measured thickness 
is about 5% greater than the calculated thickness. The measured thickness 
was used, so that the values of 47J may all be low by 5%. In addition, there 
are individual errors due to the calibration of the voltmeter, uncertainty in 
the demagnetizing factor and the like. It is estimated that the values of 
4rI are relatively correct within 3°% and absolutely within 8%. 

Before the assembly is heated, either by an alternating magnetic field or 
the heating coil incorporated, the temperature was nearly constant. During 
heating it rose almost linearly, and after heating fell fairly rapidly. The 
actual temperature rise was therefore somewhat less than that corresponding 
to the heat input, and a correction was made by extrapolating the initial and 
final portions of the temperature-time curve to the mid-point of the heating 
period. This was a simple straight-line extrapolation for the initial portion 
of the curve, but the fall was so steep in the final portion that analytical extra- 
polation had to be used. The fall was found to be exponential well within 
the limits of error, and a graphical method was devised for finding the constant 
of the curve rapidly. From observations of the temperature rise for a measured 
input to the heating element and the temperature rise produced by the mag- 
netic field in a known time, the energy loss per magnetic cycle is easily calculated. 
The time of application of the field was usually 180 seconds. The largest 
errors in the losses so found would seem to be due to the extrapolation of the 
final portion of the time-temperature curve to the mid-point of the heating 
period. The percentage error produced depends somewhat on the temperature 
interval, increasing if it is made unduly small or unduly large. For most of 
the measurements it may be estimated at 3 to 4%. The error in the calibration 
measurements will of course affect all determinations of the energy loss in the 
same way, and so produce a systematic error, estimated at 2 to 3%. Other 
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sources of error in the calorimetric measurements are small compared with 
these two; it is estimated, therefore, that the energy losses are relatively correct 
in general within 5% and absolutely within 8%. 


§3. RESULTS 


Measurements of the total loss per cycle were made with the applied field 
in the plane of the discs and having as nearly as possible the directions [100], 
[110], [111], at values of 47J from about 5000 to 15000, and at three different 


Table 1. Smoothed losses at integral values of 47J 
(Extrapolated values are enclosed in parentheses) 


Loss, ergs per gram per cycle 


Anl 3000 5000 7500 10000 12500 15000 17500 
Direction Frequency 
38 = (ty S91) 456, £040) e640) 
100 50 17. = 246 103 181 280 406 586 
[100] 82 A | 63) “AES 4049 G80) 540) ee 
121 — 81 180 316 496 787 = (1150) 
50 39 O72 185 303 461 700 977 
[110] 81 a = (O8)a 210) 1350) es 78 916 (1480) 
145 57 ~ 145° 3024 510° 809° * 4236 (irs) 
50 31 87 197 358 S87) 741 (980) 
[111] 85 (34) 108 225 405 670 (1020) — 
145 (48) 145 298 527 891 1390 = (2040) 
frequencies. From these measurements values of the total loss at integral 


values of 47J were interpolated. 


It was found convenient to use graphs of 


log (total loss) plotted against log (47/), as such plots are nearly linear, whereas 
plots of total loss against 47J are strongly curved. ‘The interpolated values 


Table 2. ‘“‘ Hysteresis” loss and ‘‘eddy”’ loss at 50~ 
(Extrapolated values are enclosed in parentheses) 


Ergs per gram per cycle 


A4nl 3000 5000 7500 10000 12500 15000 #£17500 
Direction 

Hysteresis [100] 6 20 48 79 120 140 (190) 

loss [110] 29 68 120 200 280 400 (480) 
[111] (22) 5D 140 270 350 390 (420 ?) 

Eddy loss [100] 11 26 55 100 160 270 (400) 

per cycle [110] 10 24 62 110 180 300 (500) 
at 50~ {111] (9) 32 il 90 180 350 (560 ?) 

Mean 10 27, 58 100 170 310 — 


are given in table 1. 
given in the table, but are distinguished by being enclosed in parentheses. 


The total loss per cycle for a given value of 47J is, as is usually found, nearly — 


a linear function of frequency. 


Several typical plots are shown in figure 1. 


A few values involving a slight extrapolation are also 
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The loss extrapolated to zero frequency is conventionally attributed to hysteresis, 
and the linear increase with frequency to eddy currents. The losses separated 
in this way are given in table 2. The relative errors are estimated to be not 


/000 


TOTAL LOSS, ERGS PER GRAM PER CYCLE 


0 50 /00 150 
FREQUENCY, CYCLES PER SECOND 


Figure 1. Extrapolation curves for hysteresis loss. 


more than 10%. It will be noticed that the ‘“ hysteresis’? component varies 
markedly with crystallographic direction, being almost three times as great 
for [110] or [111] as for [100], while the ‘“‘eddy’’ component is nearly inde- 
pendent of direction, at any rate for 47J less than 13000. 


§4. DISCUSSION OF RESULTS 
“ Eddy’? losses 


It seems somewhat doubtful if the “‘eddy’’ component of the losses is due 
entirely to eddy currents. ‘The heating due to this cause can be calculated 
with fair certainty; the theoretical formula for the loss per unit volume per 
second at low frequencies (flux density uniform throughout the sheet) is 


ob? /AB\? 
Se (3) ea) 


where o is the conductivity of the material, 2B/dt is time rate of change of the 
magnetic induction (the bar represents the time average over a complete cycle), 
h is the thickness of the sheet, and ¢ is the speed of light. For a sinusoidal 
variation of B at 50~ this becomes 


yy = 822 Bmax 


5 ergs per gram per cycle, Rene ot.) 
P 
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where p is the resistivity of the sheet and 4 its density, both in c.g.s. units. 
For the sheet used, b=0-040 cm., p~35000, 5~7:7, so the loss at Bax = 10000 
should be about 49 erg/gm. cycle. The observed value is 100 erg/gm. cycle, 
which is almost exactly twice as big. A discrepancy of this order is usually 
found in such comparisons (Brailsford, 1943), and has not yet been satisfac- 
torily explained. In some cases it is due, at least in part, to non-sinusoidal 
variation of B, but it may also be necessary to postulate some magnetic vis- 
cosity effect. In the present experiments B was very nearly sinusoidal, since 
discs have a very large demagnetizing factor and the field in the absence of the 
specimen was sinusoidal.* 

The ‘‘eddy”’ losses are very similar for the three directions, as 1s shown 
in figure 2. Up to 47J~12500 they are identical within the estimated limits 
of error ; at 15000 they differ by scarcely more. The separation of “eddy” 
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Figure 2. Eddy losses for the three directions at 5}0~. 
Up to 47J=12500 the losses are hardly separated on this scale. 


and ‘‘hysteresis”’ losses at 17500 is little better than guess-work, but it confirms 
the impression given by the losses at 15000 that at high inductions the “eddy ”’ 
losses increase in the order [100] <[110]<[111]. It seems doubtful if eddy 
losses in the strict sense would behave in this way; magneto-resistance effects 
should not amount to more than a few tenths per cent. This is perhaps 
confirmatory of the hypothesis that “ eddy” losses include magnetic viscosity 
effects dependent on the frequency as well as eddy losses proper, but Richer 
(1944) has advanced a theory which would suggest that eddy currents ought 
to depend on the direction of magnetization. 


‘““ Hysteresis’ losses 


As has been seen, the losses in the [100] direction are much smaller than 
those in the [110] and [111] directions, in fact only about one-third as big. 
There is no great difference between the [110] and [111] directions ; [111] 


tends to be higher in the centre of the range, [110] at the ends. The following 


* The effect of a large demagnetizing factor is to reduce the actual field to a very small value, 


in such a way as to make the magnetic induction in the specimen very nearly proportional to the 
value of the field in the absence of the specimen. 
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tentative discussion of the relation between the losses suggests itself. The 
argument is similar to that of Richer for eddy currents. 

The first step in magnetization along any non-[100] direction is to increase 
the number of domains whose [100] magnetization makes a small angle with 
the field direction, and decrease the number of domains whose [100] mag- 
netization makes a large angle with the field direction. Now, all [100] directions 
are symmetrically placed with respect to [111], and one may therefore suppose 
that all [100] directions behave similarly, and, indeed, as they would behave if 
the field were in a [100] direction, except that it will be the component of the 
field in the [100] direction that is important. An externally measurable 
value of J in a [111] direction will then correspond to an actual internal value 
of J in the [100] directions that is \/3 times as big. The curve of hysteresis 
loss against 47J might therefore be expected to coincide with the similar curve 
for [100] if the values of 4zJ are multiplied by 1/3. The postulated mechanism 
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Figure 3. Hysteresis losses for the three directions. 


for magnetization in the [111] direction is almost certainly that predominant 
for values of the externally-measured 47J less than about 1/1/3 of the saturation 
value, so that the relation between the [111] and the [100] losses should hold 
over the entire range for which the latter have been measured. 

The behaviour of the [110] direction is more difficult to predict. The 
[100] directions are not symmetrically placed with respect to [110] : two- 
thirds of them make angles of 45° (or 135°), one-third of them make angles 
of 90°. The assumption that immediately suggests itself is that the domains 
making angles of 90° with [110] make a 90° shift into one of the directions making 
angles of 45°. One might therefore expect that the curve of hysteresis loss 
againt 477 would coincide with the similar curve for [100] if the values of 
4nI are multiplied by \/2. As in the case of [111], the relation should hold 


over the measured range of [100] losses. 
In figure 3 the measured hysteresis losses are plotted against the externally 
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measured value of 47/ for the three directions,and in figure 4 they are re-plotted 
with the co-ordinates treated as the foregoing argument suggests: [100] 
plotted without change (crosses), [110] plotted with the polarizations multiplied 
by 1/2 (ellipses with long axes vertical), [111] plotted with the polarizations 
multiplied by 1/3 (triangles). The vertical lines are of a length corresponding 
to errors of 10%. It will be seen that up to 477 =15000 the curves for [111] 
and [100] agree within this amount. 

The agreement is quite impressive, and may perhaps be taken as some 
confirmation of the not quite rigorous theoretical argument. On the other 
hand [110] shows no agreement at all, and this probably indicates that the 
assumption made about the behaviour of the domains making angles of 90° 
with the [110] direction is erroneous. Qualitatively one may say that these 
domains are less effective in increasing the polarization than if they behaved 
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Figure 4. Comparable hysteresis losses for the three directions. 
For explanation of the points see text. 


as assumed. By what is probably no more than a curious coincidence, the 
assumption of the opposite extreme (that these domains are entirely ineffective) 
gives good quantitative agreement. On this assumption one would have 
two-thirds of the iron magnetized to an intensity of (3/2)./2 times the externally 
measured intensity, and the curve of hysteresis loss against 47J should coincide 
with that for [100] if the losses are multiplied by 3/2 and the values of 47J by 
3/\/2. ‘The [110] losses treated in this way are also plotted in the figure 
(ellipses with long axes horizontal), and it will be seen that they agree well 
with those for [100] and [111]. 

This explanation of the apparently anomalous behaviour of the [110] 
losses is confirmed by the tendency, postulated by Becker and Doring (1939), 
for the domains whose directions of magnetization lie nearly in the plane of the 


sheet to be favoured at the expense of those whose directions of magnetization 


make large angles with the plane of the sheet. In the present case domains 
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effective for the [100] measurements have their directions in the plane of the 
sheet, whereas those effective for the [110] measurements have their directions 
of magnetization at 45° to the plane of the sheet. This explanation might be 
regarded as qualitatively satisfactory, but the fraction one-third would be the 
result of coincidence, A little consideration shows that the [111] losses would 
not be affected if the greater part of, or even all, the domains had their direction 
of magnetization in the plane of the sheet. Richer (private communication), 
however, has suggested that the apparently anomalous behaviour of the 
[110] losses might be explained by a non-isotropic distribution of interstitial 
impurities set up in the process of rolling the sheet. 
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ABSTRACT. This paper shows how electron-optical systems can be completely 
calculated, without appeal to model experiments, by methods capable of great accuracy. 

The computation method is a numerical one. It is explained with reference to the 
example of an electrostatic lens consisting of two cylinders with diameter ratio 1°5. The 
electric potential and the field are found by relaxation methods. 

Electron paths for paraxial and non-paraxial electrons are computed by means of a 
step-by-step integration method which allows a close check on accuracy. ‘The position 
of focal and cardinal points and the spherical aberration of the lens are found to be in 
fair agreement with experimental and semi-empirical determinations by other authors. 


§1. INTRODUCTION 


HE theoretical treatment of any electron-optical system proceeds in 

two steps: calculation of the electric or magnetic field of the focusing 

system and determination of the path of charged particles in these 

fields. In some cases analytic expressions for the fields are known, but in the 

majority of arrangements in practical use this is not the case. It is then usual 

to determine the field by experiment, either with the help of an electrolytic 
tank or by other means. 

When analytic expressions for the field are known, the second step, path 
determination, can sometimes be carried out analytically. Again, in most cases 
this is not possible, and approximation methods have to be used. In a recent 
paper L. S. Goddard (1944) has given a summary of path-tracing methods. 
In that paper a step-by-step method of path computation is proposed which is used 
by Klemperer and Goddard (1944) for treatment of a magnetic lens. L. J. 
Comrie pointed out in the discussion that other step-by-step methods may be 
preferable. Goddard also indicates that the relaxation method may be used for 
the first step, the field computation. 

We have been using both the relaxation method and a step-by-step method 
proposed by L. J. Comrie during the last few years for electron-optical calcula- 
tions, and it will be shown in this paper that a complete numerical computation 
of electron-optical systems, without appeal to experiment, is practicable. 


An electrostatic lens will be taken as an example, and the agreement of the 


result of the determination of focal length with the results of other authors will be 
discussed. 

The methods used are, mutatis mutandis, equally applicable to axially sym- 
metrical magnetic focusing systems. Even saturation effects in ironclad coils 
can be taken into account. A discussion of magnetic field computation will 
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shortly appear in the Journal of the Institution of Electrical Engineers, where 
the application of the relaxation method to potential problems will be fully 
explained. 


§2. THE FINITE-DIFFERENCE METHOD OF FIELD CALCULATION 


For the purpose of the field determination, Laplace’s equation, 


mS 1 ad OH _ 

C1 ay OTS eta ss 
for the potential ¢ has to be solved with certain boundary conditions. For this 
equation a system of difference equations for the values of ¢ at points of a net 
of small but finite mesh length is substituted. The choice of the net depends 
on the shape of the boundary. In the present paper a study is made of an electro- 
static lens, consisting of two semi-infinite cylinders placed end to end and separated 
by a finite gap. The cylinders as shown in figure 1 have diameters D, and DD; 
with D,/D, =1-5, and their potentials are 200 v. and 1200 v. respectively. This 
choice fixes the boundary conditions. The voltage ratio is thus V,/V,=6, but 
when a solution is found for this ratio, a solution for any other voltage ratio can 
be obtained quite simply (see Maloff and Epstein, 1938, p. 102). Three quadratic 
nets of successively finer mesh length are shown in figure 3. Consider at first 
the coarsest net, of mesh-length a. At any point P of this net the potential ¢ 
can be expressed in terms of the potentials da, dr, ds, Pr of the four neighbouring 
mesh points, Q, R, S, T, figure 2. It can be seen, with the help of an 
expansion in Taylor series of the potential at these points, that 


$a(1- 5) + 4s(14 Fy) tbat bn 4b =0, veiae (2) 


where p=7/a and pp =rp/a, so that rp is the radial distance of point P from the 
axis. If there are m mesh points, excluding points on the boundary where 
¢=const., 2 simultaneous equations have to be satisfied. Owing to the symmetry 
of the figure, the equations on both sides of the axis are the same, so that only one- 
half of figure 1 is used for the computation. 

A solution of these equations can be found by various methods, e.g. the 
Gauss-Doolittle elimination method or R. V. Southwell’s relaxation method 
(Christopherson and Southwell, 1938). We have used the latter method, which 
we prefer because it is self-checking-in so far as mistakes made at some stage of 
the calculation are corrected later with little loss of time and effort. At any 
stage of the process, the residuals 


$a(1- 5.) +4s(1+ 7) +¢pt+dr—tpp =dp...... (3) 


are found and reduced in successive steps until they are negligible. A further 
description of the process is given in a forthcoming paper by Motz and Worthy 
in the Journal of the Institution of Electrical Engineers. | 

At the end a further test should be applied. The system of finite-difference 
equations (2) is accurate only to the second order in a, and we are only justified 
in using it when ¢ does not change rapidly in the neighbourhood of any 


32 H. Motz and Laura Klanfer 


mesh-point. It is possible to locate points or regions where second-order approxi- 
mation is not good enough by forming the higher differences of 4 in rows and 


Figure 1. 


columns. Failure of the second- 


order approximation is disclosed by erratic | 
behaviour of the higher differences, 


Consider, for example, the effect of a change 
* 


of ¢ at one mesh point by +1 on the higher differences. 
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If the fourth differences 


fluctuate by an amount larger than the one which would be produced by errors of 
+1 in the last significant figure of the result, it must be concluded that this 


336 313 283 245 

357 332 298 253 200 = 

360 336 302 256 

360 336 304 258 

382 355 317 265 

409 383 344 285 200) 

415 390 352 292 

416 392 356 297 
477, - 475 469 459 442 418 376 302 
505 503 498 489 475 453) 418 352 
508 507 503 495 483 464 435 378 
508° 507 503 495 485 471 444 405 
543 542 538 552 523 510 491 455 
568 567| 564 560} 553 543 531 520 
574 5735 571 567 562 555 546 559 
576° 575 574 569 567 564 558 549 
613 612 612 611 611 611 618 627 
635 635 635 635 636 638 646 665 
642 642 642 643} 645 649 659] 680 
646" 646 646 648 652 658 671 690 
682 685 685 688 694 701 717 747 
703 704 706 710 717 728 746 7718 
710 711 713 717 725 737 756 7288 
AS. 716 718 721 ou 743 763 794 
750 751 755 761 77a! 786 808 840 
769 | 770 774 781 791 807 830 864 
776 7177 780 787 798 614 858 871 
781° 782 784 790 803 819 844 875 
B13} 814 819 827 839 856 879 912 
830 832 837 845 857 874 898 931 
B37 858 842 850 863 880 904 956 
342! 842 844 852 867 884 908 938 
371 872 877 885 898 915 938 968 
386| 8a 893 901 914 931 954 983] 
392 895 898 906 919 956 958 987 
396. 696 901 909 922 939 959 987 
166 972 991 1025 
'78 peeves 1003 1036 
82 | 988 1006 1038 
85 991 1007 1038 
Figure 2. 


figure is uncertain by more than +1. Difference formulae, of higher than 


second-order accuracy, may then be used, or the net may be subdivided. 


It is clear that the cylinder edges, figure 1, are singularities for ¢, and it can 
be shown that near them ¢ varies as the square root of the distance from the 
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edge. Hence in the neighbourhood of sharp edges, the higher-order residuals 
cannot be neglected. Methods for dealing with this singularity have been evolved 
which are based on the behaviour of the potential near the singularity, but in this 
paper the simplest expedient has been used. It consists-in the use of two suc- 
cessive finer nets as shown in figure 1. When the net is subdivided in a particular 
region, the computation is carried on in alternative steps (i) and (ii). 


(i) Values of the potential at fine mesh points on the boundary between the 
fine and coarse net region are found by interpolation. ¢-values on this boundary 
are then kept constant and, with these values, ¢ is found at points of the fine-mesh 
region by the relaxation process. Some of those points belong also to the coarse 
net. 


(ii) ¢-values at these latter points are next kept constant and new ¢-values in 
the coarse net are found. ‘I'he fine-net values are redetermined according to 
(i), etc., until the process converges. he result of the computation with the 
coarse, medium and finest net are retained. Figure 2 shows the middle region 
of figure 1 with the medium net. The three successive d-values are recorded 
and the fourth figure is found by parabolic extrapolation to mesh-length zero. 
At points not too near the edges the convergence is good. The extrapolation 
has only been carried out in a region where electron paths have been computed. 
In any practical focusing system the wall thickness of cylinders is finite. |The 
mesh length of the finest net should be of the order of the radius of curvature of 
the edge. It will be seen later that the focal lengths of a lens with a wall thickness 
given by our finest net differ insignificantly from those of an ideal lens having 
cylinder walls of zero thickness. 


a 


§3. TREATMENT OF THE REGION OUTSIDE THE CYLINDERS 


In any practical system, electrostatic shields will be found outside the focusing — 
electrodes. We shall later compare the results of our path computation with — 
those given by Maloff and Epstein (1938) and by Spangenberg and Field (1943). — 


Neither of these papers gives information about shielding, and we have, therefore, ~ 


not included a shield. 'To treat the infinite space outside, we have carried out 
an inversion with respect to a cylinder of radius r7,=9a. A new variable, 
s=7/7, is introduced. 

Equation (1) becomes 


a4 st OS Ob 


at oe bree | SF See (4) 

‘The mesh length 6 of the s-net outside the inversion cylinder was taken as 6=} | 

(s is non-dimensional). 
Corresponding values are then : 


vk 
E 4 2 
r= 0 36a 18a 12a 9a 


The cylinders are long, but finite, compared with D, and D,. At large distances, — : 
the field due to the cylinders will behave as a dipole field. Thus at r= 00, b=0. — ) 
The effect of the end charges of the cylinders is neglected, but it is clear that these _ 
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effects must cancel at radial distances small compared with the length of the 
cylinders where we wish the potential to be accurately determined.* 
The finite-difference equations corresponding to (4) become 


16 a 16 2 SZ . 
( ayo =i") bo+ (sree 3] or’) ost ont br - (Fien*+2) gr =0. 


There is a slight ambiguity on the line r=9a, where (2) or (5) may be used. 
Both (2) and (5) are difference equivalents of (1), and we have used the mean of 
the coefficients of these equations. The result of the potential calculations is 
shown in figure 3. 


§4, DETERMINATION OR DHE PAE 


Two paraxial electron paths have been computed and one non-paraxial one. 
In figure 1 one paraxial path is seen entering the smaller cylinder at r=D,/16. 
Another path entering the larger cylinder in a direction parallel to the axis 
at the same radial distance has also been computed. ‘These two paths determine 
focal length and principal planes. ‘The non-paraxial path gives a measure of the 
spherical aberration in the forward direction accelerating lens. 

The equation for the path of a negatively charged particle is 


[ pe (5 ny dr 
dr av dz) ida la, dr\? 
Aad As WV Tea: 1+(F) | -°. sraemerree (6) 
The method of path computation is based on a process suggested by L. J. Comrie 
in Interpolation and Allied Tables (1936). 

Define i dp = dp 

Sate = av 

where p=7/a, €=2/a, anduse steps 


Pr — Pn =Fnt1 — Fn = 4: 


Equation (6) can be written 
146 14f 
v(t) av(i+4) 


F= -— = 
ee ie Nese 24 ap pee ee (7) 
The procedure consists in computing entries of the following table :— 
Table 1 
gla p Gs ae ae if F Ae AGS 
ee eee eee 

0 Po . d,”’ fo Fy A,” 
aye dye” Aj)2" 

1 Px qo a,” fi Fy § 
3/2 dsj9"" Asjo’ 

2 Pe ; da’” Se Fy, a 
dsjo d5jo/” A5)9' 


* Alternatively one might interpret the field as due to infinite cylinders and a zero potential 
hield at infinity. 
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where d’, d’’, d’” are successive central differences of p, and A’, A”, A’, ... are 


successive central differences of F. 
Assume fp, P1, P> known; in our example 


Po = Pr = ps = D,/16- 


Calculate f,, f,, f. from, for example, 


f= (2a) = 3 (2a) + 5 (24, ( oe Pe, et ee (8) 


and Fy, F,, F, from equation (7) and difference. A,’ is already found in this 


way, but not A,”. 
Las 1a 
VC Vecn whe) cra (9) 


are taken from the field plot, figure 3. (Tables of these values are prepared before- 
hand so that appropriate values of (9) can be interpolated for each value of p. 
For interpolation we use a four-point formula (Comrie, 1936, p. 685)). 

From now on the following steps are carried out: 


(i) Estimate A,”, e.g. take A,’’ =A,” and calculate d,’’ from 


t 1 1 
ae ies a C1 i ee (iv) 
(= i( f+ 173 42 40 M2 + ) eee, (10) 
(ii) Form ps; by building up from d,”’. 
(ui) Verify by calculating f, from (8) and, if necessary, make a correction. 


(iv) Now extrapolate fo=h+Fet 3 Na! ees ap he ae nee cere (11) 


(v) Compute F, and difference; 
(vi) Verify that correct value of d,” has been used in (ii) ; correct if necessary. 


The steps of z used in the path computation are conveniently chosen so that 


ao A,” is not larger than a few units of the last decimal figure of F which is 
retained during computation. It should be assumed that one significant figure 


is lost in about ten steps of the integration. It is therefore necessary, if absolute 
certainty of the result is desired, to retain during the computation more decimal 
figures than the final accuracy claimed for the path warrants. 

If, in our table 2, the last figure of F is uncertain by +1, then the 7 might, in 
any step, acquire an error of one-third of the last figure retained. ‘There are 30 
steps; the errors accumulate in arithmetic progression so that, if all the errors are in 
the same direction, the final result may be in error by 0-150 and the position of the 
focal points in error by two mesh lengths (cf. the end of the path given in table 2). 
Now it is very unlikely that all the errors are in fact in the same direction, and we 
estimate that the error is less than one mesh length. 

The process outlined above leads to table 2 for the paraxial path shown in 
figure 1. We have found it convenient to use the steps of the relaxation net for 
the pathcomputation. Thusthe first eight steps are carried out in the coarsest net. 
To start the path in the middle region medium net, two values of p are inter- 
polated at mid-points of the last two steps. The last part, again in the coarse net, 
is started by selecting the last three values of p belonging to the coarse net in the 
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middle region. At the end the path is a straight line and the first differences are 
0-147. The point of intersection is found by interpolation. The paraxial path 
was first determined in the coarse net throughout, and F,, the distance of the focal 


Table 2 
t 0 d’ qd’ dl” f F A’ A” 
—10 0-500 0-000 0-000 0-000 
0-000 0-000 0-000 
— 9 0-500 —0-:001 —0-:001 0-000 —0:002 
—0-001 0-000 —0-002 
— 8 0-499 0-000 —0:002 —0-002 —0-001 
—0-001 —0-001 —0-003 
— 7 0-498 —0-001 —0-003 —0-005 —0-002 
—(-002 0-000 —0-005 
— 6 0-496 —0-002 —0:006 —0-010 0-000: 
—0-004 —0-001 —0-005 
— 5 0-492 —0-004 —0:012 —0-015 —0-:003. 
—0-008 —0:002 —0-008 
— 4 0-484 —0-006 —0:021 —0-023 —0-001 
—0-014 —0:001 - —0-006 
— 3 0-470 —0:007 —0:035 —0-029 +0:-004 
—0-021 —0-001 —0-002 
— 2 0-449 — 0-006 —0:048 —0-031 +0-009 
—0:027 0-000 + 0-007 
— 1 0-422 —0-006 —0:061 —0-024 +0:-015 
—0-033 + 0-006 +0-022 
QO 0:389 0-000 —0:067 +0-002 —0:010 
—0-033 ~ +0-004 +0-012 
+ Ll 0-356 +0-004 —0:063 +0-014 —0-012 
—0-029 —0-001 0-000 
Se PP (hea 7 + 0-003 —0:055 +0-014 0-000 
—0:026 +0-001 0-000 
So 0-300 + 0-004 —0:048 +0-014 +0:004 
—0:022 —0:002 + 0-004 
+ 4 0-279 --+0:002 —0-042 +0-009 —0-002 
—0-020 0-000 —0-002 
== & Oevs© -+0-002 —0:038 -+0-007 +0:-001 
—0-:018 —0-:001 —0-003 
+6 0-241 +0-001 —0-035 +0-004 +0:002 
—0:017 0-000 —0-001 
=F 7) 10°294 +0-001 —0:034 +0-003 0-000 
—(0-016 0-000 —0-001 
+ 8 0-208 -+-0-0005 —0-032 +0-002 0-000 
—0:0155 0-000 —0-001 : 
== 99) 011995 --0-0003 —0:031 +0-0012 +0-001 
—0:0152 0-000 0-000 
+10 0-1773 +0-0002 —0-030 -+0-0007 0-000 
—0-0150 0-000 0-000 
+11 0-1623 +0-0001 —0-030 +0-0005 0-000 
—0:0149 0-000 0-000 
+12 01474 +0:0001 —0:030 +0-0002 0-000 
—0:0148 0-000 0-000 
+13 0-1326 0-0000 —0-030 +0-0002 0-000 
—0-0148 0-000 0-000 
+14 0-1178 —0-:0001 —0-030 +0-0002 0-000 


—0:0147 etc. 
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Table 2—continued 


Medium region computed in half-intervals 


Zu 20 d’ Oe dae if F A’ A” 

—/  10:956 —0-000 —0:028 —0-013 interpolated 
—0-016 —(-000 —0-002 

— 6 0-940 —0-004 —0:-036 —0-015 +-0-001 
—0-020 +0-002 —0-001 

= 2 Une —()-002 —0:042 —0-016 interpolated +0-001 
—0-022 —0-002 0-000 

— 4 0-898 —0-004 —0-048> —0-016 +-0-003 
—0-026 +0-001 +0-003 

— 3 0-872 —(0-003 —0-055 —0-013 —0-002 
—0-029 0-000 +0-001 

— 2 0-843 —0-003 —0-061 —0-012 +0:005 
—0-032 +0-001 -+0-006 

lO oan —(0-002 —0-067 —0-006 +0-001 
—0-034 +0-002 --0-007 

OM O77 0-000 —0-069 +0-001 —0-002 

—0-034 +0-002 +0-005 

+ 1 0-743 +0-002 —0:066 +0-006 —0-004 
—0-032 0-000 +0-001 

Se Dersilil +0:002 —0-062 +-0-007 end path start 0-000 
—0-030 0-000 +0-001 

S793) 0-684 +0-002 —0-058 +0-008 —0-002 
—0-028 0-000 —0-001 

= 4 =0°653 +0-002 —0:054 +0:007 end path 0-000 
—0-026 —0-001 —0-002 

== 5 Oxy +0-001 —0:051 +0-005 
—0-025 

On O02, ; —0-049 end path 


point from the cylinder junction, was F,=+15. Then it was redetermined 
with the more exact figures of the medium net, and finally again with the best 
extrapolated figures; the last two results for /, were the same, although the paths 
were slightly different, and we therefore considered our result as the final answer. 
The distances of the focal planes F,, F, and the cardinal planes H,, H,, from the 
plane through the cylinder junction, are given in table 3. 


§5. COMPARISON WITH THE RESULTS OF OTHER AUTHORS 


Maloff and Epstein (1938) on p. 105 give data for a lens similar to the one 
computed in this paper with D,/D, =1-5, V,/V,=6. 

There is no precise information in the book concerning the method by which 
these results were obtained. In particular, no details are given about position 
and potential of shielding electrodes outside the lens either in the electrolytic 


trough experiments determining the field or in the experimental verification 
of focal length, etc., given on p. 114. ‘It is, however, likely that one grounded 


shield was used. The wall thickness of the cylinders is not stated. All these 
data are also omitted by Spangenberg and Field (1942), p. 195, who give a chart 
for a similar lens with the same parameters. It must be concluded from the 
nature of the experimental arrangement that two shields at different potentials 
were used. The results taken from these papers are summarized in table 3. 
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It is difficult to say how much of the disagreement is due to the difference in 
shielding and wall thickness in the three cases. _ It will be seen that the distances 
between focal points agree wthin one mesh length, and this is within the accuracy 


Tables 
F, F; A, A, 
Maloff and Epstein —21°6 20 —8-8 —13-2 
Spangenberg and Field —21-68 18-11 
Our result: Paraxial — —18-56 22 —6°37 —12-01 
Non-paraxial 1925 


of our calculation. In our result the whole system of principal planes and focal 
planes seems to be shifted bodily by about three mesh lengths towards the high- 
potential electrode. 


P 


Figure 4. 


Spangenberg and Field (1942) give on p. 312 a graph for the reduction of 
length for a non-paraxial path in an accelerating lens. They state that this 
reduction curve is almost universal if percentage reductions are considered. 
We have computed a non-paraxial path (figure 1 and table 3). The reduction 
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in focal length is 87-5°%, whereas the reduction read off the universal curve 
mentioned above is 87%. 

In order to obtain the full spherical-aberration curve, several other non- 
paraxial paths were computed, and the curve obtained is given in figure 4. 


We conclude that our path determinations are internally consistent and that 
the differences between our results and those of the authors quoted are mostly 
due to errors in the field computation and the uncertainty of the geometrical 
data. We could, of course, increase the accuracy of our field determination. 
Such a further effort is not likely to effect our results materially and does not seem 
justified in view of the experimental errors. 
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A MODIFIED HELMHOLTZ LINE-ELEMENT IN 
BRIGHTNESS-COLOUR SPACE 
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ABSTRACT. Helmholtz, Schrédinger and others have considered the relation between 
the trichromatic coordinates, x, y, x and x+4x y+8y, z+ 4s, of two juxtaposed light 
patches which can just be discriminated by the eye. The Helmholtz expression, 


sx \? oy . sz \" 
a |) def eee), te = const. = 8s?, 
atx b+y cte2 


is the most readily interpretable in terms of the trichromatic theory, but to explain the 
observed variation of the hue limen through the spectrum it demands double-peak 
fundamental response curves which are at variance with other evidence. It also leads 
to a step-by-step visibility curve of wrong shape. The difficulties mentioned are removed 
by introducing different constant factors in the three terms. Such a modification is 
indicated independently by recent measurements of the liminal brightness increment, 
from which appropriate values of the factors have been derived. The factor in each term 
is simply related to the limiting Fechner fraction of the corresponding trichromatic 
mechanism, i.e., the Fechner fraction 4B/B which would be observed at sufficiently high 
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brightnesses B if visual discrimination depended on that mechanism only. The modified 
expression for ds? is applied to the calculation of hue limens, the step-by-step visibility 
curve, the Fechner fraction curves, and the general colour limen at different points of the 
colour triangle. Certain main features of the experimental results are correctly reproduced, 
but discrepancies with Wright and MacAdam’s measurements of general colour limens 
may indicate that the modified element ignores some factor which is operative in their 
experiments. 


§1. INTRODUCTION 


HE examination of the properties of the proposed line-element will be 
preceded by a brief comment on earlier work. It was Helmholtz (1891) 
who first attempted to express in mathematical form the condition for 

two juxtaposed light patches to be on the limit of discrimination by the eye. 
He suggested that any two just-distinguishable patches would have trichromatic 
co-ordinates x, y, z and x+6x, y+dy, 2+62z respectively, satisfying a relation of 


the form, ; : 
a by Oza 
36s?=( — } 4+ (| — —— |) = 30) peer Hs 
2s =(3) . (4) as (=) cs) 


where a, 6, c and F are constants. The trichromatic co-ordinates here are the 
actual quantities of a particular set of primaries required to match the light 
patch—not the actual quantities divided by their sum (the unit co-ordinates). 
It is customary to regard the quantity 6s defined in (1) as the length of the line- 
element in a non-Euclidean space of co-ordinates x, y, 2. In this so-called 
brightness-colour space, the points corresponding to any pair of just-distinguish- 
able light patches are the same elementary distance F apart. If a change is 
made to a different set of primaries and the relation (1) is expressed in terms of 
‘the new co-ordinates of the light patches, the simple form of the relation is lost 
and it is no longer the case that each term depends on one co-ordinate only. 
Thus, if the original form is valid at all, it is valid for only one set of the infinitely 
many sets of primaries which would serve equally well for the statement of the 
results of colour-matching measurements. The Helmholtz line-element implies, 
in fact, a particular set of primaries which, on the usual physiological interpre- 
tation of the trichromatic theory, we should expect to be the fundamental set 
such that each primary stimulates one, and only one. of the cone mechanisms. 
The quantities of the fundamental primaries required to match a monochromatic 
patch of unit energy intensity define three functions of wave-length—the spectral 
sensitivity curves or fundamental response curves of the three mechanisms. 
Helmholtz took Kénig’s measurements of the trichromatic co-ordinates of 
the spectrum colours and sought to express them in terms of a set of primaries 
for which the relation (1) would be valid. The test was that the values of the 
hue limen through the spectrum as derived from (1) should agree with the measure- 


ments of this quantity made by Kénig and Dieterici (1884). Helmholtz actually 
used the simpler form, . 


\2 2 2 as 
35% = (*) + (~) a (=) 22378) ) =e ee a2) 
x y g 


to which (1) reduces at sufficiently high intensities when a, b, c can be neglected 
compared with x, y, x respectively. He was able to find a set of primaries which 


x ‘The factor 3 is introduced so that F has a simple physical meaning, as explained below. 
Relation (1) is sometimes referred to as the threshold condition. 
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satisfied this test fairly well. But the corresponding fundamental response 
curves (reproduced as figure 8 in Peddie’s Colour Viston) each have two pro- 
nounced maxima and are quite unlike the Grundempfindungen proposed by Kénig 
from a consideration of the properties of colour-blinds. Moreover, as Schré- 
dinger (1920) pointed out, they lead to an impossible double-peak form for the 
step-by-step visibility curve (the curve showing the relative energies of neigh- 
bouring spectrum colours which appear equally bright). 
Schrodinger raised a further objection to the Helmholtz line-element. Although 
the step-by-step visibility curve can be derived from the line-element and its 
associated fundamental response curves, some additional assumption is necessary 
to determine when the eye will judge as equally bright two juxtaposed patches 
which differ widely in colour. Schrédinger assumed that two patches would 
appear equally bright if any change in the intensity of one of them would increase 
the minimum number of just-perceptible steps needed to pass from one to the 
other. ‘The minumum number for any two light patches A and B is proportional 
to the integral 
-B 
| ds 
A 


taken along the geodesic in b.c. space. ‘Thus the assumption is that A will 
match B when the intensity of B is adjusted to make this integral a minimum. 
While mathematically attractive, this definition of equality of brightness is not 
very plausible from a physiological standpoint. Applied to the Helmholtz 
element, it makes brightness a non-additive property of lights, in conflict with 
Abney’s law. Schrodinger proposed an alternative line-element, 
ne E jaar | SHU peat (3) 
ax+By+tyz x y z 

which, on his assumption, makes brightness proportional to ax+f/y+ yx and, 
therefore, strictly additive. An adequate criticism of this line-element will 
not be attempted, but it appears to the writer inconsistent with the threshold 
measurements considered below. 

Since Schrédinger’s theoretical work, several investigators have determined 
by experiment sets of pairs of just-distinguishable patches on various lines and 
surfaces in b.c. space. Most of these measurements were made at fairly high 
intensities and were confined to a particular class of just-distinguishable pairs, 
namely, those in which the members of each pair appeared equally bright (colour- 
limen type*). The observations on which the modified line-element of this paper 
is based were of a rather different nature, and are best introduced by considering 
the application of the Helmholtz element to another special class of just- 
distinguishable pairs of patches (liminal-increment type). 


§2. APPLICATION OF THE HELMHOLTZ -LINE-ELEMENT 
MOLPRAIRS OF JUST-DISTINGUISHABLE PATCHES 
OF LIMINAL INCREMENT TYPE 


A pair is of this type if the physical stimulus in one patch is obtained from 
that of the other by the addition of a further stimulus. Consider the special 


* Measurements of the hue limen amd the minimum perceptible colorimetric purity (Priest and 
Brickwedde, 1926) are of this type, as well as the more general colour limens recently determined 


by Wright (1941). 
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case when one patch is monochromatic, of wave-length » and energy intensity 
W,,, while the other is a mixture of wave-lengths » and A (A in general different 
from p), of energy intensities W, and U, respectively. ‘The threshold condition 
(1) takes the form, 


Uarg Z | Ung : | U jb, é 
—_—, —2___ | + | ————| | =3F"*, ...... 4 
li | + | ba Wag, aan (4) 
where 7’, g;', b;' are the fundamental response functions of the’'x, y, 2 of redie 
“ greeny’, “blue”? mechanisms respectively. Equation (4) may be written 
Love 
G ) = Fr RaW, JP + Lei G.W,)P + (Bin BuWae ee (5) 
or : 
i} 2 ] 2 ii 2 ( 1 ) 
ral) pe ee ee 6 
ie (a) 3 ew Un °) 
where 
= 1 
uf (P) = 14+9p 
and 
R= 11 |aFV3, R,=7;'/9a, on = 7179(RLW,), 
% 1 | 
g1=81 [DF V3, Gi=g, /9b, Tg, 8GuMs), > ae (7) 
= 1 
b,=by |cFV 3, B,=6,' /9e, Ty =O BuW). 
Ab 


} 


U;,, Ujg amd Uj, may be regarded as the values which U, would assume if, in 
turn, each of the three mechanisms were acting alone. Since 7(0)=1, then, since 
= =r, when W,,=0, and 7, is the reciprocal of the absolute threshold (in energy 
Ar 

units) of the red mechanism for light of wave-length A. Also, since 7(1)=0-1, 
R,, is the reciprocal of the energy intensity W, of light of wave-length required 
to raise U,,, the difference threshold of the red mechanism, to ten times the 
absolute threshold (1/r;). Similar interpretations apply to gj, Gy, b;, By. 
For A=p, we have . 


at eee (8) 


V2 a 


For fixed A and y it is convenient to illustrate relation (5) by a plot of log (1/U,) 
against log W,, and to show plots of log (1/U;,), log (1/Uj,) and log (1/U;) 
in the same diagram (figure 1). The latter curves (component curves) all have the 
same shape, but their respective positions in the diagram are determined by the 
values of 7, 2, 6, (position on the ordinate axis) and R,, G,, B, (position on the 
abscissa axis). If =A, the component curves all tend asymptotically at high 
values of W,, to the line : ordinate =log (1/F/3) — abscissa, while the resultant 
curve (log 1/U,) tends to the line: ordinate=log(1/F)— abscissa. If uA, the 
component curves may intersect so that the resultant curve exhibits a more or 
less pronounced “ change of law” as illustrated. 
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§3. BASIC MEASUREMENTS 


Given the seven functions 7(p), 7... ... B, the resultant curve for any 
combination (A, ) can be drawn at once. An attempt at the converse problem 
of determining the seven functions from a number of experimental curves of 
log (1/U,) against log W,, has recently been made. The main points of this 
work (Stiles, 1939) will be briefly indicated. The present application was not 
in mind and the conditions of observation differed considerably from those 
implied in the foregoing discussion. 


(a) The subject viewed, not two small juxtaposed patches but a large (10° diam.) 
patch of wave-length and energy intensity W,,, at whose centre a small (1° square) 
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Figure 2. 
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Figure 1. 


patch of wave-length A and energy intensity U;’ was applied in flashes (0-063 sec.) 
as an additional stimulus. ‘Thus within the 1° square during the flash the energy 
intensity equalled W,,+ U,’, while that of the surrounding field remained at 
W, throughout. The subject fixated the centre of the 10° patch, and great care 
was taken to ensure that the observations corresponded to foveal and cone (rod- 
free) vision of the additional stimulus. By repeated exposures of the latter at 
different energies U;’, W, being kept constant, the critical intensity U, was 
determined at which the subject had a 50% chance of seeing the 1° patch. 


(6) By measuring U, in this way for a series of values of W,, from zero upwards, 
a curve of log(1/U,) against log W, was obtained for various pairs of 


46 W. S. Stiles 


wave-lengths A and y.. These experimental curves corresponded to the broken- 
line curves in figure 1. For example, with A=480 mp, »=540 mp, the experi- 
mental curve closely resembled the broken-line curve of figure 1(b). It was the 
characteristic form of such experimental curves which first suggested that they 
were really resultant curves derivable from component curves of fixed shape. 
This view was supported by the fact that when A was kept constant at 480 mu 
and p was varied, the experimental curve was modified as though (referring to 
figure 1(b)) the component curves log (1/Ujy) and log (1/U%) were displaced 
parallel to the axis of log W,, while remaining in the same position along the 
axis of log (1/U,). On the other hand, when » was kept constant and A varied, 
the experimental curve was modified as though the component curves were 
displaced parallel to the axis of log (1/U,). These regularities in the results may 
be called the ‘‘ displacement rules ”’ 

From nearly 100 curves of log (1/U,) against log W, for the writer’s left eye, 
it was concluded that all such curves for this eye could be regarded as resultants 
of three component curves of a common fixed shape whose positions in the 
diagram depended on A and p. ‘The most satisfactory fixed shape to suit all the 
results was chosen by trial, and the positions in the diagram at which component 
curves of this shape had to be placed to yield a resultant curve agreeing with the 
experimental curve were determined for the various combinations of A and pu 
studied. Actually no more than two component curves contributed to the 
experimental curve for any particular combination, but to cover all combinations 
in accordance with the displacement rules three components were necessary. 


(c) The common shape of the component curves was represented as a function 
€(p) (figure 2) defined so that &(p)=1 for p=0, &(p)=0-1 for p=1. Thus &(p) 
played the réle of the function 7(p) above. It is apparent from figure 2 that the 
difference between the two functions is small and confined to values of p 
below 1. The expression fora particular component, say the ‘‘ red’’ component, 
then took the form, 


log (1/U,) =log [r:€(R,,W,,)] 


where for given A and py the constants 7, and R, fix the position of the curve in the 
diagram. Conversely, when the position had been determined as explained above, 
the values of r, and R, were obtained. Although this method yielded values for 
each of the functions 7;, R,, ga, etc., in a part only of the visible spectrum, the 


ranges were extended by assuming the ratios Ri : Ga Bi to be approximatel 

Pi) Wi. aby PP ¥ 
independent of A. This appeared to be nearly true in regions where both numer- 
ator and denominator functions could be followed. The tentative determina- 
tions of the six functions 7, to B,, are plotted as the continuous-line curves of 
figures 3 to 5. ‘The absolute values depend on the unit in which energy intensity 
s expressed. Here and elsewhere in this paper, by energy intensity is to be 
understood the flux of radiation entering the eye, expressed in ergs per sec, 
per square degree of the light patch. ‘The values given refer to the case when the 
pupil of entry at the eye is small and concentric with the natural pupil. Rather 
different values apply when the pupil is large or eccentric, on account of the 
directional properties of retinal sensitivity. 
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R, G; B 
(d) On the assumption that es and - are constants independent of 4, 
Ne) A 


the values of these ratios were found to be (very approximately) in the proportion 

Ry Ga B 

—: 2: 2420-78: 1: 4-46. 

Yr), Si b, 

(e) The resultant curve was related to the component curves by the following 

rough rule: for a given log W,, log 1/U,... equals the greatest of log (1/U,) 

log (1/U,,) and log (1/Ux), except when the greatest and next greatest have 

nearly the same value. In the latter case log (1/U,) exceeds the greatest by a 
small amount, of the order of 0-1. 
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(f) A breakdown in the analysis of the results on the above lines occurred 
when the 10° patch was orange or red and of very high intensity (about 30,000 
photons) and when the additional stimulus was blue (limited conditioning effect 
of red light on the blue mechanism). Colour-match equations also break down 
after adaptation to intensities of this order (Wright, 1936). 


Setting aside the effect noted in (f), the experimental results differ from 
the predictions of the line-element and associated primaries of eke Faroe in the 


following respects :— 
(1) The function 7(p) is replaced by a slightly different function &(q). 
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‘ 


(ii) The spectral sensitivity functions have the same general form as Konig’s 


Grundempfindungen. 
(iii) The ratios nae ay @ are not equal: Bi ig materially greater than 
Ty ep 103 b; 


the other two. It follows that for A=, the three component curves 
in diagrams such as figure 1 no longer tend to a common line at high 
intensities. 


(iv) The effects of the three mechanisms add up, if anything, to a smaller 
extent than required by the sum of squares relation (6). 


I 


aeericr X +417 7] 


45 


A_OR L_IN MILLIMICRON, 
ACORN MILEIMICRON, 


Figure 5. 
For A=y, the ratios Un Ug Un are the Fechner fractions * of the respec- 
lW. A : W, : W, P 

; 1 
tive mechanisms. Since &(p) tends to ® as p becomes large, it is clear that 
Uy, 1 R 
ie War t(RW,) approaches sles at high intensities W,;. ‘Thus the ratios 
R, G, B, 
Tee are proportional respectively to the limiting Fechner fractions of the 


Historically it would be more accurate to call these ratios ‘“‘ Weber fractions’, but it has 
seemed preferable to follow the more usual practice. 
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three mechanisms, and (iii) implies that, contrary to the assumption of Helmholtz 
the different mechanisms have different limiting Fechner fractions. 

‘The precision of threshold measurements is low, and during a long investi- 
gation the properties of the eye may undergo systematic changes. There can 
be no question of a precise determination of the unknown functions by the 
methods used. All that can be claimed is that the derived functions reproduce 
the main features of the experimental curves. 


§ 4. MODIFIED LINE-ELEMENT 


Following the hints given by the basic measurements, we may construct 
a tentative line-element for application to any pair of just-distinguishable light 
patches, not necessarily of liminal increment type. Suppose the patches P’ 
and P have absolute energy distributions W;’dA and W,dA respectively. Put 


((p)=94(p),  R=[WiRdd,  3R=|(Wi- WR ......, (10) 


and define G, 5G, B, dB similarly, where €(p), R;, G;, B, are the experimental 
functions* plotted in figures 2to5. Introduce the quantities p, y, B proportional 
to the limiting Fechner fractions of the three mechanisms, 


hap ieee Oh BS 
Be area gee eo: Leesan (11) 
and satisfying 
ia sl 1 = 
B oF Z te eo OS ae (12) 
Be 0612 = = 0 309, : 
so that ao = j z 200, Fe =(0-0185. The proposed line-element - is 
then 
: dR e dG ‘ dB Z 
ast =| = 208) + [Sua] +| Fae) | SW wr Ne (13) 
p y B 
which reduces at high intensities to 
13k? 1oG|? 16B]? 
Se? =] - pies ee || a ae 
s - Z| +[ = = | +35 | | cits an Me MD Aeon (14) 


where F is a constant. 

If the two patches have the same relative energy distributions, their inten- 
sities may be specified by their total energies T’ and T (or by their total energies 
weighted according to any function of wave-length, such as the visibility curve). 
Their Fechner fraction is then given by 


Of eal E> 8 R. 0G. 0B 


Thus, by (14) and (12), a =F at high intensities, or F is the limiting Fechner 


fraction for any pair of patches of the same relative energy distribution whatever 
form that may take. The factor 3 introduced into (1) gives Fa similar meaning 
for the Helmholtz line-element. Experimentally, the value of F depends on the 


* The definitions of R,, G,, B, are modified later. 
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precise conditions of observation. It is round about 0-01 for the conditions 
normally used in determining colour limens, and this value will be used in the 
calculations which follow. 

For just distinguishable patches of liminal increment type, 6R, dG and dB 
are all positive or all negative. In applying the line-element generally, it is 
assumed that the signs of 5R, 5G, 5B are immaterial, so that, for example, slight 
excess stimulations of the red and green mechanisms make the same contribution 
to the limen whether the excesses occur both in the same or one in each patch. 


§5. COLOUR MATCHING 


Light patches whose representative points P’ in b.c. space lie within the 
small ellipsoid centred on P and defined by (13) are indistinguishable from, 
and therefore in colour match with, P. In colour-matching measurements 
the objective is to bring the variable colour to the centre of this ellipsoid, 
and the condition for an ideal match may be taken to be 'R=dG=6B=0, or 


fwrr, ar=|WAR; dX, and two similar equations.* R,, Gj, B, were, in fact, 


obtained by a method completely independent of colour-matching, but their 
introduction into the line-element relates them at once to the distribution 
coefficients *, y, % derivable from colour-matching measurements. In terms 
of these coefficients, the conditions for a colour match are 


| Ho De | Wen 


and two similar equations. Thus R,, G,, B, should be linear forms in &, y, 2. 
The distribution coefficients for the writer’s eye (a normal trichromat) have not 
been determined, but using the coefficients of the C.I.E. Standard Observer 
in the Standard Reference System (C.I.E. 1931) an approximate linear repre- 
sentation of R,, G,, B, is obtained as follows :— 


(Ri) 6-61 . 102 #+ 1-26 . 103 — 1-12 . 102g, 
(G;) —4-38. 102% 41-62. 1037 +1-23. 108%, Fo... (15) 
(Bi) 7-08. 104% +0. 5 +417. 1078, 


These linear forms are plotted as the broken lines in figures 3 to 5. The agree- 
ment between the broken and continuous curves is fairly satisfactory when it is 
recalled (a) that the experimental uncertainties in determining R,, Gj, B, are 
considerable—much greater than for x, y, 3 ; (b) that results for a single eye and 
for an average eye are being compared. The basic measurements (§ 3) 
could be represented almost as well by using the expressions (15) in place of 
R,, Gy, Bj, and by making the same substitutions in the definitions of R, 5R, etc., 
the line element would lead automatically to the correct colour-matching equations 
for the C.J.E. standard observer. From now on, therefore, R;, G;, B, are supposed 
defined by the expressions (15). We are, of course, still dependent on the 
basic measurements for the coefficients in (15) and for the relative values of 
ps ¥> B. 


* Even if an ideal colour match is not achieved, (13) ensures that these equations shall be 
approximately true, i.e. 5R will be sinall compared with R, and so on. 
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§6. STEP-BY-STEP VISIBILITY CURVE 


In the experimental determination of the visibility curve by the step-by-step 
method, the two juxtaposed patches have slightly different colours and the 
intensity of one is varied until the “difference” between the two patches is a 
minimum. If the colour difference is appropriately chosen, the two patches 
are just indistinguishable at this minimum setting of the intensity. The visibility 
curve is derived from the line-element (14) by a corresponding procedure. The 
unit co-ordinates (in the fundamental system) of the two patches, 7’, j,k’ andz, j,k, 
are required, where 


a ne Bike oe His mk 16 
RBGPERORE Pe) eo? 1 RAGES ose eo. 
andj’, k’ ; j, k are similarly defined. 
In terms of these, the line-element (13) takes the form 
elie te a es S 
r= ¢(z) | + twosimilarterms. ...... (aa 
pP 


By varying &’, keeping z’, 7’, k’ constant, the intensity of patch P’ is varied without 
altering its colour. The minimum value of 5s? under such a variation is reached 
when 


Dene Seal ee orocimnilarrerne 
p 
or, reverting to the earlier notation and retaining only first-order quantities, 
when 
dRER : rake if a15 Ih 3 
= —|—C(R)| + a1 -—&G)] +s1sc(b)]. ...... 18 
0= 8) Fam] + Gl ou@] +3] Ge (18) 


Strictly, the colour difference of the two patches should be chosen so that at the 
minimum setting 5s= F, but it suffices both in the theory and in the experimental 
method for this condition merely to be approached. . 

The equality-of-brightness condition (18) applies to any pair of patches of 
nearly the same colour. In the special case of monochromatic patches of wave- 
lengths A, A+ 5A and energy intensities W;,, W;,+6W,, (18) becomes 


BWa _5,[Ce dR: , CPdGs | Cy? dB,), 
WW, i, | CRIN I ER 


SO ee 
© [Baa] +[San] +[72@] 


and C,?, C,? are similarly defined. By successive applications of (19) poten 
determine the energy intensities of a series of monochromatic patches w on 
wave-lengths increase in small steps from the blue to the red, and each of De 
matches in brightness its next neighbours in the series. ‘The reciprocal of the 
patch energy plotted against wave-length represents a step-by-step eae eye 

whose shape will depend in general on the intensity level, specified, say, by the 
; Hes 


where 
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value of W, atX\=555 mp. At sufficiently high intensity levels C,?, C,?, C,? tend 


. olf 
to 1/p?, 1/y2, 1/8? respectively, and the visibility curve approaches the limiting form 


const. 2 2 2 
V;, = W = const. R,ve Go BP © — evemec renee us (21) 
A 


As usual, the constant is adjusted so that V, has the maximum value unity. 

The circle points in figure 6 represent relative visibilities calculated from 
the expression (21). Except in the blue, these points lie very close to the full- 
line curve, which is a plot of the C.I.E. visibility function. The discrepancy 
in the blue is greatly reduced if the comparison is made with the mean visibility 
curve obtained by Gibson and Tyndall (1923) for 52 observers (38 on ends) 
using the step-by-step method.* In the several investigations of the visibility 
curve in which large groups of observers have been used, the variations between 
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Figure 6. Step-by-step visibility curves. Figure 7. Fechner fraction curves. 


observers are markedly greater in the blue (say 455 my) than in the red (say 655 
my). On this general ground any theoretical derivation of the visibility curve 
~may be expected to meet greater difficulties in the blue than in the red. 

The broken curve in figure 6 is the visibility function for the Helmholtz 
line-element, but assuming the response functions (15). It is derived from (21) 
by putting 

1 PS iaret 

pt Ge Re Be 
Equally unacceptable curves are obtained by applying to the Helmholtz element 
Konig’s Grundemfindungen or the double-peak response curves of Helmholtz. 
The Hecht response curves (Hecht, 1932), each of which is a fair approximation 
to the visibility curve, would make the Helmholtz or the modified line-element 
consistent with the visibility curve, but would render it incapable of explaining 
the basic measurements (§ 4). 


* I am indebted to Prof. Hecht for drawing my attention to this point. He also points out 
that the mean curve of Coblenz and Emerson (1917) shows a hump in the blue, and that a 
protanope studied by Hecht and Shlaer gave a curve similar in this respect to a normal observer. 
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The standard photometric brightness level at which the C.I.E. visibility 
curve applies may be taken as 25 equivalent metre candles seen through an eye 
pupil of 10 sq. mm. area, or 79-5 photons. A monochromatic patch of wave- 
length 555 my has this brightness when its energy intensity W, equals 3-66 x 10-4 
erg per sec. per square degree. The limiting form (21) of the derived visibility 
curve is barely reached at this value of W;;;, and the curve must be derived 
from (19) by a method of successive approximation. The results of the calcu- 
lation are shown as the cross-points in figure 6. The agreement with the C.L.E. 
curve is substantial, although in the red it is rather less close than for the high- 
intensity curve. As the intensity is reduced below the photometric level, the 
derived curve undergoes some further modification, to assume finally the form 
shown by the triangle-points of figure 6, when the absolute threshold is reached. 
The second limiting form (22) can be obtained from (19) or directly from the 
line-element (13) by using the principle that at the absolute threshold all light 
patches are equally bright : 


V,=const. Re), a (Sy ER (3) 4 ot aveueee (22) 


The main effects of lowering the intensity level are to increase slightly the visi- 
bilities in the orange and red and in the blue. It must be emphasized that the 
line-element (13) is intended to represent foveal and strictly cone vision. The 
small changes in the visibility curve just considered are quite distinct from 
the very much larger changes (Purkinje effect) which occur when, by the use of a 
large or eccentrically fixated matching field, the rods participate in the matching 
process. 

No determinations of the foveal cone visibility curve by the step-by-step 
method have been made at low intensities. Walters and Wright (1943) used 
a fixed red comparison patch and worked with the full-colour difference. They 
found that at low intensities the visibility curve showed a slight shift of the 
maximum towards the blue, which they think may have resulted from some 
admixture of rod vision, and a “lump”? in the curve at about 600 my, which 
they attribute to a tendency of the red mechanism to retain its sensitivity over a 
greater intensity range than the green or blue mechanisms. 

Relation (18) is the differential of the function 


P(R, G, B) = P,(R) + P(G) + PB), 


pany [28 [ep 
oes ep 

and P,, P, are similarly defined. All lights (not necessarily monochromatic) 
for which P(R, G, B) has the same value can be linked in small steps by series of 
intermediate lights, adjacent lights in any series having the same brightness. 
We may say that all such lights have the same ‘“‘small-step ” brightness. But 
P,, P, and P, are not linear functions of their arguments, so that “ small-step” 
brightness is not additive. ‘To estimate the magnitude of the breakdown in the 
additive law, we will use for P(R, G, B) the form to which it reduces at high 
intensities, namely, 


where 


log, [RU Gur Bue", 
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If lights 1 and 2 have the same small-step brightness, 
Ril? G2 BYP = Re GY BP, 


and the equality of brightness still holds if their intensities are changed by the 
same factor, « say, since 


(Rya)¥e" (Giq)u" (Bya)u = (Roo)Ve" (Gaa)4” (Boa)uPr. 


But if lights 1 and 2 each match light 3 in brightness, the mixture of 1 and 2 
will not in general match light 3 increased to double its original intensity, since 


(Ry + Ra)" (Gy + G,)!"(B, + Bo)?" = (2Rg)"" (2G,)¥"(2Bs)" 


is not in general true. ‘The discrepancy may be expected to be greatest when 
the colours of 1 and 2 are most widely different, and it is estimated that about 
the worst case arises with monochromatic lights in the blue (470 my) and red 
(680 mu). The colour of light 3 is immaterial. Calculation for this case shows 
that 


(Ry + R,)¥¢ (G, + G,)¥” (B, + Ba) !® = (2-63R,)"% (2-63G,)"" (2-63G,)UP™. 
The breakdown in additivity corresponds to the difference between 2 and 2-63, or 
about 27%. As a second example, take a white of energy distribution W,dA 
appropriate to a black body of temperature 2060°K. and compare its small-step 


brightness with the sum of the small-step brightnesses of its component wave- 
lengths. Additivity would require 


; F , 1/e* fy? 1/8? 
{ RA GAY BAW dA = { | WARd { { W.G.ar| { | W,Bar} , 


but calculation shows that the equation is in error by about 8%. This second 
case is the one used by Ives to test the additivity of brightness as assessed in the 
flicker photometer. He found additivity held good within the experimental 
error (1 or 2%). No similar check of additivity for small-step brightness has 
been made, and it cannot be said with certainty that discrepancies of the kind 
indicated by the line-element would not be observed, although there is no positive 
evidence for them. The results of brightness matching by flicker photometry, 
like those obtained by direct-comparison photometry with full colour difference, 


cannot be derived from the line-element without additional assumptions and will 
not be discussed. 


$7. FECHNER_ FRACTION CURVES 


Although the modified line-element assumes that the red, green and blue 
mechanisms have different limiting Fechner fractions at high intensities, it still 
leads to a common value F for the limiting Fechner fraction of a pair of just- 
distinguishable and similarly-coloured light patches, whatever that colour may be. 
This is in accordance with Kénig and Brodhun’s classic measurements of the 
Fechner fractions of white and monochromatic light. These measurements 
were made with large patches, each 3° x 44°, and it is only at the higher intensities 
studied that the results can be accepted as referring to cone vision. The wide 
divergence in the observed values of the Fechner fraction at lower intensities 
is caused by the intrusion of the rods. The change with intensity in the Fechner 
fraction of monochromatic patches as predicted by the modified line-element (13) 
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is shown in figure 7. The Fechner fraction, plotted as ordinate, is the just- 
distinguishable difference in patch intensity divided -by the intensity of the 
feebler patch, and the intensity level, plotted as abscissa, is the intensity of the 
feebler patch. For A=555muy, this intensity is expressed in photons and for other 
wave-lengths in units such that all colours have the same brightness when com- 
pared on the basis of the calculated step-by-step visibility curve for the photo- 
metric level (figure 6). The Fechner fraction curves for different wave-lengths 
are not very different, but for the extreme blue the Fechner fraction has a rather 
smaller value at low intensities than for other colours. 


§8. HUE LIMEN 
The hue limen, or the difference of wave-length of two just-distinguishable 


monochromatic patches of equal brightness, is obtained from relations (13) and 
(18) by the substitutions, 


i W,R,, dbR=R,W,+ WS), Etc, 


On eliminating W;, the equation for the hue limen 8A takes the form : 


: G21 dR, Gee dG, C,2 dB, 2 (R 2 p. 1 
ec : “ <i rs 7 lar t 
(oy, H i FAN 4 eh Bea 5 oR) +two similar terms 


At high intensities this reduces to 


£6, pe 2 
F2=si2 a - a ae a8 a 1/p? + two similar terms J, 
D A A 


The equation used by Helmholtz is the special case of (24), when p, y, 8 are made 
equal. 

The values of the hue limen through the spectrum have been calculated from 
(23), assuming a common brightness for all colours on the basis of the calculated 
step-by-step visibility curve for the photometric level (figure 6) and an actual 
brightness of 80 photons for A=550 mp. 


‘Fable 1. Values of hue limen at the photometric level derived from the 
modified line element 


r OA X OA r OA 
(mp) (mp) (mz) — (mp) (mz) (mp) 
410 5-7 510 2:06 610 1-43 
420 3-9 520 2:74 620 1-88 
430 2-22 530 2-88 630 2-42 
440 1-75 540 2-87 640 3-2 
450 1-41. 550 2-62 650 4-7 
460 1-63 560 2-03 660 6:8 
470 1-14 570 1:72 670 11 
480 0-82 580 1-44 680 18 
490 0-97 590 1:27 690 29 
500 1-45 600 1-26 


The calculated values are plotted in figure 8 together with five experimental 
curves obtained by Wright and Pitt (1934) at an average intensity level (480 to 
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650 mp) of 70 photons. An outstanding difference between the theoretical and 
experimental curves is that in the former the limens in the blue are relatively low. 


9 
8 \ 
\ meas exam CALCULATED FROM 
N7 ---=D.H.2 
x + C. 8.0, Line ELEMENT (13) 
6 ab o-0 W.D.W. 
| » ee FH.G.R 
(4 
‘ Ud 
2 \ 
5 4 ‘ ee / 
= Zz Xx ‘ a ; 
2 
y 
zl y 
$o0 440~460~C480~CS00~C20 ~~ S40~«S~C«CS~*«~OS*C‘iSC*C«C EO 


WaVE-LENGTH (Imp). 


Figure 8. Observed (Wright and Pitt, 1934) and calculated 
values of the hue limen. 


In the determination of the experimental curves it was necessary to use lower 
intensities in the extreme blue, the total range being about 20 to 1. This may 
account for some of the difference. 


§9. GENERAL COLOUR LIMEN 


If two patches P and P’, not necessarily monochromatic, have the same 
brightness and are just distinguishable, they define a general colour limen. The 
general colour limens for a given common brightness level are conveniently 
displayed in a C.I.E. (x, y) chart (a) by marking pairs of just-distinguishable 
colours on various lines joining spectrum colours to other spectrum colours or 
to purples (Wright, 1941), or (b) by drawing small closed curves (ellipses) about 
various points in the chart, the colours corresponding to all points on each such 
curve being just distinguishable from the central colour. ‘The modified line- 
element will now be applied to calculate the g.c.l. (general colour limen) ellipses 
for the 25 colours used in the recent investigation by MacAdam (1942). Using 
an ingenious optical system, MacAdam arranged that, by turning a single control 
knob, the colour of a light patch was varied along a line in the C.I.E. chart while 
its brightness (C.I.E.) was automatically held constant. Repeated colour matches 
were made with a second juxtaposed patch of the same brightness and of a fixed 
colour on the givenline. The standard deviations of the settings were determined, 
a setting being specified by the distance in the chart between the variable and 
fixed points. ‘The standard deviation is closely related to, and may be taken as 
proportional to, the diameter of the g.c.l. ellipse at the fixed colour in the direction 
of the given line. MacAdam determined for one subject the complete ellipses 
for 25 fixed colours distributed fairly evenly over the domain between the spectrum - 
locus and the line of purples. 

Denoting by ¢,,,(m, n=1, 2, 3) the coefficients in the forms (15), the R 
co-ordinate of a patch of absolute energy distribution W,dA is given by 


R=tyx' +b" + t32', Ree Sn WS) 
where aos iz dx, Vee | Wiy dd, = [wzar, 
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and similar equations hold for G and B. The C.LE. unit co-ordinates GY, 
are then x=2'/w’, y=y'/w', x =2'/w', where w' =x’ +y’+2'. The quantity y’ is 
proportional to the C.I.E. brightness of the patch and, in fact, recalling the units 
in which W,d) is expressed, we have : 


patch brightness in photons =2-19 . 105 W,% da, 


where the mechanical equivalent of light has been taken as 0-0015 watts per 
lumen. MacAdam worked with a C.I.E. brightness of 15 millilamberts seen 
through a pupil of 2-6 mm. diameter, that is, with a brightness of 271 photons. 
The calculation will be applied, therefore, to the case, 

ZEA 


jE a) ty 104 TO 
[wasan yf awy= ze = 1-24. 10 


The & co-ordinate of the just-distinguishable patch is R +5R, where 
OR = 1,62’ + tyydy’ +4352", 


and similar equations hold for 5G and 6B. In MacAdam’s observations Sy’ =0. 
The differentials 5x’, 52’, expressed in terms of dx, 5y, are then : 


Thus, 
OR =x . (411 — ts) + By Fe (%ty3 — ty — ts), 


(ay en ee, ae re (26) 
y y 


and, inserting these expressions, the line-element (13) takes the form 


a = alts ch + {Pa acy} + {ea By} 


(2 Y 


2 G))\2 B))\? 
+ 28x dy [ Pup} + Pubes} + Pada | ] 
2 2 Z 
Ly dy? Res (Rb {ee @} Ee {Pe «ph | seme retsgsutihe (27) 
where 
a 7 (i ~ ts) pum sle ty3 — ty — tis)? 
pu= (tn — tog); Pro = 2 (x tpg — to1 — toa) 
ie o (t51 — ts2); P32= 7 (# t33 — ta1 — 33)» 


For a given colour «, y, the g.c.l. ellipse is completely determined by the equation 
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(27). The semi-axes and orientation of the ellipse have been calculated for 
MacAdam’s 25 colours and are as follows :-— 


Table 2 
Angle 
Colour between major 

MacAdam’s C.LE. Semi-axes axis and 

figure co-ordinates aieton oe 

number (degrees) 

3 y z Major Minor 

23 0-160 0-057 0-783 0-235 51052 O13 ORs Ome 85 
24 OF1S7e OLS SFO0r695 0-62 0-285 iD 
DS 0-253) 10-1250 7622 O27 0-272 53 
26 0-150 0-680 0-170 2°84 1-26 106 
27 0-131 0-521 0-348 2-60 1-04 101 
28 0-212 0-550 0-238 Desi lel Dy) 
29 0-258 0-450 0-292 1-90 0:97 86 
30 0-152 0-365 0:483 1-86 0-76 96 
31 0-280 0-385 0-335 1-67 0-85 78 
32 0-380 0-498 0-122 153 1-04 83 
33 0-160 0-200 0-640 1-10 0-44 92 
34 O222350225 050-522 1-24 0-57 76 
35 0-305, 0°323" 0-372 1-43 0-72 68 
36 (ests Wess) (evn 1-47 0-83 63 
37 0-472 0:399 0-129 23 0-78 54 
38 O51: 70-3505 O23 1-14 0-61 44 
39 0-475 0-300 0-225 1-27 0-52 43 
40 0-510 0-236 0-254 ieily/ 0-31 35 
41 0:596 0:°283 02121 1-02 0-35 33 
42 0-344 0-284 0-372 1-30 0-60 57 
43 039090237 50°373 1-18 0-44 45 
44 0-441 0-198 0-361 1-13 0-284 36 
45 0:278 0:223 0-499 1-14 0-49 62 
46 0-300 0-163 0-537 0-91 0-33 50 
47 0-365 0-153 0-482 0-99 0-24 39 


In figures 9 and 10 the calculated and MacAdam’s ellipses are plotted in the 
C.I.E. chart. There is some similarity. The orientations of the ellipses 
correspond rather well (figure 11), and their areas show a correlation (figure 12). 
The computed areas exceed MacAdam’s on the average by a factor of 20-7 
(4-5?), but MacAdam’s refer to standard deviations, which he estimates to be 
about one-third the corresponding limens. The axis ratios of the ellipses cover 
a range of about 2:5 to 1 and have average values 2:3 ( calculated) and 2-9 
(MacAdam’s), but they show no correlation (figure 13). 

MacAdam reports a limited number of observations for a second subject. 
The ratio of the g.c.l.s of the two subjects varies, for different points and direc- 
tions in the chart, over a range of about 0-8 to 2-5... For colour No. 35 sufficient 
results are given to construct the ellipse for the second subject (figure 14*), and 
in table 3 the characteristics of the ellipses of the two subjects are compared. 


* The plotted points in this figure were read off the smooth curves of MacAdam’s figures 8to 19, 
and there will be some reading error. 
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Figure 9, G.C.L. ellipses calculated from the line-element (13). 


Table 3 
Inclination of 
Major Minor Ratio Area of major axis to 
semi-axis semi-axis of axes ellipse axis of x 
(degrees) 
Principal subject OP 10 Otho. 10-4 veo O:7el Ome 70 
(P.G.N.) 
Second subject 0°31 0-165 1:88 16-1 44 


(D.L.M.) 
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In Wright’ 


On various line 
col 


S investigation (1941), the subjects determined equal colour steps 
sin the C.I.E. chart, and in figure 15 the logarithm of the observed 
our step is plotted against distance from one end of the line for four of the 
principal lines studied. Wright’s colour steps were not actual limens but 


Figure 10. MacAdam’s “ standard deviation ” ellipses for P.G.N, 


“equal small colour differences ” valie 
maintain a reasonably steady criterion. 

three times the limen. In reporting his 
four subjects to about the Same averag. 
factors : figure 15 shows the results befo 


found that after Practice a subject could 
On the average, the colour step was about 
results, Wright reduced the steps for the 


€ values by applying suitable constant 
re the application of these factors. 
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No direct calculations from the line-element were made for comparison with 
Wright’s results, but a fair idea of the variations along the several lines was 
obtained by interpolation from the 25 computed ellipses of figure 9. A similar 
derivation was made from MacAdam’s ellipses for P.G.N (figure 10). In figure 15 
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Figure 11. Inclination of the major axis of Figure 12. Log (area of ellipse). 
ellipse to the axis of x. 


the computed colour limens times 3 are shown as the thick-line graphs and 
MacAdam’s standard deviations times 9 as the broken-line graphs. These 
factors allow for the average systematic differences between small colour step 
and limen, limen and standard deviation, as estimated by Wright and MacAdam 
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respectively. Admitting that these factors may differ for different subjects, 
and that F may also differ, it is fair to suppose the four line-element curves in 
figure 15 all shifted up or down by the same amount when comparing with the 
experimental curves of each subject. Even then there is no close agreement with 


6 W. S. Stiles 


2 (Q) Line 460 70 570. 
7 


eo W.0.W ar 
x HV. W. (weianr 
@ R.A.L 
7S Milas W. 

@ P.G.N. MacAoam 
© FRom Line 


ELEMENT ( 13) 


(82). 


LOGARITHM OF THE COLOUR STEP 


“2 “4. 76 
DISTANCE ALONG THE LINE IN THE C.1.E.(X,¥) CHART. 


Figure 15. 


the curves of any of the five subjects. It may be noted that the experimental 
curves of different subjects, compared in the same way, show differences which 
are hardly less marked. 


§10. DISCUSSION 


The line-element may be regarded as a mathematical model of the mechanism 
underlying a particular group of visual properties (increment threshold, colour- 
matching relations, visibility curve, etc.). In the Helmholtz and modified 
elements, the respective terms correspond to three component mechanisms, with 
each of which is associated a unique spectral sensitivity curve. A difference 
between adjacent light patches will be perceived by the action of one of these 
component mechanisms, say the red, if the product, 


|dx| (x3 8) telmholz element) or |dR| (2) (modified element), 


exceeds Ff. The first factor in these products is the difference of stimulation 
of the red mechanism in the two patches, evaluated on the basis of the red spectral- 
sensitivity curve. ‘The second factor defines the sensitivity of the red mechanism 
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in the state to which it is brought by the total stimulation, x or R, which may be taken 
here as approximately the same in the two patches. The product, on this view, is 
a measure of the response of the red mechanism to the difference of stimulation. 
The only interaction between the red, green and blue component mechanisms 
contemplated by these line-elements is a certain degree of summation of their 

8x1/3 


responses, eS 


, etc., according to the sums of squares relation. 


The Schrédinger element (3) is an example of the wider class represented 
by 
6s°=[6x : Sal, y) =)? a [dy &y(%, ys %))? a [dz 5 EAX, y; z))? =f" 
Seaver (0) 


For line-elements in this class, the component mechanisms preserve their indi- 
viduality—with each is associated a unique spectral sensitivity curve *—but they 
are so linked that the state of each as indicated by the sensitivity factors g,, gy, g- 
depends on the stimulations of all three mechanisms. ‘This linkage is addi- 
tional to the summation of responses according to the sum of squares relation. 

For the above elements, the axes of the small ellipsoid defined by the line- 
element at each point of f.c. space are parallel to the co-ordinated axes. In the 
more general element, 


OSs ON One — Niet lols) B= cake (29) 


m,n 
the axes of the elementary ellipsoids at different points are not necessarily parallel 
to the axes of the ellipsoid at any one point in the space, but different rotations 
are in general required for different points. A rotation of the axes corresponds 
to a change of the primaries so that on this model the spectral-sensitivity curves 
of the three mechanisms are progressively modified as the stimulation (defined 
by the position in b.c. space) is varied. 

To what extent does the simple model defined by the modified line-element 
(13) fit the facts? The matter is complicated by the very considerable differences 
in the experimental results of different subjects. ‘To be satisfactory, the line- 
element (13) with a particular numerical determination of Rj, Gy, By, p, y, B, 
F and &(p) should reproduce the colour-matching values, the step-by-step 
visibility curve, the increment and colour limens and possibly some other visual 
properties of a particular eye, or, more precisely, of a particular rod-free area of a 
given retina. Changes in the visual properties in passing from one rod-free area 
to another in the same or a different eye should all be accounted for at a stroke 
by a change in the numerical determination of the line-element. For the com- 
parisons made in this paper, a single numerical determination has been used : 
for the distribution coefficients x, y, z the C.I.E. values, for the fundamental 
response curves the function defined in (15) for the limiting Fechner frac- 
tions numbers in the proportion 0-78 : 1 : 4-46, and for F the value 0-01. 
While certain main features of the experimental results are reproduced, 
satisfactory quantitative agreement with the measurements of all the visual 
properties for some particular subject has not been demonstrated. In fact no 
such complete set of measurements is available. For the restricted group of 


* Except in the degenerate case when any two of Bry Sy 82 stand in a constant ratio. 
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colour-limen measurements, the comparisons made (figures 8 to 15) show differ- 
ences which could no doubt be reduced by a different numerical determination 
of the line-element for each subject. But in the case of MacAdam’s subject 
P.G.N. it can be shown that complete agreement could not be reached by any 
change in which the fundamental response functions Ry, Gy, By remained always 
positive, linear forms in the C.I.E. distribution coefficients. (See Appendix.) 
Our conclusion must be that the modified line-element, constructed to suit 
measurements of increment limens, leads to the right kind of step-by-step 
visibility curve, but is in difficulties when applied to measurements of colour 
limens. 
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APPENDIX 


The Gaussian curvature of the surface of constant C.I.E. brightness according 
to the modified line-element 


Following Silberstein (1943), the unit C.I.E. co-ordinates x, y for colours of a 
given constant C.1.E. brightness can be regarded as curvilinear co-ordinates on 
a surface embedded in a three-dimensional Euclidean space, the distance apart 
of neighbouring points on this surface being defined by the expression for the 
line-element 6s. ‘This expression, for the modified line-element, is the square 
root of the right-hand side of equation (27). Silberstein uses an empirical 
expression for 6s derived from MacAdam’s measurements for the subject P.G.N. 
(figure 10), and shows that the Gaussian curvature K of the surface is not only 
not constant, but that it is positive in some areas, negative in others. 

The Gaussian curvature is invariant to transformations, linear or non-linear, 
of the co-ordinates x, y, and it is permissible to employ R, G in place of x, y. For 
the modified line-element, an expression for K is obtained quite easily if we 


accept for ¢(p) the approximation 


y er 
1+ 9p’ to which it reduces when p exceeds 1. 
This corresponds to using curve 7(p) instead of curve €(p) (figure 2). K is then 
given by 

Ke p’y’B(1+9R)(1 + 9G)(1 + 9B)v,v,v,(9y’ + ¥,+ Ug + V5) 
[oo 21 + OR E+ Pe + 9G + ROB 


where y’=v,R+uG+v,B=const. ‘This condition is obtained by solving 
the equation (25) for x’, y’, 2’ ; it expresses the fact that the colours considered 
all have a constant C.1.E. brightness (equal to 2-19 . 10° y’ photons). 

Since for real colours R, G and B are all positive, it is clear that the Gaussian 
curvature will have the same sign at all points of the surface of constant C.LE. 
brightness. 

Computation gives 


C= 2710 Ome 
Ug= 4:07 x 10-4 
Vy = —4-76 x 10-5 
UptUg+o,—= . 6-29 104 
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Thus v,vgv,(9y' +v,+vg+vs), and hence K, is negative for any value of y’, 
i.e. for any C.I.E. brightness. 

It follows that the modified line-element is inconsistent with MacAdam’s 
results for P.G.N., which for different colours yield K values of opposite sign. 
This inconsistency persists however p, y, 8 be chosen and whatever change be 
made in the fundamental response curves R,, G,, B; provided these remain linear 
forms in the C.J.E. distribution coefficients, and lead to all positive co-ordinates. 
R, G, B for real colours. 
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ABSTRACT. ‘The main purpose of the paper is to discuss the effects of errors of centring 
on the optical performance of plate-mirror systems. Such errors may occur during the 
aspheric grinding of the plates and also in the lining up of the system. It is first shown 
that if a grinding technique is used which builds up the asphericity at a proportionately 
equal rate all over the surface, the result of a large number of small centring errors at 
different stages of the grinding is very nearly equivalent to a simple decentring of the 
asphericity on the surface, combined with a small amount of primary astigmatism. Next, 
general formulae are obtained for the effects on the Seidel errors of a centred system 
of decentring and tilting its components, and it is shown that in the case of a plate-mirror 
system the formulae can be used to estimate the practical tolerances for disturbances of 
this kind. Lastly, the general formulae are applied to discuss the lining up of a two-sphere 


one-plate Schmidt-Cassegrain camera. 


§1. EFFECT OF ERRORS IN GRINDING 


nN the construction and use of plate-mirror systems, problems of centring 
are met-with at two stages. The first is during the actual aspherizing of 
the surfaces ; the second during the lining up of the system. An error in 
centring at the first stage will in general lead to an asymmetrical error in the 
figure of the surface. ‘This will result in a slight change in the centring chosen 
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for the surface when the complete system comes to be lined up, for the centring 
chosen will be that at which the optical effect of the asymmetry appears as small 
as possible. Thus the practical effect of the asymmetry depends not so much 
on its actual amount as on the extent to which it cannot be compensated by 
lateral shift of the surface which bears it. 

Broadly speaking, the most successful methods of figuring aspheric plates 
are those in which the asphericity is built up at a proportionately equal rate all 
over the surface. One such method, well suited to the resources of the amateur 
astronomer, has been described in an ably written paper by H. W. Cox (1939). 
It depends on the use of a flexible grinding lap consisting of lead facets arranged 
in petal-shaped areas on a Sorbo rubber backing, operated by a mechanical 
traverse. ‘The lap is traversed in short, fairly rapid strokes as the plate slowly 
rotates. With such an arrangement, asymmetry arising from random variations 
in the stroke is eliminated for all practical purposes, and inaccuracy of centring 
of the plate in its holder, or of the plate-holder on its spindle, becomes a potentially 
important source of error. 

To estimate the practical effect of errors. arising in this way, we observe that if 
the grinder is adjusted to generate a true curve, with only small relative errors, whose 
depth increases at proportionately equal rates all over the surface, then the effect 
of having the plate wrongly centred for part of the time is to contribute a weak 
‘“‘ Schmidt plate”’ with displaced centre to the figure instead of adding its strength 
to the main (i.e. truly centred) figuring. 

Suppose that m errors of this type are made and that the plate is centred by 
trial in the final assembly. ‘The plate figure is then the sum of m figurings, each 
with aslight decentring, whose strengths together make up the full plate-strength B. 
We can investigate the properties of such a plate by applying the methods of 
Conrady (1919). 

The deviation (Ax, Ay) of the point of intersection with the Petzval surface 
of the ray through the point P with polar co-ordinates (S, £) in the entry pupil 
is given in the case of pure primary spherical aberration by 


Ay BS? cose, Ay = BS sin E, 


where the coefficient B measures the amount of the aberration. From adecentred 
spherical aberration of centre (€,, 7,.) and amount B, we therefore have deviation- 
contributions 

An= BS 20s k- Ay=B,S,3 sin E,, 


where (see figure 1) 
COs LE, =Sicosh cy S,sin£,=Ssin E—», ; 
S77 = S?—2S(E,.cos E+, sin E) + (€,2 +7,2). 


If €,, 7, are so small that their squares can be neglected, the deviation-contributions 
can be written 


Ax = B,[S? + (£2 +,") —2S(€, cos E +7, sin E)|(S cos E— é,) 
= B,.S° cos E — B,€,,S?(2 + cos 2E) — B,n,.S? sin 2E, 
Ay = B,S* sin E — B,€,.S? sin 2E — B,n,.S*(2 —cos 2E) 


oe 
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and the total deviations effected by the plate are accordingly 


Ave ( E B,) S® cos E — ( > B,e,) S?(2 + cos 2E) — ( s Ban) S? sin 28, 


r= r=1 


r= 


Ay = ( > B,) S8 sin EB — ( 3 Re.) S?sin2E — ( : B,n,) S2(2 —cos 2B), 
: r=1 


r=1 | 


Wn 


i.e. are those corresponding to primary spherical aberratioa of amount B= © B 
r=1 


together with a coma contribution which is measured by the sum of the vectors 
B,(é,, 7,)- ‘This coma will be annulled if the plate is given a lateral displacement 
(=, H), where 


Yr? 


mm 1 

es fade ne Soyer, o> Tee (2) 
r= r=1 
when the system is assembled. Since the plate will in practice be centred by 
trial, such a displacement will be given automatically. The largest errors then 


Figure 1. 


remaining are the second-order terms (neglected above), and, as will appear later, 
these give rise to astigmatism, in second-order small quantity, uniform over the 
image-field. 

In certain important cases* this astigmatism also can be compensated by 
adjustments of the components perpendicular to the axis ; it follows that in 
these cases a first-order centring error in aspherical grinding, provided it 1s of the 
type specified above, leads only to errors of the third order of small quantities in 
the finally adjusted system, i.e. to extra terms in the expressions for the displace- 
ments Ax, Ay in the final system whose coefficients are third-order small quantities. 
Such terms have no perceptible effect on the physical image in many practical 
cases. In general, a second-order amount of astigmatism remains after the coma 
has been removed by empirical plate-adjustment. 

The general conclusion is that, provided the grinding is arranged so as to 
generate the asphericity at proportionately equal rates all over the surface, the 


* For example, in the monocentric Schmidt-Cassegrain. systems and in the two-plate Schmidt- 


_Cassegrains. 
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effect of the centring errors, which must always be present in some degree, is far 
less serious than might appear at first sight. It seems likely that this circumstance 
plays a large part in determining the success or failure of a grinding technique. 


§ 2. EMPIRICAL LINING-UP OF A CENTRED OPTICAL SYSTEM 


Consider the ‘‘ equivalent plate system”’.* Let B,, B,....B, be the plate- 
strengths ; 0), o,....0, the respective distances of the (real or virtual) plates 
in front of the entry pupil ; (&, 7), (€ 72), ----(En» Mn) their respective 
lateral displacements. For the present, we suppose the components of the system 
to be free from errors of figure. 

Then, as before, the displacement-contribution from the first plate to the 
intersection with the Petzval surface of a ray, originally parallel to the axis, 
through the point of the entry pupil whose polar co-ordinates are (S, £), is 


Ax = BiS,* cos E,, Ay = B,S,? sin E,, 
where 
S, cos LE, =Scos F—&,, S,sin £, =Ssin E —y,, 
Sy? = S? + (€? +77) —2S8(& cos E +7, sin £). 
Thus 


Ax = B,[S? + (€)? +717) —2S(€, cos H+, sin Z)|[S cos E — &] 
= B,[S* cos E — S*(E,(2 + cos 2E) +7, sin 2E) 
+S((38)? + 7,7) cos B + 2& yn sin E) —& (1? + 3”)], 
Ay = B,[S° sin E — S*(€, sin 2E + y,(2 —cos 2E)) 
+ S((? + 3m?) sin B+ 2&n, cos E)—14(:? +71”)]- 
If we write - 
fi +i1=pie™, 
these equations can be given the more compact form 
Ax +iAy = Bye'?’[See2— $2) — $25, (2 + e2i(E—$1)) 
+ Spy2el®-ON(2 4+ e293) — p, 3], 


For the complete plate-system we therefore have, in Seidel approximation, 


T= 


n Fi * 
Sil r=1 7T=1 


n n n ; 
+ § (20 = Bp, +e—% E B,p ete) 2 2 Bp, Rompe = (3) 
t= r= t= 


where B, denotes the strength of the rth plate and p,e#r=£,+7n, its lateral 
displacement. Here Ax, Ay measure the lateral displacements in the Petzval 
surtace andB,. Bae neue measure the plate-strengths, both in suitable units. 
The connection between the two sets of units can be expressed by saying that if for 
the unit of plate-strength we choose the strength of the Schmidt plate of a concave 
spherical mirror of paraxial focal length f,, then the actual linear displacements — 
5x, Sy in the Petzval surface are obtained on multiplying Ax, Ay by f/8f,3, where 
f is the focal length of the system. : 


* In the sense of the method of plate-diagram analysis. See C. R. Burch, Mon. Not. R. Astr. 
Soc. 102, 159-165 (1942) ; E. H. Linfoot, ibid. 104, 48-64 (1944), for an account of this method. 
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The first three terms on the right-hand side of equation (3) are easily recognized, 
on separating them into real and imaginary parts once more, as arising respectively 
from spherical aberration, coma, and astigmatism combined with change of 
focus. ‘The fourth term, being independent of E and _S, represents a lateral shift 
of the image as a whole. 


nr 
We shall suppose the system designed so that © B,=0, i.e. spherically cor- 
- r=1 
rected. ‘Then the conditions for the axial image to be error-free are, if we do 
not permit a change of focus, 


mB pe =). pape or), 2 Dep Us 0 ee ene (4) 

r=1 r=1 r= 

In practice, however, the focus is usually determined empirically and the term 

2Se” x B,p,2 will not show itself as a defect in the image. It follows that the 
f 


conditions to be satisfied are 


n n nw 
2B 0, » B,p,€%r =0, ae bi Pigs =A Oe eee (5) 
r=] r=1 r=1 
We can express this result as follows : Jf the (real or virtual) elements of the plate- 
system associated with an optical system are given lateral displacements v,, v2 .. 
(v a vector) satisfying the conditions 


HPL, Dye She he ee (6) 


the axial image will not develop spherical aberration, coma or astigmatism. It 
will be laterally displaced by an amount proportional to XB,|v,|?v, and will be 


moved along the axis, away from the system, through a distance proportional to 


2>B,|v,?|. 


Equation (3) gives a simple proof of the property, mentioned in §1, that 
when a plate suffering from errors of figure of the type there described is centred 
empirically, the residual error is effectively primary astigmatism, in second- 
order small quantity. It is only necessary to set o,=0,.... 0, in the present 
discussion (¢,,+3, ---- %, remaining arbitrary) and to observe that the operation 
of centring the plate empirically is expressed by adding to each of p,e™, .. 
Pmé?m the complex constant 


n ; Me 
c= » B,p,€?r 3, Des 
pt pasil 


After this operation, the coefficient of S? in (3) vanishes and the third: term 
becomes the leading term. 


Off-axis images 
To determine the effect on the off-axis errors of a system of decentring its 
plates, we consider a parallel beam coming in at an angle V with the axis, the 


azimuth of this angle being y. The principal ray meets a plane distant o in fron 
of the entry pupil in the point 


X =oV cos x, Y-ol sim y, 
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and the effect of the plate is, therefore, the same (in Seidel approximation) as if it 
were situated in the entry-pupil but with a lateral, displacement (—oV cos x, 
—oVsiny), and the pencil were parallel to the axis.* It follows that the off- 
axis errors for laterally displaced plates are obtained on replacing (€,, 7,) by 
(€,—o,V cos x, n,—0,V sin x) in (3). Set 


E,—0,V cosy=P,cos®,, n,—0,V siny =P, sin®,, 


which gives 


Pe*+**r _ p,ex'?r —o, Vex, 


and we obtain in place of (3) the equation 


Ax +iAy= See 3 B, - s+(23 BP +e 3 B,Pier) 


fa r=1 


Ui nr n 

+ S (20 >» BP Pe B, Pte) + 2B, BP i ana 
r= r=1 pil 

The four terms of this new equation have the same significance as before, and 

they now specify the spherical aberration, coma, astigmatism with change of 

focus, and lateral displacement} at any point of the image-field of the system 

with decentred plates. 

We assumed above that the optical components were free from errors of 
figure. ‘The same analysis suffices, however, to cover the more general case 
where the plates suffer from errors of figure of the type discussed in § 1. It is 
only necessary to represent each plate as a cluster of mutually decentred, error- 
free component-plates held rigidly in contact, and to apply (7) to the system of 
component-plates so obtained. ~ 

An inspection of the terms on the right of (7) enables us to draw one or two 
general conclusions about the effect of these errors of figure on the performance 
of an empirically lined up system. 

Since the first term S%e” 25, , which represents the spherical aberration, does 


not involve V or the AEG os Pps Py, it follows that the spherical aberration of the 
system is unaffected by the errors and remains constant all over the field. 
The second term, which represents the coma, falls into two parts : 


— S22DB, er + 2E DBP ein) 
-- S2(20B, per + e#2 SB peor) 
r Up 
+ SV (2 + ER) EB io eee, (8) 


The second part represents off-axis coma of the usual kind, and vanishes identically 
if, as we suppose, the system is designed to satisfy the aplanatism condition 
2 B,o,=0. The first part does not depend on V or x, and so represents a coma 
; 


uniform over the whole field, which can be eliminated by a lateral shift of any one 


* The effect of the small relative tilt between plate and off-axis pencil is zero in Seidel 
approximation. 

} The last term is associated with the distortion of the system, but does not by itself measure 
the distortion. See Mon. Not. R. Astr. Soc. 104, 49 (1944). 
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of the plates in the plate-system, or of several plates together, so as to secure the 
fulfilment of the condition 2B,P. 6 =0. 
; 


The third term, 
S(2e? XB, P+ e— TB, P2e%?r) 


= S(2e XB,” +e 7k =B,p,2e%r) 
y 7 


r 


—2SV(2e" = B,p,0, cos ($, — x) +e UB, p,.0,e4?r +X) 
r n 


SV 2 2eior eX IB) Big Fo (9) 


represents astigmatism and focus shift.* It falls into three parts, in accordance 
with (9) ; the third part contains the ordinary off-axis astigmatism of the system, 
which vanishes in systems (the so-called Seidel anastigmats) in which ©B,o,2= 0. 


; 
The remaining parts represent the effect of the decentring ; we shall consider 
them in more detail, in the next section, in the particular case of a three-plate 
system. 
Finally, the fourth term, 


DIY be eae = UB,P,e?| Pe [2 
r 7A 
= wb e ae VIB, p,20,(2e'% + en) 
r : 
ar V?XB,p,0,2(2e%r ae e2*X— tbr) 
e 
OOS Bee eee ee (10) 
ip 


- which represents the plate-contribution to the distortion of the system, remains 
independent of S and F when the plates are decentred, and so does not affect 
the sharpness of the images. But it does depend on V and on the quantities 
P+» Py, SO that, as was pointed out by Conrady in the paper already referred to, a 
system lined up to give satisfactory images may still contain centring errors 
which appreciably modify its distortion. 

The above analysis applies in principle to ordinary lens systems with spherical 
surfaces ; their associated plate systems, to which we can apply (7), consist 
entirely of virtual plates. But the application is not very satisfactory in the case 
of wide-aperture lens systems, since in these the Seidel theory is far from giving 
a true picture of the image-errors, especially in the outer parts of the field. In 
its application to plate-mirror systems of the Schmidt-Cassegrain type, however, 
the theory supplies information of practical value even at nominal apertures as 
highas F/2. This different state of affairs can be attributed to the much smaller 
size of the fifth-order error-coefficients which results when mirrors are used in 
place of lenses as optical power-surfaces, together with the circumstance that 
vignetting restricts the field diameter of an F'/2 Schmidt-Cassegrain to something 


less than 15°. 


* An incidental advantage of the complex notation used in this section is that the terms 
representing focus-shift appear separated from those representing astigmatism. The former are 
the terms in Set/, the latter those in Se—?£, 
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§3. APPLICATION TO SCHMIDT-CASSEGRAIN CAMERAS 


As a practical example, we consider the lining-up of a Schmidt-Cassegrain 
camera consisting of two spherical mirrors and one aspheric plate. The general 
appearance of such a system is shown in figure 2. 

The optical lining-up of a correctly figured system of this type can be carried 
out as follows : : 

(a) place the primary mirror M, symmetrically in its mounting by measure- 
ment ; 

(b) place the plate P at the measured distance from M, and centre it by 
measurement ; 

(c) move the secondary mirror M, along the axis, keeping it roughly 
centred, until spherical aberration is zero, when P is moved laterally 
so as to annul coma, and thus allow the spherical aberration to be 
observed ; 

(d) tilt M,, readjusting P laterally each time, until the astigmatism is zero 
on axis ; ‘ 

(e) move P along the axis until the coma remains zero for small image-dis- 
placements off axis. 


ey meee re ee meet ee es ee 


Se 


: m1 —> 
Figure 2. 


The system is then in “ Seidel adjustment’. This adjustment must now 
be altered, by an amount which depends on the field to be covered, in order to 
balance the higher errors by suitable amounts of opposing Seidel error over the 
required field. 

In the monocentric Schmidt-Cassegrains, the change consists in moving M, 
and P along the axis to balance the lateral spherical aberration against primary 
spherical aberration, while keeping the coma zero. This gives a system which 
is extremely close to a monocentric system. The last operation is always to 
adjust P laterally to bring the coma to zero on axis. In the flat-fielded aplanats, — 
the compromises to be made are of a more complicated nature. 

In a practical case, however, the plate-figure usually contains a little primary 
astigmatism arising (in the manner already described) from centring errors 
during grinding. ‘The amount of this astigmatism can be kept within tolerable 
limits by the use of proper precautions during the aspherizing process. For 
example, the plate of an F/2 Schmidt-Cassegrain system might be regarded as 
satisfactory if the astigmatic confusion-circle corresponding to its astigmatism 
of figure did not exceed 1/559) inch in diameter. But the presence of this astig- 
matism complicates the task of lining up the system ; if the above procedure be 
adopted it results in a decentred line-up which Seidel theory and experiment 
agree in condemning as unsatisfactory in most cases because of its effect on the 
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images in the outer parts of the field.* It is therefore better to begin by placing 
the two mirrors M,, M, in correct alignment by mechanical means (with pro- 
vision for varying the position of M, along the axis) and to confine the empirical 
adjustments to operations (c), (e) and the annulling of coma by lateral plate- 
shift. 

The question then arises, how accurate the mechanical alignment of M, and 
M, needs to be. We shall investigate this by applying the general formulae of 
the last section. To be of practical value, the investigation must take account 
of the higher aberrations as well as the Seidel properties of the system. In 
the case of wide-aperture lens systems, where the higher aberrations have to be 
carefully compensated over a relatively large field, a general theoretical discussion 
of centring tolerances would seem to be hardly practicable at the present time, 
since it would at once involve us in the formidable analytical complications of 
fifth-order optics. We should need to rely instead upon conclusions drawn 
from experiment, and these conclusions would only remain valid for systems 
very close in design to those on which the experimental observations had been 
made. In the case of the Schmidt-Cassegrain systems, comprising only three 
or four effective optical surfaces and working over a smaller field, high definition 
is a consequence of the intrinsic smallness of the higher aberrations of each 
element rather than of their mutual compensation.t A small error in the lining 
up of the system therefore results in only a small percentage change in the amounts 
of its higher aberrations. Since the Seidel coma and the primary spherical 
aberration are under continuous control, t it is the size of these higher aberrations 
which determines the actual confusion-spreads in a Seidel anastigmat (or near- 
anastigmat) in correct adjustment. It follows that small centring errors will not 
markedly affect the confusion-spreads of a Schmidt-Cassegrain system pro- 
vided that their contribution to the Seidel errors is small compared with the 
image-spreads of the correctly adjusted system. A more convenient formulation, 
practically equivalent to the first, is that the centring errors can be tolerated if 
their Seidel contribution to the displacement Ax+zAy is everywhere small 
compared with the allowed upper limit to the size of this quantity for the system. 

We can obtain a reasonable estimate of the amount of decentring which can 
be tolerated in the line-up of the secondary mirror of a Schmidt-Cassegrain 
camera by laying down that the resulting final disturbance 5x +75y of the Seidel 
displacement-vector in the Petzval surface must nowhere exceed one fifth of 
the allowed image-spread for the system. Ina practical case, the allowed image- 
spread might be 0-001 inch, and the tolerance would then be given by the condition 

[5x -4iBy|-<zaqm inch 
over the whole field. 

* The monocentric Schmidt-Cassegrains are to a certain extent an exception to this rule. 

+ Higher spherical aberration (zonal aberration) appears as an exception to this statement, 
inasmuch as it is usually fully compensated in the plate-figuring. But this aberration also is 
intrinsically very small, as can be shown in a striking manner in a particular case by removing 
the secondary mirror of an F/2 two-sphere aplanat, placing the plate in contact with the primary 
and Foucault-testing the pair with a pinhole at an adjustable distance along the axis. When the 
pinhole-distance is suitably chosen, the residual spherical aberration is found to be scarcely 


detectable with the knife-edge. 
t They can be varied at will by altering the spacings of the components. 
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The disturbance of the Seidel properties due to centring errors in the line-up 
can be calculated from the general formulae of the last section as soon as the errors 
have been expressed in terms of displacements of the elements in the associated 
plate system. ‘To do this, we image the surface-figurings, aspheric plates and 
missing Schmidt plates of the system through into object-space in the usual 
way, taking account of the errors in tilt and in lateral adjustment of the mirrors 
when calculating the images according to Gaussian optics. ‘The associated 
plate system thus obtained has its elements laterally displaced and tilted. The 
tilts of the elements have no effect on the Seidel aberrations and are therefore 
disregarded ; their lateral displacements are a set of parameters which describe, 
sufficiently fully for the present purpose, the errors in the line-up of the original 
system. 

A Schmidt-Cassegrain with both mirrors spherical is a ‘‘ three plate system” ; 
that is to say, its plate-diagram consists of three plates and an aperture stop. 
Two of these plates are virtual ; they are the missing anastigmatizing plate P, 
of the primary M,, and the image P, in M, of the missing anastigmatizing plate 
of M, (both with strengths reversed, since the plates are missing). ‘The third plate 
P, is the glass aspheric plate P of the system. 

In the aplanats, the plate-strengths B,, B,, B, and the distances o,, o2, o3 of 
the plates in front of the aperture-stop satisfy the relations * 


B,+8,4+ B,;=0, Bo, + Bao, + B30, =0, Got boc (11) 
while the off-axis astigmatism is proportional to 
Byo\?-+.Boae i Bea... ane) eee (12) 


Taking M, as defining the axis of the system by means of the normal through 
its pole, we set p,=0 in the general formulae (8), (9), which then become, for an 
off-axis beam of azimuth y, 


Coma : 

— S*[2Bpp,e%?: + 2Bspze"*s + e*"( Bop pe"? + Bapse—#2)]. ss. s, (13) 
Astigmatism and focus-shift : 

S[2e"( Bop,” + Bsps”) + e—”(Bopa2erbe + Bsp3°e""?s)] 
—2SV[2e'"{ Bopoa9 cos (2 — x) + Bgp3a3 Cos (¢3—x)} 
+e“ 1X! Bopyone"?: + Bapzoze*}] 
ES PAF a eth) Bog Boa, oe) ee (14) 

as already remarked, the terms Se represent the focus-shift and those in Se~' 


the astigmatism. 


In the lining-up procedure described above, the decentring (ps, d3) given to 
the plate P, is determined by the condition that the coma shall vanish, i.e. by the 
equation 

Bopoe'* + Bapzes =0, 
* The notation here and below is that of a previous paper (Mon. Not. R. Astr. Soc. 104, 48 
(1944)). Since the unit of plate-strength is there defined as the strength of the Schmidt plate 


of M,, it follows that expressions (13) and (14) below must be multiplied by //8f,? to obtain their 
contributions to the actual displacements in the Petzval surface, and that the value of B, is —1. 
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which is equivalent, since B, and B, are both positive, to the relations 


$3=7+ do, B55, —4s.0,: eer Ch) 
Expression (14) then reduces to 


; ot 
S(2e% + ribet) e pe 
3 


— 2SV (a =-o3)[2e cos (6, — x) + eh 92+ X) By, 
+ SV%(2e™ + eX) Bip, + Boo,” + Bsps?), 
in which B, = —1, B,+B,=1. 

If the plate and mirrors are free from errors of figure, the reduction of the 
observed astigmatism in, the axial image to zero by tilting M, (operation (d)) is 
equivalent to setting p, =0 in (16), i.e. to the correct alignment of M, relative to 
M,. It then follows from (15) that p;=0 and the system is correctly centred. 

The practically important case, however, is that in which the secondary 
mirror M, is slightly out of alignment in relation to the primary ; that is to say, 
the virtual plate P, is decentred and the consequent decentring of P, is given 
by (15). The decentring of P, may represent the result of applying operation (d) 
to a system whose plate is not quite free from astigmatism of figure, or it may 
represent an inaccuracy in the mechanical positioning of M, relative to M,. In 
either case, we only need to consider lateral shifts of P,, since the longitudinal 
adjustment of M, is made empirically.* 

To fix ideas, we consider the special case of the two-sphere one-plate astro- 
graphic camera 


G0 30) A=; I=), DP =0-7399, oa 21447. vase ek (17) 


constructed with f, =f,=20 inches. We suppose that the aperture stop, 10 inches 
in diameter, is situated 20 inches in front of the pole of the primary mirror and 
that the system, of focal length approximately 28 inches, covers a field of 6° 
diameter at a nominal aperture of F'/3. Then inthe plate-diagram of the system 


Bi=-1, B,=@(2—¢?=0-2601, B,=l=0-7399, 
Sag) 


—1=0-5882, Geo — | = 114475 0 were (18) 


and the unit of length f, = 20 inches. 
To express the effect of decentring P, in the azimuth $,=0 and readjusting P., 


to annul comaall over the field, we set 
By 
db: =9, Ps San. 
Expression (16) then becomes 
S(2e® + e-®) x 03515 p,2 + SV(2e cos x + e+) x 02895 py 
+ SV2(2e + eX-iF) x 00596, 
* An error in the longitudinal adjustment of M, would cause a change in the strength of P, as 


well as in its position. . 
+ The stop is placed about half-way between the plate and the primary mirror in order to 


obtain uniform lighting of the field with the least possible increase in the diameters of plate and 
primary over that of the stop. 


p= 0-3515p,, pg=7. 


ww 
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since 21-9) 
B,o,2-+ Boo? + Byos? = a =0-0596. 
3 
The actual displacement-contribution in the Petvzal surface for the ray through 
the point (S, £) of the aperture stop is therefore, 


Sx +i3y = on (Ax + iAy) 
1 


=. a a [22+ e-#) & 0-3515 ps7 /fz7 


+ (2e” cos x + e~ +X) x 0:2895 Vps/f, 
+ (2 + eX) x 00596.) 7], | eee (20) 
where V denotes the off-axis angle of the incident pencil and x its azimuth. 


In the case py =0 of a correctly lined-up system, the contribution of the off-axis 
astigmatism to the displacement is 


(2e# + e'x—-*#) x (0596 SV?2/8(1 —q). 
Since the absolute value of this expression never exceeds 


3 x 0:06 x 5 a: f 10 _ 0.0004 i h 
x ex 0 CG Pa inch, 


we see that the system is for all practical purposes anastigmatic. 
In the case py>0 we have at the edge of the field: 


(Change in 6x+726y due to decentring) 


= a pleet re) x 0-00072 py 
+ (2e# cos y + e~#+*X) x 0-00072 pa], 


where ps, S are measured in inches. Since S<5, the absolute value of this 
expression never exceeds 


0-001(2:35 ps? + 1-94 py) inch. 


A decentring in which p, does not exceed 0-1 therefore disturbs the contribution 
from astigmatism and focus-shift to the Seidel displacement 5x +78y in the Petzval 
surface by an amount which nowhere exceeds 0-00022 inch in absolute value. 
Its effect on the image-quality will therefore be harmlessly small. 

There is, however, another consequence of the decentring which needs 
to be considered, namely a disturbance in the distortion of the system, called 
“dislocation of the image” by Conrady. Only if this also is harmlessly small can 
the decentring be tolerated in a system whose images are to be used for measure- 
ments of position. ‘To settle this point, we use expression (10), which in the 
case ¢2=0, ¢5=7, Byp.=B3p3 becomes 


(Bop2° — Bsps°) — V(Bop,?o + Bzps?03)(e% + 2 cos x) 
+ V°(B3p20" — Bsp303”)(2e'X cos x + 1) — V3eX UB,o,3. 


The first term merely represents a lateral displacement of the whole field, while 
the last term is the ordinary plate-contribution to the Seidel distortion of the 


On decentred aspheric plates 77 


system. The effect of the decentring is represented by the terms in V and V2, 
whose contribution to the actual displacement 5x +73y in the Petzval surface is 


f B 
SF — VByp.?| o.+ B, os) (@ +2cos x) 
+ V?Bop(a9?—03")(2e'X cos x + 1) | 


1 ; 7 
= 3940 [—0-2576 Vp,?(2e*x + e-*X) —0-2508 V2p,(2e + 1)]. 


: 1 
Since V<59> the absolute value in inches of this contribution never exceeds 
(17-55? + pg) 10-, 


which for p,;<0-1 is utterly negligible. 

It follows that a decentring for which p,<0-1 causes no appreciable dislocation 
of the image, nor any serious deterioration in the image-quality. Centring of 
this accuracy is readily obtainable by mechanical methods.* 

Finally, we remark that the compensation of an axial astigmatism 


dx + 1dy =4Se— x 0-00022 . 
due to error in the plate-figure can be shown by the same methods to involve a 
decentring, of approximately five times the above tolerance limit, which introduces 
astigmatism near the edge of the field several times greater than that removed 


at the centre. The use of an intentionally decentred line-up is therefore not to be 
recommended. 
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ABSTRACT. Equations are derived expressing the frequency of vibration of rectangular 
wood and plywood plates in terms of elastic constants, dimensions and density. For 
supported edges, the solution is exact and gives the complete series of overtones. For 
clamped edges, the fundamental frequency is estimated by the approximate Rayleigh 
method, but evidence is presented which indicates that the errors introduced by the use 
of the Rayleigh method are, in fact, negligible. A method of measuring the frequencies 
is described; comparison of calculated with observed frequencies shows that the former 
are on the average about 23% high for clamped edges and about the same amount low 
for supported edges. ‘This discrepancy is attributed to lack of correspondence between 
the experimental edge’ conditions and those assumed in the theoretical treatment. Finally, 
an empirical method, based on the rough proportionality between elastic constants and 
density of wood, is suggested and justified for obtaining approximate estimates of the 
fundamental frequencies in terms of plate dimensions only. ; 


INTRODUCTION 


ONSIDERABLE advances.have recently been made in the elastic theory 
of wood and plywood. ‘This theory has been used to calculate stresses 
and deflections under lateral loading and critical buckling stresses of 
plywood plates ; the agreement with experiment is reasonable though not exact. 
The frequency of vibration can be obtained theoretically by closely allied methods, 
and one of the objectives of the present investigation is to extend the comparison 
between theory and experiment to this frequency. A second objective is to 


gain some information regarding the frequency, which does not appear to have 
been studied hitherto. 


§1. DEFINITIONS AND BASIC EQUATIONS 


As is usual in problems of this type, it is assumed that the material is homo- 
geneous and perfectly elastic and that, in the two-dimensional case, its elastic 
behaviour is determined by the existence of two perpendicular symmetry axes. 
For the particular purpose of the present investigation, it is further assumed 
that the symmetry axes (i.e. the directions parallel and perpendicular to the grain) 
are parallel with the sides of the plate. The dimensions of the plate are a, 5 
and h, where 

a=length of side parallel with Ox. 
Da gry Sees be oo OE 
h=thickness of plate, parallel with Oz. 
(A is small in comparison with a and 5.) 
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Under these conditions the differential equation governing the small 
deflections w in the Oz direction is (March, 1936, 1942) 


O*w O*w O4*w 
Diva se Oipre eb? farce ad oa Retevenrie (Eek) 
where 
VOR 
opie 
ais oe 
D= SiON , \ a, Oeeters (122) 
oh. E\95,h3 Gi 
De iN 6 
= 1), +2D;-: 


E, = Young’s modulus in bending in Ox direction ; 


> 
E,= ” »? ” » ” Oy 


” b] 


G,.= Rigidity modulus for shear stresses in Ox and Oy directions ; 


> 
contraction in Ox direction 


on tenslOnin Oy direction 
extension in Oy direction 


A= 1 —0530,.~ 0-99 for wood and plywood ; 
P=pressure acting in Oz direction. 


The differential equation governing the vibration of the plate can be found by 
substituting ph(d?w/dt?) for — P in (1.1) : 


O4w o*w 07w 


0*w 
EE ee ord Pees WE Drs ere + phases =i il} a aateree (1:3) 


where p = density of material (weight at test/volume at test) and ¢= time. 


The boundary conditions for a plate clamped at the edges are (Timoshenko, 
1940, § 22): 


w=( at:— eh RSG SS SS 

Ow 

IMG te eee We ae (1.4) 
Ow 

Arte mae =0, =b; 

ae 0 at y y 


while for a plate simply supported at the edges they are (Timoshenko, 1940, 
§§ 22, 26) : 


w=0 at. x=), x= a, y=0, ao | 
2 
oe =0 at :— Oe aes 
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The potential energy of bending (V) is (March, 1942 ; Goland, 1942) : 
: 0?w aw\? 0?w 0? w O?w } 
v=5[ [4p (F) +D,(54) +2D 5m aa + 4Ds (<a) dxdy, 


and the kinetic energy T'is 


§2. DERIVATION OF THE FREQUENCY EQUATIONS 


In the case of supported edges, an exact solution can be obtained by expressing 
w as 


NMmx . nary . (2 1) 


wW=W,) SIN - See Caer 54s 


where p=2z x frequency (v), m and n are integers, gnd wp» is a constant. Sub- 
stitution of (2.1) in (1.5) shows that the boundary conditions are satisfied, and 
substitution of (2.1) in the differential equation (1.3) leads directly to 


7 


v {5 [D,m'*b4 + Don*a4 + Dyin} eer ee: (2.2) 
where the suffix s indicates supported edges. 

The fundamental is given by m=n=1, and the overtones by higher values of 
m and n. 

In view of the complexity of the problem involved in a complete solution, the 
fundamental frequency for clamped edges has been derived by the Rayleigh 
method (see e.g. Timoshenko, 1937, § 16), supplemented to some extent by the 
use of the more accurate Ritz modification. For comparison both with the | 
exact solution (2.2) and with the results for clamped edges, the Rayleigh and 
Ritz methods have also been applied to the case of simply supported edges. 

The Rayleigh method consists in working out the potential energy V and the 
kinetic energy 7 from the equation 


w= f(x, y)g(t),, ee eee (2.3) 


which expresses the deflection w as a function of time ¢ and of position on the 
plate. If the system is conservative, then 


Vas = Teas ware ene (2.4) 


and the frequency of vibration is found from this condition. 

The virtue of the Rayleigh method is that in general it is possible to use an 
approximate expression of f(x, y) in (2.3) and still obtain estimates of the fre- 
quency sufficiently accurate for all practical purposes. In the present instance, 
the simplest algebraic expressions consistent with the boundary conditions have 
been chosen to represent f(x, y). Assuming a sinusoidal vibration, this leads to 


w= Kx y(a—x)(b—y)* sin pis ea 0 ee ee (225) 
for clamped edges, and to 


w= K,xy(a—x)(b—y)(a? + ax — “2B? + by — _ 42) idee Brak (2.6) 
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for supported edges, where K, and K, are constants, the suffix indicating the 
boundary conditions. It should be noted that, although (2.5) and (2.6) satisfy 
the boundary conditions, they are not solutions of (1.3). 

The potential energy can now be found for clamped edges by inserting w, 
as given by (2:5), in (1-6) : 


V = Clip 


es 1575 | Dit + Daa + Dot | sin pt) tas. (2.7) 


7 
and for supporting edges by inserting (2.6) in(1.6) : 


62a°b®K? , ral ce ed Wes 
Vi= eae | Ds + Dyat+ 734 P34 D2 eine Pree aes (2.8) 


The kinetic energy is given similarly for clamped edges by (2.5) and (1.7) as 


_ phaeb?K.2p? cos? pt 


aoa as ce (2.9) 


and for supported edges by (2.6) and (1.7) as 


ie 961 pha°b?K 2p? cos? pt 
ae a OS) ae Ce en ee UT ie vieley eiiayc 


In each case the potential energy can be written V=MK* sin? pt and the kinetic 
energy T=NK?p? cos? pt. The maximum V occurs when sin?p¢=1 and the 
maximum T when cos?pt=1. Hence from condition (2.4), p?(=47?v?)= M/N. 
_ Thus from equations (2.7) and (2.9) 


Ly 426 4 3 
Ve= alg | Dib + Daa + 7 Dab: || = | Saban anil i 
and from (2.8) and (2.10) 
fs(e756 ee COT oe |i 212 
v= ape ok | Dib! + Daa + 734 Ps4 b eT Seer (CANA) 


Comparison of equation (2.12) with the accurate solution (2.2), putting 
m-=—n=1, shows that the form of the equations is identical, and further, the 
numerical coefficients are very close in the two cases. ‘The coefficient outside 
the brackets in the accurate equation (2.2) is 7/2 = 1-5708, while the corresponding 
coefficient in (2.12) is m14/(756/31)=1-5719 ; the coefficient of the term in 
D,a2b? is 2 in equation (2.2) and 867/434 (= 1-998) in equation (2.12). 

In order to obtain an estimate of the degree of approximation involved in the 
clamped-edge solution (2.11), the Ritz modification of the Rayleigh method 
has been applied to the present problem. The procedure, with an oe 
on the one-dimensional case of the vibrating string, is given in Timoshenko s book 
(1937, § 62). The application of the method to the two-dimensional case of the 
vibrating anisotropic plate does not present any special difficulties, although it 
leads to somewhat laborious algebra and arithmetic. For the above reasons, 
only the final results, and such description of the process as ts required for their 
presentation, will be given here. 
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The forms assumed for w are generalizations of (2.5) and (2.6) : 
(1) Clampededges, w=X?*Y?(K,+K,XY)sin pt ; 
(2) Supported edges, w=XY(a? + ax —x*)(b? + by —y?)(Ky + K,XY)sin pt, 
where X =x(a—«x), Y=y(b—y). 

If K,=0, the above equations reduce respectively to (2.5) and (2.6). An 
essential feature of the Ritz method is that the values of K, and K, are adjusted 
at one stage in the process so as to make the frequency a minimum. With two 
constants, this leads to a quadratic in the frequency, one root of which gives a 
more accurate estimate of the fundamental frequency than that given by the 
Rayleigh method. 

For the purpose of comparing the various estimates, it is convenient to confine 
attention to square plates, and, in the case of wood and plywood, to particular 
values of the elastic constants. For square plates, equations (2.2), (2.11) and 
(2.12) are of the form 


V=Ci/07V pl. = = ere (2.13) 


Values of the constant c, have been calculated by the Ritz method for plates of 
birch 5-ply and birch veneer having the following “‘ flexural rigidities” D (C.G.s. 
units ; values based on elastic constants given in table 2 below). 


dD, D, Ds 
Birch 5-ply 0-191 x 101° 0-071 x 101° 0-044 x 101° 
Birch veneer 2:97 x 108 0-21 x 108 0-69 x 108 


For an isotropic solid D, = D, = Ds, and it is therefore possible, by introducing 
‘this condition, to reduce any of the equations so far given to the isotropic case. 
The exact value of c, for a square isotropic plate with supported edges is known 
(Timoshenko, 1937, equation 214) to be 7/D,, where D, = EA?/12(1 —o?), and 


Table 1. Calculated values of c, in equation (2.13) 
(Square plates) 


Clamped edges Supported edges 


Method Birch : Birch 
Isotropic Isotropic —|—————___, 
5-ply Veneer 5-ply Veneer 


Rayleigh} 5-730 VD, | 1-914 x 105 | 0-6755 x 10° | 3-1429 VD, | 0-9300 x 105 | 0-3357 x 105 
Ritz 5°729 VD, | 1-912 x 10% | 0-6750 x 10° | 3-1420 VD, | 0-9294 x 105 | 0:3355 x 105 


Exact 5°727/ Dy, — — 3°1416 VD, | 0-9293 x 105 | 0-3354 x 105 
a a  ( 


Tomotika (1936) has given an estimate of c, for a square plate with clamped 
edges which, though not strictly exact, gives a highly acturate lower limit. The 
various estimates for c, are collected in table 1. 

Even the largest difference in the above table is only just greater than Os 
and it can be concluded that the fundamental frequencies estimated by the 
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Rayleigh method are mathematically sufficiently accurate for the purpose of 
comparing theoretical with observed frequencies. The Rayleigh equations 
(2.11) and (2.12) have, therefore, been used in all subsequent calculations. 


§3. EXPERIMENTAL 


The general features of the method used to determine the frequencies can be 
seen from figures 1, 2 and 3. The wood or plywood plate was held in a suitable 
frame and two small iron plates A and C (figure 1) were attached by means of 
soft wax, one (A) at the centre of the test plate, and the other (C) a convenient 
distance away. A small iron-cored solenoid (B), fed by alternating current 
from a beat-frequency oscillator, was placed near to plate A, while a telephone 
receiver, stripped of its diaphragm and case, was placed so that the magnets 
were close to plate C. Any currents. induced in the coils of the telephone were 
fed to a cathode-ray oscillograph, and the experiment consisted in varying the 
frequency of the driving current supplied to B until the response, as indicated 
by the oscillograph, was a maximum. ‘The telephone receiver D was connected 
to the vertical plates of the oscillograph, and the driving current from the oscillator 
connected to the horizontal plates, in parallel with the coil B. In this way, a 


saea| 
Ke 


FROM 
OSCILLATOR 


OSCILLOGRAPH 


Figure 1. Experimental arrangement. 


Lissajous figure was obtained on the oscillograph screen, and it was found that, 
with the present experimental arrangement, the plates tended to respond to 
integral submultiples of the resonant frequency, in agreement with results, so 
far unpublished, on wooden strips. ‘The most powerful response in all cases 
occurred when the frequency of the driving current was one-half the true resonant 
frequency, as is shown by the fact that the Lissajous figure indicated two vibrations 
of the plate for each cycle of the forcing current. ‘The results given in table 2 
are all based on this response mode. 

The apparatus can, of course, be used to measure overtones, and a number of 
observations were actually made of the frequencies at which these occurred. 
Owing to certain difficulties, mainly theoretical, a satisfactory interpretation of 
these measurements is not available at present, and they will not be discussed 
further. 

The frames for holding the plates were made of wood and are shown in 
figures 1, 2 and 3. The frame for clamped edges (general view, figure 2; cross- 
section, figure 1) consisted of two halves bolted together, with the plate sand- 
wiched between them. ‘The frame for supported edges (figure 3) consisted of 
four separate grooved members held together by brass bolts and wing nuts, with 
the edges of the plate at the base of the grooves, as shown in the section diagram, 
figure 3 (8). 

6-2 
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For convenience of lining up and adjustment, the apparatus was assembled 
on an optical-bench base. ‘The frames were held at opposite corners by ordinary 
laboratory clamps attached to vertical rods ; the solenoid B and telephone receiver 
D were also mounted on rods, all of which were held in optical-bench stands. 

The plates were originally cut to 45-8 cm. squares, and were progressively 
reduced in size in order to provide measurements over a range of lateral dimen- 
sions. It can be seen from equations (2.11) and (2.12) that if D, and D, differ, 


Figure 2. Clamped edge frame. 


the frequency of an oblong plate of given dimensions depends on whether the 
direction of greater flexural rigidity coincides with the longer or shorter side of the 
plate, and provision was made for the testing of plates fulfilling both of these 
conditions. For example, two 45-8 x 45-8 cm. plates of obeche veneer were 
originally available (table 2). The first reduction in size was accomplished on 
plate 1 by shortening to 30-6 cm. the side a, lying parallel to the greater Young’s 
modulus, and on plate 2 by shortening to the same length the side 4, lying parallel 


Figure 3. Supported edge frame. 


(a) General view ; (6) cross-section. 


to the smaller Young’s modulus. On the second reduction in size, side 6 on 
plate 1 and side a on plate 2 were shortened to 30-6 cm., thus giving square 
plates again. This procedure was repeated once more, finally leading to the 
smallest size of plate, 20-4 x 20-4 cm. The dimensions given refer to the sup- 
ported-edge condition ; the clamped dimensions are smaller owing to the overlap 
of the frames (see figures 1 and 2). 

Young’s moduli in bending (£), and rigidity modulus (G), on which, by 
equations (1.2), the “‘ flexural rigidities” depend, were determined respectively 


85 


Frequency of vibration of rectangular wood and plywood plates 


EE DOE SNES SAE REE a Ee ee ee ee ee Le eee 


8L7é 96 VLE O06£€ OCr 8CC T6L 981 CBT C8L 817 O06T v8T “UD 4-06 ==? 
891 LOL = 00C 187 €Cl >= 90T (Go) 60T es +6 => +:07=9 ‘9-0 =” 
= a OCE =. = my 991 3 =a c= 181 = 98T 9:0€=9 ‘b-0T=” 
Lop OCT OST 061 L8T 80T +6 v8 C6 88 SOT HI T8 “Wd 9:0 =G=” 
06 jhe oF 8 OST cs <= es OS 89 a 09 = 9-0€ =9 ‘8-Sp=P 
— a TST = = = S6 = aa ore +6 = 76 8-Sh=4 ‘0-0€ =v 
= cae 19 L9 Ble ca ev Iv SE ve Or 6£ SE US 3°S7—=—g—v 
sasps paytoddngs *(-0a8/saTo49) sorouanbasy [eJUsUIepUNF paAtasqg (¢) 
S19 $S9 OS8 Of8 086 tvs HIS OES vCS OTS 9SS COS try WW 0:41 —¢=2 
C8E CLE ae 98P OS9 SVC ra Ive VCC 6CC 77 107 ee 0-41 =@ ‘7-47=P? 
a = 90L = rae ae Orr ce a ra 9647 co 90S C:¥C=9 ‘0-4-1 =? 
8L7é OLT 9CE OE COE OG SGG 00¢ COC 9ST COT Tol 681 ON CoC aa 
COT SST — Lol HSC 96 re S6 C6 96 aes ; 88 ra C-¥7=49 ‘€-6€ =P 
== a vSC me ae a C8I = oat os col ra GLaY €-6€=4 ‘T-47=” 
= = Srl Orl O9T aa C8 68 8L OL OL C8 L9 TW §-6§—G=—0 
sasps podurya *(-998/saTd£9) satouenberz [ejUSUIepUN; paartasqg (7) 
€S-0 | vLL:0 | TIZ-0 | 9IL-0 | 069-0 | 614-0 | 659-0 | 029-0 | £6€-0 | 66€-0 | O0F-0 | 06-0 | OFE-0 99/"ums) J 
Clr-0 6Lt-0 SCS@ cES-0 909-0 so) 60£-0 60£-0 STE-0 OTE-0 SLE-0 CG6E20 1672-0 ("wD) y 
06:0 +S-0 £60 88-0 06-0 61-1 98-0 8-0 Ir-0 St-0 85-0 bE-0 05:0) neers [ ed 
Sc:0 S¢c-0 Ler) LE-0 LE-0 Le) €£-0 €€-0 CGO) 61-0 61-0 L1-0 L1:0 soudp } BOLT 
80:-T 08-0 Gal ol 1 SVE 8-0 $S-0 $s-0 vI-0 61-0 O£-0 S0-0 LY-0 | or-OT X | “7 
v8 €-6 TAG 0-Cl £-6 LU. Coe SONS]! LS 624 Soe VL 6:9 UGE 
S}UBISUOD ITyseTy (T) 
Ajd-¢ | <Ajd-¢ (Z) (1) Ajd-¢ yog (Z) (1) (€) (Z) (T) (Z) (1) 
If adedQg, yorg i 
sersnoqd Ajd-¢ yosirg yooog uoogersy sypeqqQ, 
spoomAyg : sI90u9 A, 


86 R.F. S. Hearmon 


by flexural and torsional vibration experiments (Hearmon and Barkas, 1941) 
on approximately 2:5-cm. wide strips of the material cut parallel to and near 
the sides of the largest plates. ‘The Poisson’s ratios (c) were not measured, but 
were obtained from published values of o for the appropriate solid wood. In 
some cases, no relevant measurements could be found, and it was then assumed 
contraction along the grain 
extension across the grain 
across the grain, had the value 0-025. As in much work on wood and plywood, 
the contribution of the terms involving the Poisson’s ratios is small, and the 
exact values adopted are not very important. 


that the Poisson’s ratio for solid wood, for tension 


§4. RESULTS AND DISCUSSION 


The complete experimental results are given in table 2 ; the side a is taken 
as parallel with the direction of F, and 6 with E,, where £, >). 

In order to compare the observed and theoretical frequencies, it is convenient 
to express the equations for the latter in terms of elastic constants, rather than 
flexural rigidities, by substituting equations (1.2) in (2.11) and (2.12). Thus, 


1-044 - a a a 2 $ 
os @,/p Ey + (5) Ey + @ (0°57Byo0, +1°13G) | Joa) 


and 


0-456h a\* a\2 : 
alo Es (5) Ey + (5) (2-000, + 3-966.) | ee: (4.2) 


where A has been put equal to 0-99, and the numerical coefficients are given to 
three significant figures. : 

The observed frequencies are plotted against those calculated from (4.1) and 
(4.2) in figures 4and 5. ‘The lines in these figures are not necessarily the ‘“‘ best”’ 
fit, but are drawn through the origin with slope determined by the average 
value of the ratio (observed frequency/calculated frequency) for all the obser- 
vations. Numerically, this average is 0-772 for clamped edges and 1-238 for 
supported edges. ‘Thus, the calculated frequencies with edges clamped are 
consistently higher than the observed frequencies, while the reverse is true for 
the supported edges. This suggests that the actual conditions at the edges of 
the plates are always intermediate between the conditions appropriate to strictly 
clamped and supported edges (equations (1.4) and (1.5)). The difficulty of 
realizing experimentally the exact boundary conditions assumed in the theoretical 
treatment of this type of problem is well known, and has been encountered, for 
‘example, in buckling tests of plywood panels in shear (Norris, Voss and Palma, 
£945). 

Ideally, of course, the points in figures 4 and 5 should lie on a line passing 
through the origin, and having a slope of 45° ; the general tendency of the points 
to depart from this line is a measure of the effect just discussed. In addition, 
there is also a certain amount of random scatter, which is partly accounted for by 
casual experimental errors in the frequency measurements, or by variation in 


the material between the test plates and the strips used in the determination of 
elastic constants. 
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Superimposed on all these effects, however, there is definite indication that 
with clamped edges (figure 4) the best line through the points is curved, and not 
straight ; the curvature is such that the ratio (v observed/v calculated) becomes 
smaller at higher frequencies. The physical interpretation of this result is that 
the frame used for clamped edges becomes less- effective with thicker plates and 
smaller lateral dimensions, and under these conditions departures from the ideal 
clamped-edge conditions tend to be higher. 

In order to obtain further information on this effect, a mild steel frame was 
made to hold the smallest square plates in the clamped-edge condition. This 
frame was of the same dimensions as the corresponding wooden frame ;_ the 
two halves, however, were held together by twelve 4-in. bolts instead of the 4-in. 
size used with the wooden frame. The observed frequencies in the steel and 
wooden frames and the theoretical frequencies (equation (4.1)) are given in table 3. 
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Figure 4. All results, clamped edges. 


In each case the frequency is higher in the steel frame than in the wooden 
frame, and is, therefore, closer to the calculated frequency. ‘The average ratio 
(observed frequency/calculated frequency) is 0:67 for the wooden frame and 
0-83 for the mild-steel frame. ‘These results confirm that a considerable part 
at least of the discrepancy between observed and calculated frequencies can be 
attributed to the conditions at the edges of the plates. 

With the exception of the opepe 3-ply, the observed frequencies in the 
steel frame are lower than those calculated. This suggests that, even in the steel 
frame, the restraint at the edges is not sufficient to realize fully the ideal clamped- 
edge conditions. 

Inspection of figure 5 (simply supported edges) does not reveal any marked 
tendency to curvature, though statistical analysis (see table 4) shows that 
there is in fact a slight tendency to curvature such that the ratio (observed 
frequency/calculated frequency) increases somewhat as the frequency rises. 
It is rather surprising that the supported-edge frame should prove more 


(pazepnoyes) 
086 OS6 OLCI O9CT O6rT O€Z 062 O8Z O18 OLL 008 OLL Or Aouonbes.y 
S (euue.1y [293s) 
S OZ 096 OLTT OLTT OCCT LY9 8c9 LE9 £19 £6S cv9 $8S FIs Aouonbes.y 
S a ees 
= (surely Uspoom) 
n 819 *S9 | 0s8 0£8 086 VES Vs) SOE SE ycc. | 0] sun) noccee We cacm tare Aouenberg 
ca pad, a SRE ree Airy See a 
me a hd-g | © (T) Ayd-¢ aes (Z) (1) (€) (Z) (T) (2) (1) ea: 
se[snog odode) Ajd-¢ yorrg Pra yov0q uwooqey eyo0qQ 
es eee 


88 


W9 )--[=G=v + podurjo saspy ‘998/saqo40 ut sarouanbary *¢ Ie I, 


Frequency of vibration of rectangular wood and plywood plates 89 


satisfactory in this respect than the clamped-edge frame, because the former has 
quite apparent defects. There is a small area near each corner which is un- 
restrained, but, perhaps more important, the degree to which the wing nuts 
holding the frame together (figure 3) are tightened is a matter for personal judgment 
in obtaining a sufficiently firm fixing of the test plate without applying appre- 
ciable compressions in its plane. Nevertheless, the experimental results indicate 
that, although the conditions at the edges do not constitute true simple support, 
they are to a first approximation independent of the dimensions of the plate. 

In view of the rough proportionality between elastic constants and density 
of timber, it is possible to derive certain empirical relations which enable the 
frequencies to be estimated even if the elastic constants and densities are not 
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Figure 5. Ali results, supported edges. 


known. If it is assumed that elastic constants are exactly proportional to density, 
then equations (4.1) and (4.2) take the form 


Vie Cee Biccne eee (4:3) 
where c, is a numerical factor which depends only on the edge conditions and on 


the ratio a/b. - 

Two related consequences follow from (4.3) : 

(1) For given edge conditions, and given ratios a/b and h/a®, the frequency 
should be constant irrespective of the values of elastic constants and density. 
Examination of the experimental results for the veneers in table 2, all of 
which have a thickness of about 0-3 cm., shows that this prediction is approxi- 


mately verified. 
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(2) For given edge conditions, and a given value of a/b, a plot of v against 
hja? should give a straight line passing through the origin. The theoretical 
frequencies for square plates are plotted on this basis in figures 6 and 7, and it 
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Figure 7. Square plates, supported edges. 


will be seen that, apart from three exceptionally low frequencies in each figure, 
the points lie very close to the lines 


¥,=47 x 104h/a® (clamped edges), aianorae (Hea) : 
v,=23 x 10*h/a (supported edges). RAB COS fs), 
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The exceptionally low frequencies all refer to the opepe 3-ply, and serve as a 
reminder that this method of estimating frequencies must necessarily give results 
somewhat in error in those cases where the assumed proportionality between 
elastic constants and density does not hold. 

The theoretical frequencies for oblong plates have been analysed in the same 
way and have given the following equations : 


Clamped edges : 


e Glo 2/3, vp~45 x 10*h/a? ; 
QO = 3/2, v-~95 x 10*h/a?. 

Supported edges : 
a/b = 2/3, vz LUX LOT ae ; 
Bip == 312s v,~29 x 104h/a?. 


It will be remembered that a is the length of side parallel to the direction of 
greatest Young’s modulus. ‘The scatter of the points for a/b =2/3 is greater than 
that shown for square plates in figures 6 and 7, and the uncertainty in the predicted 
frequency is therefore greater in this case. 

The experimental results for square plates (table 2) have been analysed 
statistically on the above basis, both as a whole and in groups corresponding with 
the different sizes of plate. ‘The regression equations for predicting the experi- 
mentally observed frequencies in terms of y(=10*h/a?) are summarized in 
table 4. 


Table 4. Summary of regression equations 
(v=observed frequency, y = 10*h/a?) 


Clamped edges Supported edges 
a=b=39-3 cm. | v=42-9y— 7:3 | a=b=45°8 cm. | v=—23-6y+ 3:8 
a=—b—24-2 cm. |) v—34-4y--14-9. | a@=b=—30-6 em. p= Oey = Oeil 
a=b—14-0 em. | v=27-7y--/0°2. | a=b=20-4 cm. | v=31'8y—44-0 

All results VO ONTO All results v=28-3y— 2-9 


For the largest sizes, the equations are not very different from (4.4) and 
(4.5), which are derived from the theoretical frequencies. With clamped edges 
there is a noticeable decrease in the slope of the lines and an increase in the 
intercept on the y axis with decreasing size of plate. ‘These effects are reflections 
of the tendency to curvature already discussed in connection with figure 4. The 
supported-edge results show a slight tendency to curvature in the opposite 
direction when analysed in this way, although, as already mentioned, this ten- 
dency is not obvious on inspection of figure 5. In other words, the edge conditions 
in the frames tend to approach each other at the higher frequencies, and there 
are indications that at the highest frequencies the “clamped ”’-edge frame may 
even impose less restraint at the edges than the “supported ”-edge frame. 
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THE OPTICAL SINE-CONDITION 
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ABSTRACT. ‘The Sine-Condition is formulated on the basis of optical path differences, 
leading to an expression which is valid in the presence of spherical aberration. It is 
shown to contain results obtained previously by means of a “‘ ray-intersection ”’ analysis. 
The “next terms” are estimated and serve as an approximate measure of the error 
involved ; or, alternatively, as an indication of the size of field over which the sine- 
condition is a valid measure of off-axis aberration. 


system ensures freedom from the aberration of coma for objects that are 

but a small angular distance from the axis.. Because of this fact, the sine- 
condition has been used extensively by lens designers as a check on the off-axis 
performance of both telescope and microscope objectives. There are, however, 
respects in which this usage has been unsatisfactory. 

In the first place the condition is a valid measure of off-axis aberration only 
when the angular distance off-axis of the object point is so small that its square — 
can be neglected. There is an uncertainty in the size of field over which this 
approximation is permissible, given a required order. of accuracy in the result. 
Secondly, the physical meaning which can be attached to the violation of the: 
condition is—in my experience—generally not appreciated. The permitted 
departures from exact fulfilment of the sine-condition are consequently 
either wholly arbitrary or based on experience of “ what works”. It is true that 
Conrady (Applied Optics and Optical Design, vol. i, p. 393) gave a “ Rayleigh 
tolerance ” for coma, but this was simply quoted without proof. 

In what follows, an account of the imagery of off-axis object points is given 
which leads to a simple interpretation of departures from the sine-condition 


IE is well known that fulfilment of the so-called Sine-Condition by an optical 
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and also provides an approximate rule giving the size of field over which the 
condition is a valid measure of off-axis aberration. 

In figure 1, O, and Q, are object points (real or virtual), the one situated on, 
and the other a small angular distance off, the axis, O,Oq’, of an optical system 
having A refracting surfaces. The diagram shows a meridian section containing 
the off-axis object point Q,. ST is a section of a spherical wave-front with 
centre O,; : ‘ST Q is a section of a spherical wave-front with centre Q,. Thus 
O,Po, Q:RoPy and O,E,, Q,E, are rays along which elements of the waves proceed 
from the object points O,, Q,. Ry is on the sphere SgT'g. 

Ox’, Qy’ are the Gaussian images of O,, Q, respectively; S’T’, So’ T9’ 
are spheres centred on O,’, Qy’. The forms of the wave-fronts emerging from 
the system and proceeding towards O,’, Q,’ can be specified by their departures 
at a given moment from the appropriate one of the spheres centred on O,’, Qy’. 
Such a procedure has been used before.* 


\ST SURFACE. X™ SURFACE 


Figure 1. 


When the departure from the appropriate reference sphere is zero (or constant) 
for the whole aperture, the corresponding image point is free from aberration. 
The departure at.a given point on the wave-front will be referred to as the wave- 
front aberration. It is a function of the position of the wave element on the 
emergent wave-front, and will be measured, not as a geometrical, but as an optical, 
path-length. 

The incident wave-fronts and the reference-spheres of the emergent wave- 
fronts are conveniently chosen so that they intersect in Ej, E,’ the axial points of 
the entrance and exit pupils. In the diagram (simply for the purposes of clarity) 
these latter have been taken to be external to the system. ‘Thus Q,E, is the 
initially incident principal ray from Q,. 

An element of the wave from O, proceeds along the ray O,;P)P,Po...... Py Py11, 
intersecting the O,’-reference sphere in P,,,. Suppose the final wave-front 
aberration of this element is W,’, measured positively when the element is in 
advance of the O,’-reference sphere, and such that this latter and the emergent 
wave coincide at E,’—that is, W,’ is zero for the element of wave which proceeds 


* For a derivation of the theory of the aberrations of lens-systems on this basis, see the writer’s 
“ Monochromatic Lens Aberration Theory ”’, Phil. Mag. (publication impending.) 
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from O, along the axis 0,0)’. Then the wave-front aberration of the wave- — 
element considered is 


Wy’ =(E,Ayhy. 22. Ag En) (Poi Pee ee (1) 


where the brackets denote optical path lengths. 

A second element of wave proceeds from Q, along the ray QyRoP)R,Ro-....- | 
R,Ry4, intersecting the Q,’-reference sphere in Ray,. If this element has a | 
wave-front aberration V,’, measured positively when the element is in advance | 
of the Qy’-reference sphere, and such that this latter and the emergent wave 
coincide at B,4,—that is, Vy’ is zero for the element of wave which proceeds 


from Q, .along the principal ray E,B,B,..... B,B,, ;—then the wave-front | 
aberration of the element considered is 
Vo = (i Bi Bacar By By) —ERo Pp Rae tee ee (2) 


where B,+;, Ra, 1 are the points of intersection of the emergent rays with the ! 
Q,’-reference sphere. Rewriting (2), 


Vy’ =(E,B,B,..... ByBy+1) —(PoRiRe.-.-- RaRa+iR) + te aR —RoPo), 


where R is the point of intersection of R,R,, with the O,’-reference sphere, it 
follows, using (1), that : 


V2 Wx (Rig RORS Ee) Ce (3) 
where 
CS Bahar ByByy1) —(ExAyAg----- A,E,’)} 
—=(( PRG Ra dees Ri Rie) aCPsP Pee PyPya)¥. 


It can be shown that C, is zero (or negligible) for small angles of field. 

From P, draw P,h,, P,h,’~perpendicular to P)P,, P,P., and meeting P)R,, 
R,R, (produced if necessary) in h,, hy’. Then the following are differences in 
optical path length : 

(PpR, — PpP,) =N,. 4, sin], +N. Dp (sec x, — 1) 
(R,R,— P,P.) =N,’ . J,’ sin I,’ —N,’ . 1, sin I,’ + N, . Dy(sec x, — 1) 


(RR —P, Pasz)= Ny’. Asin y+ Nya, - Dy (seexa4,—1), 


where/, = P,R,,/,=P,Rg, ..... ete:5 %,=R, P,P. fh ete.; D> = P,P, De Pee 
ets A Dy= PaPy413 Ny, Neg=Ny',......etc., are refractive indices; and J, r’ 
are the angles of incidence and emergence of the ray P)P,P,..... P, at the 
successive surfaces of the system. Addition of these equations gives, the angles 
x being assumed small, 


(MEARE. RaRag RR) (PoPiFoen oe PxPaon)) 
= $4" NyDo+ $%27. NoD, +222: $40h 4 Naa eles (4) 
the other items cancelling because of the law of refraction between J and 7’. 
By a similar argument applied to the principal ray it can be shown that 
{(E By Boreas ByBy+1)—(E1AiAg...... A,E,’)} 
=ty,?. Nid) +42: Nodi+...... “hye a. Nyt) cee =t5} 
where dg= Ha Ay dpe eee etc.; and4¥,=B,E,A), and so on. 
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From (4) and (5) 
C=4y2 Wald, Dy) +4922 Ned Di)... 49941. Nigal—Dy), 
where a “paraxial”? approximation has been made in putting %,=V2, %;=Ys5 


Se etc., and where x;=y,=w1; w, being the angle of field Q, E ,0,, which 
is small. If, now, retaining the paraxial approximation, there is written 


(62) (6) (62) 
Tae V;, ora Vo, 000. Oa! Vata 
LNT) Zo on ae V4; are angles from a traced paraxial principal ray, then 


2 
C,=4 (7) Na deal oN (de Deine 
1 
NM dD a ea in eee (7) 


which gives an approximate value for C,. For small angles of field, powers 
of w, higher than the first are neglected, and C, is regarded as zero. Equation (7) 
then provides an estimate of the most significant error arising from this approxi- 
mation. On the other hand, given a “safe” limit for C,, (7) provides an expres- 
sion for w,, the size of field over which the approximation is sufficiently valid. 
C, is on an optical path length, and is one of the principal errors involved in 
V,’, when C, is assumed zero in equation (3). 

Attention has been confined so far to rays from Q,, lying in the meridian 
section. It is, however, easy to show that the result (4), and therefore (7), hold 
equally for skew rays from Q,. For, in this case, the elements of refracting 
surtaces at/P;; P...<.- P, are regarded as plane. Planes are constructed through 
each of these points perpendicular to the rays incident and emergent at the 
points. The differences in optical path length then become 


(PoRy — Po Py) =N,./, sin L,. cos ¢,+ Ny. Dy {sec ¥, — 1) 
(R,R, — P,P.) = N,./, sin I, . cos ¢2—N,.1, sin I,’ . cos ¢, + Ny.D,(sec x, — 1) 


and so on, where 4j, ¢p, ..... etc., are the azimuths of the points of incidence 
| ea Ry. On addition, terms on the right-hand side again cancel, giving 
the result (4). Thus C, is zero for both meridian and skew rays, providing 
powers of w, higher than the first are neglected. Furthermore, equation (7) 
still provides an estimate of the error involved in the approximation. 

Construct, then, rectangular co-ordinates having O, as origin, O,O,’ as 
x-axis, and Q, in the XOZ-plane. Let O,E,=7, Q,E,;=p ; and let Py, (on the 
incident wave from O;) have co-ordinates (x, y, z). Then if Q, is the point (0, 0,%), 
the length Q,P, is given by 


Q,Ppe=at ty?+(s—7)?=p?— 22m. 
In the diagram 7, and w, are shown negative, following the usual optical 
convention of signs. 


Thus, since O,E,Q, =, is assumed small, 
I) eee BOI ia 
RyP) = Q,P) -—p = —-— —-3 . 
oo =Qilko—P ? Pp 
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3 


2 
or, avitisteuahies = _! (1 a "), 
Pp Y 


Replacing z by polar co-ordinates (p,, ¢), where ¢=0 is the positive direction 
of z, and p, is the perpendicular distance of P, from the axis, this becomes 


cos 
R,Py= — APY cos p— BEY cos. } tan® Ww, (sane, eee) Sorters (8) 
1 18 


where the suffix 1 has also been given to y, w and 7. 
Similarly if Qy’ be joined to P,,, (and continued), the intercept between the 
two reference spheres is very closely equal to Ry, ; and, therefore, . 
a , , / , iY cos 
Ry R=— ENS . cos ¢ — TAP cos . tan? wy,’ (| tan aw,’ — eee’) ) 
Tr Tr a Tr 


where p,’ is measured at the exit pupil. 
Subtraction of (8) and (9) gives for the optical path difference 


os N. N: Te ea: 
(RaziR —RoPo) = (“a - Meee) cosh-KOs yt ieee (10) 
nN 
where C,, a second-order term, is 


cos 
Cle Nii .cosd. stan? w,? (san w1— ane) 
% Ty 


N- / / , fe cos 
= AMPA | cosd. 1 tan®an’ Cit eg ate : 
Tr Tr 
It will be seen that C, is zero in the sagittal section, 6=7/2. Inthe meridian 


section it has the value 


N. / 
C,= BP 1 tan?w, (tan w,— 2% 


N- , y Le , 

A EN anton janen eae wy ae (11) 
Tr Tr 

~ This is zero if powers of w, higher than the first are ignored. (11) then estimates 

the error involved in the approximation. 


Substitution of (10) in (4) gives, ignoring C, and C,, 


Vy! =W,! + |e “ “ane isk (12) 
Fy Tr 
where 7, 7’ have signs such that, in figure 1, 7, is positive and 7’ negative. This 
will be so if 
Ue ES = beer } 


aN = lor, bee U yarax., 
with an obvious notation. Thus the wave-front aberration associated with the 


wave proceeding to Q,’ consists of the spherical aberration W,’ of the wave 
proceeding to Oy’, together with a coma-type term, which is zero in the sagittal 
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section (6 =7/2) and has a maximum value in the meridian section (p=0). The 
latter is merely reversed in sign in the section 6=7. In the section o=0; (12) 
becomes 

Wee |e aN Pa i oe (14) 

l pr, “parax, l Dan l parax) 
and this can readily be given forms suited to special cases. 
Thus, for telescopic systems in air (N, = Ny’ =1), 
Vi’ = W,' +[p,.tanw,—p)’.tane’], — ...... (15) 


where p,, pa,’ are measured in the (paraxially defined) entrance and exit pupils 
along the marginal ray. If h,, 4’ are the initial and final incidence heights of a 
paraxial ray trace, the magnification is 

h, _ tanwy’ 

Te tan Wy ; 
so that 


oe hi ? 
Vi = Wp +[ a Fay’ | tan ete (16) 
A 


Using the notation of figure 2, in which 8m’, the final angular aberration is 


MARGINAL RAY 


-o. 
y~ = 


\m SURFACE. 


EXIT PUPIL. 
Figure 2. 
shown as positive, 
Pei Ca Stan ony ee (17) 
From this it will be seen that, in the presence of spherical aberration (6u,0), 


pa’, and hence the coma part of Vy’, will vary with /’,,, the position of the exit 


pupil. 
The case of telescopic and photographic objectives in air gives 


ti / 
V>'=W,' + P= mv | tae re doo (18) 
l Pry P paraxy 


where use has been made of the relationship 7)’ = F' . tanw,, F being the equiva- 
lent focal length of the system. In figure 3, m’ is the angle between the finally 
emergent marginal ray and the axis, Sua’, ds,’ are the angular and longitudinal 
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spherical aberration respectively of the marginal ray PPOs 
angle 5’ being small, and negative as shown, 


Evidently, the 


py =7y' . Sin (ua’ — buy’) 
=r! {sin uy’ + Os’. sin uy’. cos uy’ /7y"}, wee ee (19) 
since, with the given sign convention, 
Se 05x’ . on Ur’ 
Tr 
Using this in (18) gives 
: { sy’. ‘ 
Var=W.+ E — Fsin Ua se aoe] tal wy) Sip h eae (20) 
l PIA l paraxy 
It will again be seen that. the coma part of V;’ is independent of shifts of the 
diaphragm in the absence of spherical aberration (6s;’=0) : on the other hand, 


Qu 


» 
MARGINAL RAY. 


’ 
Oo, Ss\ Corn 
PARAXKIAL MARGINAL 
RAY FOCUS Ray FOCUS 


EXIT 
PUPIL. 


Figure 3. 


V,’ is independent of the spherical aberration when /’,,, = co—that is, with a 
diaphragm in the anterior focal plane (telecentric stop). With the exit pupil at 
the posterior nodal point, F= — (Zor, — E'paraxy) = —7n, and 

Vy = Wy! + [py —(F'— 8s)’ . cos Up’). sin U)’]tanw,, —....... (21) 


or, since in this case w, =a’, 
f F—6s,’.cos U,' 
Lex = Wy’ a frees = = 
where w'parax, Tefers to a paraxial trace having the same initial incidence hei ght 
as the marginal ray. 
For systems in which both object and image distances are finite, (14) becomes 


sin uy | eee (22) 


Ve = W,’ ae | Non sin US Nn ay sin U,' F (1 ae ree), 
l pr paraxy 


(23) 
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use having again been made of the result (19). Again it will be seen that V7,’ 
is independent of 8s)’ when /,,’= co ; and independent of Pom when 6s,’ =0. 
These results were obtained ee Conrady (loc. cit. p. 375) using a “ray-inter- 
section” investigation. (23) can be put into a different form, thus 


ae E U, _ sin Pa 4 § Osx’. cos Uy on cos U,' ) ee ee (24) 
uy Un Lom E paraxy 
prnere Al Nia = IN yw ai = =N,'u'n (the Lagrange Invariant) ; 


and small w’s refer to paraxial angles. When 5s,'=0, or lpn, = 00, this gives 
for the absence of coma the classical (but very restricted) result 
sin U, _ sin Uj’, 


Us Uy’ 


oe (25) 


The above formulae will give the coma part of the wave-front aberration 
even in the presence of large amounts of spherical aberration. The Rayleigh 
tolerance states simply that the total V,’ must not exceed A/4 ; Conrady’s toler- 
ance applies only to the coma part of V,’. » Given the same tolerance on Wj’, 
N,’ might thus reach 4/2 at some parts of the aperture. 

The coma-part of the expression for V}’ includes all the coma terms depending 
on the first power of the field, i.e. those in which the wave-front aberration is pro- 
portional to the cube, fifth, seventh..... powers of the aperture. The accuracy 
of the results can be easily checked by noting the values of Cy and Cy. If, as ina 
good many telescope objectives, the incidence heights of the principal rays are 
all small, there can be written in (7) 


VeVi Nas 
and so on, giving the useful result 
C, =4N,,"[(dy — Do) + (4, — D) + (4, — Dg) +... + (d, —Dy)]. vee (26) 


‘The correction term C,, depending on powers of w, higher than the square, will 
generally be smaller than C,. 

Although the above demonstrations have assumed the (more usual) system 
consisting of centred spherical refracting interfaces, it will be appreciated that 
essentially the results apply to other systems. Moreover, it will be seen that a 
similar analysis yields a ‘“‘ Sine-Condition ”’ (based on an accurate principal ray- 
trace and another ray from the same object point) which will give the wave- 
front aberration associated with image points in the neighbourhood of an off-axis 


image point. 
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HERSCHEL’S CONDITION 


By H. H. HOPKINS: 
W. Watson and Sons, Barnet, Herts 


MS. received 12 October 1945 


ABSTRACT. The change in spherical aberration introduced by a change in object 
distance is discussed. Herschel’s condition is shown to be valid for originally corrected 
conjugates. It follows as a special case of the general formulae derived. 


HE satisfying of Herschel’s condition by an optical system which is free 
from spherical aberration for given conjugates ensures freedom from 
spherical aberration for neighbouring pairs of conjugates. Attention 

does not seem to have been called to the invalidity of the usual form of the con- 
dition in the case when the system is not corrected for the original conjugates. 
It is, however, possible to formulate the condition in such a manner that it not only 
applies to such uncorrected systems, but such that the resulting generalized 
condition also measures directly the change in spherical aberration consequent 
upon a small change in the conjugates. The treatment resembles closely that 
which has been used elsewhere in connection with the Sine Condition (Hopkins, 
Proc. Phys. Soc. 58, 92 (1946)). 


In the figure, O, 1s an object point (real or virtual) situated on the axis O,O,’ 
of an optical system consisting of A refracting surfaces. Oy,’ is the Gaussian 
image of O, after A refractions. Spheres ST, S’T’ centred on O,, O,’ serve as 
reference surfaces, with respect to which the initially incident and finally emergent 
wave-fronts can be specified (vide Hopkins, 1945). These are conveniently 
chosen to intersect the axial points of the (paraxially defined) entrance and exit 
pupils. The latter are, respectively, Es and E,’.* | : 

Assume the incident and emergent wave-fronts to coincide with the O,- 
and O,’-reference spheres, respectively, at the points E,, E;’. Let the element 
of wave from O, which is associated with the ray through Py be in advance of 
the O,-reference sphere by an amount W,, which is an optical path length. This 

* In the presence of spherical aberration the outer parts of the wave-front close up (or open out). 
more rapidly than the central parts. Because of this, the angular distribution of wave-front 


aberration depends upon the radii of the reference spheres. These are conventionally chosen 
in each case, to refer the aberration to the entrance and exit pupils. , 
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is, then, the initial wave-front aberration of the wave-element. Let W,’ be the 
amount this element is in advance of the O,’-reference sphere upon emergence. 
This is, then, the final wave-front aberration of the wave-element. Since optical 
path lengths along rays between corresponding wave-fronts are equal, 


Wha Via (PGP Pe eecre pha SR pAG As. AREA), 2 tees (1) 
where the brackets denote optical path lengths. 
Suppose, now, the object point be moved to Q,, whose Gaussian image is 
Q,’. Let a wave-element, having initial wave-front aberration W,, proceed 


from Q, along the ray RyPpR,R,.... RaRay, (where Ry, Ray, are on the spheres 
SeTa, Se’Te’, which are centred on Q,, Q,’. There can then be written 


Ve AVA RAPE RS Rowe ts, Rye) = (Ey Ay Aen AGE 2) eee (2) 


where V,’ is the final wave-front aberration of the wave-element on emergence 
from the system, measured with reference to S’g'T'g’, the Qy-reference sphere. 
On subtraction (1) and (2) give 
Ve Ver CP GP) Paces Da ate | aig Uigktke acess RRs) 


=P EPPN eS ae. Py Pe Cea Ra Ry R) 
+ (Ra+1R — RoPo). 


The first two brackets have been shown elsewhere to be equal, and they 
cancel, providing the angle O,R,Q,=AU, is sufficiently small for powers of 
AU, higher than the first to be ignored. ‘The error implied in this approxi- 
mation has also been discussed elsewhere (Hopkins, /oc. czt.). ‘To this accuracy, 
the change in the wave-front aberration of the element considered is, therefore, 


dW =(RyeR — Rabe = a2! Sea. (3) 


Construct, then, plane rectangular co-ordinates having QO, as origin, x-axis 
along the line O,Oy’, and y-axis in the plane of the diagram and positive above the 
optical axis. Let P) have co-ordinates (x, y), and let E, be the point (7,, 0). 
Then 

PoQy? =(x— Ax)? +9" 
SSP, INK, 


since P, is on the circle x*+y?=7,?. Ax,=O,Q, is assumed small compared 
with r,, and powers of Aw,/r, higher than the first are therefore ignored. ‘Taking 


square roots, 


aire, 
P)hQ= 14 - 


fl 
to the required accuracy ; and further, 


OR, = 0,5 =7, Ax. 


Thus 
P,Q, — QURo = Ax,(1 = =) 
Lal 
or 
PeRp ae (lcosUsy kame (4) 


where 6, is the angle P,O,E,, O, being the paraxial object point. 
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If, now, Oy’ be joined to P a4, (and continued), the intercept between 
the spheres S’T’, Sq’Tq’ is very closely equal to Ra;,R. —_ Using, therefore, 
an expression for Rawk similar to (4), the expression (3) becomes 

dW=Ny,, . Ax,’(1 —cos 84’) — N, . Ax,(1 —cos 9,). 
Thus, using the (paraxial) longitudinal magnification formula 
Ny Ay = NC Ad ee (5) 


the change in spherical aberration produced by the change of conjugates denoted 
by Hy, is given by 


iW eos mas ee). i. oo 6) 


where uw’, u, refer to paraxial angles associated with the original conjugates. 
This formula can be given forms suited to special cases. 

Thus, in the absence of spherical aberration, 0,=U,, )’ =U)’, where Uj, 
U,’ are angles from a marginal ray trace. (6) becomes 


_jifsmetx ¥ ne) | * 7 
sor (SOE (OEE 0 


which leads to the usual form of the Herschel condition. For dW will be zero, 
providing 


singU,’ — singU, 
Tp PAR pn) 2 
If the system is uncorrected for the original conjugates, 0,’ 4U,’, 0,4U, : 
(7), (8) are then no longer valid. 
If, now, 


ds, . sin U, ‘ 85, . sin U3’ 


du, = ai ) dun ite [’ 


i? 
lor, a lparax, PIX l paraxy 


are the initial and final angular spherical aberrations associated with a given ray, 
defined when the axial parts of the incident and emergent wave-fronts are at 
E,, Ej’ respectively, and with the notation of an earlier paper (loc. cit.), then 


0, = U,—8u,, 0,’ = U,’ — du,’ 
and 
ds, sin U, 


cos 6,=cos U,+sin U,. j I 5 
DIi 9 “parax, 


ds,’ . sin U,’ 
cos @,’=cos U,’ + sin Uy’ pee ets 
l PIA = paraxy, 


Substitution of these in (6) gives, after some reduction, 


sin $U,'\2 1—cos U,’ 
aw-| (“Se ).(1 1-cosUy 5, 
Un [2 A Lom —E parax, 38 


7 a) 1—cos U, 1 j 2 
( orp 1+ ree ye btrg Moreen 88) 


rere a 
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In telescopic systems, the value of Ax is given in angular measure. Thus 
it h,, h, are the paraxial incident heights at the first and last surfaces of any system 
and a change of conjugate Au, results in a change Au,’ in the object space, then 
generally 

INE IE ge ALSO oe 


Ba ait a 


and, therefore, the longitudinal magnification formula becomes 
Nyx hy Au’ = N,h,Au, = Hy, ay yh © Wcaoiteneen (1 1) 


It is this form which is used in telescopic systems. Let p,, pa’ be the inter- 
section heights of a given ray in the paraxially defined entrance and exit pupils ; 
then for small angles Uy’, U, (a “nearly” telescopic system), 6,’, 6, are necessarily 
small, and 


ji Pr. ; hy 
a ON Seas 
Lom — U paraxy, l paraxy 
In the limit E paraxy > <0, the ratio 
US 
ur’ hy 
Thus, for telescopic systems, 
en ’ \ 2 
He OSE Oe ny re (12) 
Un? 2 hy 


with a similar expression in the object space. Using (11), (12) in (6), the change 
of spherical aberration, consequent upon the change of conjugates Hy, is 


dW = {| (S) A (2) | cine Maer (13) 


It has been shown before (Joc. cit.) that 
pa = Y, = (on, — X;) . tan Ou; , 


with the notation of the earlier paper. In the absence of spherical aberration, 
Herschel’s Condition for a telescope is thus 


1—cos@ 
It will be apparent that use of one or other of the forms of ( a. ) used 


in (10) and (13) will give results suited to photographic objectives (infinite object, 
finite image distance) and eyepieces (finite object, infinite image distance). 

If it is assumed that for given conjugates a system is both free from spherical 
aberration and satisfies the sine-condition, there follow immediately results 
which have a bearing on a number of questions. 

Thus, consider a corrected microscope objective of N.A.=N, sin Uj, having 
a primary magnification m. A change of tube length Aw,’ gives, N,’ being 
that of air=1, 


2 
Hy, =u, . Ax,’ = (“) EAN, 


104 H. H, Hopkins 


The sine-condition is satisfied, so that 


sind Uy Uy ssi Oy 


pee tues Yee 2 
u,’/2 eo u,/2 Boel es 


or, the light coming from the short conjugate, u,’ is small, and 


sin $ x) (= 4 ) -(- —cos =) 
u,//2 tia uy , 


Substituting in (7) gives 
_ 2 
i joe PAR. ee (15) 


Thus for a oes objective having N.A.=1-20, m= x 100, and using an oil 
for which 1:53, 


dW = —0-0000167 . Ax,’, 


i.e. a change of wave-front aberration of $ wave-length of sodium D light follows 
a change in tube-length of 17-5 mm.* With an objective of this type “ good 
observers”’ claim to estimate tube-length to this order of accuracy, although the 
effects are generally unnoticeable at the “‘best”’ focus. 

In a corrected telescopic system the sine-condition takes the form 


which is identical with equation (14). ‘Thus, in this case, Herschel’s condition 
and the sine-condition are simultaneously satisfied. 

Consider a corrected photographic objective, satisfying the sine-condition. 
Let the long-conjugate change from oo to Ly, and let h, be the radius of the 
entrance pupil. ‘Then by (11), 


2 
ie 


H- —— 
L i 


assuming N,=1. ‘The sine-condition gives 


or, Ff being the equivalent focal length, 


Gl i (at) = 1—cos U,’ oe: 
| h, u,//2 h, 


4 Since writing this, Mr. A. G. Sabin has given me a similar value for a }-in. versalic objective 
obtained by accurate ray-tracing methods. For N.A.=1°16, m= x98, a change in aberration of 


1-6 wave-lengths was found to be produced by a change in tube-length of 50 mm. The above 
approximate formula gives 1:3 wave-lengths. 


+ This is also true at unit magnification for a copying system corrected for spherical aberration | 
at that magnification. 
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The change in spherical aberration is, therefore, 
Fe 
aW= 4(1~c0s ui) = ot OE ages (16) 
Thus for an #’/1-0 objective, 


fie 
dW =0-0056 .— 
aoe 


and, for a 1-in. focus lens, a change in aberration of one wave-length occurs for 
L,=242 inches. Further, for a given relative aperture, the value of L, for a 
given change in spherical aberration is proportional to the square of the equivalent 
focal length. For a given focal length it is approximately proportional to the 
fourth power of the relative aperture. 

It will be seen that an identical analysis yields a ‘‘ Herschel condition” applying 
to extra-axial object points. For first-order changes in the conjugates, the 
aberrational change is symmetrical about the principal ray (i.e. a “spherical 
aberration”’) and greater changes in the conjugates are required to introduce 
changes in the comatic and astigmatic aberrations. 


THE-THERMAL CONDUCTIVITY OF TEXTILES 
By S. BAXTER, 


Wool Industries Research Association, [Leeds 
MS. received 18 September 1945 


ABSTRACT. A review is given of measurements of the thermal conductivity of natural 
clothing fibres and of formulae that have been proposed for the overall. thermal con- 
ductivity of fibre-air mixtures in terms of their density and fibre conductivity. 

Measurements of the conductivity of wool felts ranging in density from 0:01 to 
1-0 gm./cc. have been made, using Hercus and Laby’s method for determining the con- 
ductivity of air. "The conductivity of solid keratin (horn) has also been measured by 
the Lees disc method, and it is suggested that wool fibres and horn have identical con- 
ductivities of 4:62x10-*c.g.s. units at 0°-7% regain and 5:32x10-*c.g.s. units at 
10:7% regain. The data given for wool felts and horn are used to obtain an empirical 
relation between conductivity and bulk density. 

Results obtained with cotton fabrics indicate that published values of the thermal 
conductivity of cotton fibres are too high. 


Sil INPRO IDNEKE AB IKON 


two components: the heat transmitted through the air spaces and the 

heat transmitted through the textile fibres. ‘Typical wool fabrics, as used 
for clothing purposes, have densities ranging from 0-1 to 0-3 gm. per cc., Le. 
the air space in the fabrics varies from 92°/, to 77% of the volume of the fabric. 
Clearly this means that in a fabric the thermal conductivity of the air will be the 
controlling factor and the conductivity of the fibres will exert only a secondary 
influence ; in particular the conductivity of a fabric of low bulk density will be 
practically independent of the fibres used. ‘This result has been shown by many 
workers investigating the thermal insulation of fabrics. 


T's thermal conductivity of fabrics or masses of fibres is made up of 
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If any differences in the thermal conductivity of various fibres are to be shown, 
the thermal conductivity measurements must be made on fabrics of high density, 
where the conductivity of the fibre becomes increasingly effective. 

Few attempts have been made to determine the thermal conductivity of textile 
fibres from measurements of the conductivity of fibre-air mixtures. It was 
well known that low-density fibre-air mixtures had thermal conductivities 
approaching the value for air, and Speakman and Chamberlain (1930) have shown 
that the conductivity increases with the bulk density, depending on the type of 
fibre used. 

The present work is an attempt to find an empirical formula for the wool 
conductivity-density curve from measurements of the conductivity of felts of 
increasing density and also of the conductivity of solid keratin. A satisfactory 
formula has been obtained, and the conductivity of any natural fibre can be 
predicted from its conductivity-density curve. Also, if the thermal conductivity 
of the fibre can be determined directly, as is the case with most synthetic 
fibres, then the conductivity of a fibre-air mixture of any density can be 
predicted. 


§2. HISTORICAL 


According to Schuhmeister (1877) the first investigations of the conducting 
power of cotton and wool fabrics were made by Peclet and reported by him in 
his Handbook of the Science of Heat. He concluded that the conducting power 
of various fabrics was little different from that of air. Later, Forbes (1872) 
found that the conducting power of cotton decreased with compression, a result 
which Schuhmeister interpreted as meaning that the conductivity of cotton 
was less than that of air. ‘The present results show that it is quite possible, under 
certain conditions, for the conductivity of a fibre-air mixture to decrease as the 
density increases, although the fibre conductivity is greater than that of air. 
The problem of the conductivity of textile fabrics was first attacked seriously by 
Schuhmeister, using the method developed by Stefan for investigating the 
thermal conductivity of gases. In this method, the fibres are packed into the 
space between two concentric cylinders, the centre cylinder acting as an air 
thermometer. Schuhmeister showed that his results could be expressed by the 
formula 


k= hi tab) (1) 


where &,,, is the conducting power of the mixture, k, the conducting power of 
the gas, P the weight of textile fibre present, and a is a constant. 


Taking the volume occupied by the fibres into account, the formula reduces 
to 


ie = ky af bp, 


where p is the bulk density of the fibres and 4 is a constant. For wool the value 
of the constant is 2-71, for silk 2-50 and for cotton 7:35. The range of values of 
bulk density used was small, the highest value being 0-15 gm./cc. In calculating 
the thermal conductivity of the fibres, Schuhmeister assumed that the fibres 
were uniformly distributed in all directions, or, for the purposes of calculation, 
that one-third were plane slabs parallel to the direction of heat flow and touching 
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both plates, and two-thirds were slabs at right angles to the direction of heat 
flow. On these assumptions, 


kik 
y= Main the)+§(Ge—), a, (3) 


where k,, k, and 7, vp are the conductivity and fractional volume of air and of the 
fibres, so that 


Using equation (3), and assuming the conductivity of air to be unity, Schuh- 
meister gave the value for cotton as 37, for wool 12 and for silk 11. 

Rubner (1895), using the Stefan calorimeter and Schuhmeister’s method of 
calculation, gave the conductivity of wool as 6-1, of silk 19-2 and of cotton and 
linen 29-9, taking the value for air as unity. Rubner measured the thermal 
conductivities of a wide range of keratin fibres and found no significant difference 
between them. ‘The same result was also found for a wide range of cellulosic 
fibres, and for natural silks. His absolute values of conductivity are given in 
table 1, but it should be noted that they are not consistent with his ratios of the 
conductivities. 

Table 1. Rubner’s values 


Thermal conductivity 


SHS ete k (c.g.s. units) 
Air 0532) —10—- 
Mammal hairs 4-791 
Silk 8:870 
Plant fibres 14-199 


Rubner also investigated the increase in the value of k,, with increasing 
density, but again he investigated only the region of very low densities, his 
maximum value being 0-12 gm./cc. His values for the constant in equation (2) 
were 3-65 for wool, 3-62 for silk and 5-9 for cotton. 

One of the main reasons for the discrepancy between the results of Schuh- 
meister and Rubner was their choice of working in the region of low bulk density, 
where, as will be shown later, the conductivity of fibre-air mixtures behaves in 
a peculiar manner. One might have expected, however, that the agreement 
would have been better since Rubner’s work appears to be an exact replica of 
Schuhmeister’s. The error in estimating k, from equation (3) must be very 
large when low bulk densities are used. The higher the bulk density used to 
determine ,,, the smaller will be the error in estimating ky. The validity of 
equation (3) is discussed later. 

The theory of the conductivity of heterogeneous media is complex and has 
only been attempted for certain simple cases. Maxwell in his Electricity and 
Magnetism has worked out the conductivity of a medium formed by embedding 
spheres of conductivity k, in a medium of conductivity k,, the spheres 
being a large distance apart. He also worked out the cases where the media 
were separated by planes both parallel and at right angles to the direction 
of heat flow, and the case of right prisms of rectangular cross-section, one dimension 
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being large compared with the other, and the length of the prisms being large 
compared with the larger cross-sectional dimension, placed either with their 
edges parallel or perpendicular to the direction of heat flow. These cases 
have limited application on account of the geometry implied. Lees (1900) 
has also worked out an expression for the conductivity of heterogeneous media 
from considerations of the heat flow through a medium formed by an equal 
number of equally long prisms of square cross-section having conductivities 
k, and k, and arranged in draught-board pattern, the planes of constant tem- 
perature being drawn through the diagonals of the cross-sections of the prisms. 
For this case Lees shows that the effective conductivity k,,, is given by 


Ry hyhsts ee eee (4) 


From this result he reasons that the thermal conductivity of a compound medium 
in the general case where the constituents are in the form of long prisms with their 
axes perpendicular to the direction of heat flow is given by 


Fey = Ryn + oe (5) 


where v, and vy are the fractional volumes of the media of thermal conductivity 
k, and k, respectively. It must be noted, however, that, even in the general 
case, equation (5) only holds when the prisms are square in cross-section and 
are arranged in some form of draught-board pattern with their edges touching 
along their length. 


§3. EXPERIMENTS WIEFH THE LEES DISC 


The experimental method first used by the author to determine the absolute 
conductivity of felts and fabrics was that of the Lees disc (1898). The apparatus 
is sufficiently well known only to require brief description here. Its essentials, 
shown in figure 1, are the brass plates C and U, which are screwed together with a 


q— 
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Figure 1. Lees’ disc apparatus. 


heating coil between them. The fabric to be tested is placed between this hot 
plate and a lower plate M, and the whole block is then suspended in an en- 
closure of constant temperature until the readings of heater current and voltage 
and the plate temperatures, which are measured by thermocouples, become 
steady. The calculation of the results is somewhat tedious, and for details 
the original paper should be consulted. 

The conductivity of layers of thin cotton, wool and rayon fabrics were 
measured by this method, and in each case it was found that the thermal con- 
ductivity increased with the number of layers used. Figure 2 shows the values 
obtained using 1, 2, 4 and 8 layers of cotton fabric, both with the heater above 
and below the fabric. The results can be explained by assuming a constant 
air layer of thickness d,/2 at each plate, the rest of the space being assumed to 
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have a constant thermal conductivity k, figure 3. The measured conductivity K 
will be a function of the thickness d, and will be given by 


kk,d 
hata. © oe ee (6) 
From equation (6) it is seen that Kk as doo. Curve (1) of figure 2 is for 
equation (6) with k,=0-6 x 10-4 c.g.s. ; kR=2-72x10-4 and Gp lel 10. cme 
Curve (2) of figure 2 is for R=2-2 x 10-and d,=1:1x10*cm. Although the 
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Figure 2. Conductivity of layers of cotton fabric as measured by Lees’ disc apparatus. 


curves fit the experimental points, the large value of d, in comparison with the 
single fabric thickness 0-04 cm. gave rise to doubts as to their validity. In Lees’ 
calculation the edges of the material tested were assumed to have the same emis- 
sivity as the metal plates. This equality can be achieved with a solid sample 
merely by varnishing the whole block after assembly with, say, shellac, but when 
using textiles such a course is impossible. Also, when working with textiles, 
heat can be lost from the edges of the fabric by convection of the air; i.e. by 
hot air leaving the centre of the fabric and diffusing out at the edges. Such 
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Figure 3. Boundary resistances. 


effects cannot easily be allowed for in the calculations, and almost certainly con- 
tribute to the results of figure 2. In view of these difficulties, the use of the 
Lees disc to measure the conductivity of textile-air mixtures was abandoned. 
An attempt was made to estimate the thermal conductivity of textile fibres 
by placing them in liquids and measuring the thermal conductivity of the mixture 
using Lees’ disc. If this were done using liquids which had both higher and 
lower conductivities than the fibres, it was hoped that the conductivity of the 
fibres could be estimated. ‘The results obtained using cotton and wool were 
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somewhat erratic, depending on whether or not the liquid was absorbed by the 
fibres, but the value for the conductivity of wool appeared to be 3-5 x 10~* c.g.s. 
units and for cotton 6-0 x 10-4 units. No great accuracy is claimed for these 


results. 


§4. EXPERIMENTS WITH THE APPARATUS OF HERCUS 
AND LABY 


As mentioned earlier, textile fabrics contain a very high proportion of air 
by volume, and in measuring the thermal conductivity of a fabric, it would 
appear better to use an apparatus designed to measure the thermal conductivity 
of gases rather than one designed to measure the thermal conductivity of solids. 
The types of apparatus used to measure the conductivity of gases fall generally 
into two classes : the hot-wire type, in which the heat losses froma heated wire 
in a cylindrical tube are measured, and the parallel-plate method, where the heat 
transmitted through the material from one plate to a colder plate is measured. 
An apparatus of the second type was most suitable for our purpose, and Hercus 
and Laby’s (1919) design has been used. ‘The apparatus, which is shown in 
figure 4, was made in two halves, one unit comprising the top plate, centre plate 
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Figure 4. Apparatus based on that of Hercus and Laby. 


and guard ring, the other unit being the cold plate. Each plate was made of two 
sections which were screwed together. ‘The top sections in the case of the 
top plate, centre plate and guard ring had spiral grooves cut into their lower 
faces, the heating coils being wound into these grooves. ‘The bottom section of 
each of the three plates had a radial groove cut in its upper face to take the thermo- 
couples. ‘he top section of the bottom plate also had a radial groove cut into 
its underside for thermocouples. The cold plate was kept at a constant tem- 
perature by circulating cold water (26°c.) through it from a thermostated 
tank. The apparatus could be used with the hot plate above the fabric and the 
cold plate below or vice versa. Full details of the experimental technique and 
method of calculating the results will be found in Hercus and Laby’s paper 
(1919). Briefly the method is to heat the guard ring and top plate to a temperature 
slightly higher than the centre plate, so that all the heat input into the centre plate 
is transmitted throught the material to the cold plate. Corrections are applied 
for the small amount of heat gained from the top plate and the guard ring. All 
the outer surfaces of the plates-were silver-plated, and when radiation corrections 
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were applied, the emissivity of the plates was taken as 4-1 x 10-12 cal./sec./em?/deg4 


The e.m.f. generated by the thermocouples, which had previously been calibrated, 
was measured on a Cambridge potentiometer. 


(a) Results using air 


As a check on the apparatus, the thermal conductivity of air was determined. 
The air was room air. The results obtained were reduced to ky, the value of the 
conductivity at 0°c. by use of the formula, k;=k,(1+0-003 2), as used by Hercus 
and Sutherland (1934). The values obtained are given in table 2. 

The values obtained with the heater below the air gap confirm the well-known 
result that convection plays no appreciable part until the thickness of the air gap 
becomes greater than about 1 cm., when the effect becomes quite large, as is seen 
by the last value of the table. ‘The mean value obtained for the conductivity of 
air may be compared with some of the most accurate measurements made to 
date. Kannuluick and Martin (1934), using the hot-wire method, give 
ky as 5-76 x 10-* c.g.s. units, and Hercus and Sutherland (1934), using the parallel- 
plate method of Hercus and Laby, give ky as 5-72 x 10-°c.g.s. units. The agree- 
ment of the present results with these values serves as a check on the calibrations 
of the thermocouples and electrical measuring instruments and shows that no 
systematic error was present. 


Table 2 
Hot plate above Hot plate below 
Thickness Thickness 
of layer Ry (c.g.s. units) of layer Ro (c.g.s. units) 
(cm.) : (cm.) 
0-150 5:44 x 10-8 0-171 syste x HOKE 
0-153 5:64 0-184 5°64 
0-450 5-61 0-476 5-84 
0-458 EoD 0-966 5-83 
0-953 5-82 
1-946 6:09 1-968 12-48 * 
Mean 5°75 Mean 5-72 


* Not taken into account in determining the mean yalue. 


(b) Results using fabrics 


The thermal conductivity of layers of thin fabrics made from wool, cotton 
and viscose rayon were determined. ‘The wool and cotton cloths were made as 
nearly identical as possible, but the rayon cloth was thinner than the others. 
The wool and cotton cloths had weights of 4-3 oz. per sq. yd. and the rayon cloth 
weighed 3-1 oz. per sq. yd., the thicknesses under a pressure of 5 gm./cm: 
being 0-037 cm. for the wool fabric, 0-035 cm. for the cotton fabric and 0-023 cm. 
for the rayon fabric. The measurements of conductivity were done with the 
fabrics at room regain,* and with the hot plate both above and below the fabric, 


* Regain is the moisture content expressed as a percentage of the dry weight. 
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Teh2 
Table 3 
Number Bulk density Thickness 104 x conductivity Mean temp. 
of layers (gm./cc.) (cm.) (cal./cm./sec./°Cc.) (E)) 
Heater above fabric : Cotton 
8 0-450 0-252 1-70 27:9 
8 0-450 0-252 1-69 28-0 
17 0:467 0-515 1:74 32-1 
Mean 0-455 Mean 1-71 Mean 29:3 
Wool 
9 0-475 0-310 12S 29°3 
18 0-480 0-612 1-26 32°3 
Mean 0-477 Mean 1:26 Mean 30°8 
Rayon 
1D) 0-520 0-229 1-45 28:7 
16 Oras 0-310 1-45 28-9 
16 0-515 0-310 1:48 27°8 
Mean 0:517 Mean 1-46 Mean 27:5 
All-wool greatcoat material 
il 0:295 0-241 1:05 28-1 
it 0-295 0-241 1-02 29-1 
2 0-285 0-501 1:07 29-4 
2, 0-285 0-501 1:07 29°5 
3 0-287 0-740 1-09 31:6 
4 0:287 0-986 1-13 32°8 
5 0-296 1-201 1-2 33-6 
Mean 0-290 Mean 1:08 Mean 30-6 
Heater below fabric : Cotton 
9 0-460 0-276 1-62 29-4 
17 0-460 0-521 Ae7A 32-6 
17 0-460 0-521 1-65 ies 
Mean 0-460 Mean 1:66 Mean 29°8 
Wool 
9 0-482 0-307 1:17 29-2 
9 0-482 0-307 1:18 27°1 
18 0-468 0-630 1:26 32-0 
18 0-468 0-630 1:28 Silzg/ 
Mean 0:475 Mean 1:22 Mean 30:0 
Rayon 
10 0-531 0-189 1235 29-0 
DS 0-524 0-287 1-40 30-1 
Mean 0-527 Mean 1:38 Mean 29-6 
All-wool greatcoat material 
3 0-290 02735 1-08 33-8 
5 0-287 1-239 1-09 33-3 
Mean 0-289 Mean 1:09 Mean 33:6 
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in order to determine whether or not convection losses through the fabric were 
possible. The values obtained are given in table 3. 

The results show that there is little change in conductivity with increasing 
thickness, and that in most cases where an increase is shown, the higher value 
has been obtained at a higher temperature. The values obtained with the 
heated plate below the fabrics show that no convection took place through the 
fabrics. ‘The bulk densities of the thin fabrics when in the apparatus and at 
room regain were very similar ; the mean values expressed as gm. of fibres per 
cc. were 0-476 for the wool fabric, 0-458 for the cotton fabric and 0-522 for the 
rayon fabric. Since the weaves of the fabrics were similar, it must be concluded 
that the wool fibre has the lowest conductivity of the three materials investi- 
gated. 


(c) Variation of conductivity with density 

In order to be able to make an estimate of the conductivity of the wool fibre, 
it is necessary to determine the curve of conductivity against density for wool 
at a constant regain and then to extrapolate the curve to find the conductivity 
corresponding to a density of 1-3 gm./cc. In order to do this, a length of felt 
was specially prepared with a density of approximately 0-1 gm./cc. and thickness 
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Figure 5. Conductivity against bulk density for wool. 


approximately 3inch. Values at higher densities were obtained by compressing 
layers of the felt together at a temperature of 100°c. and allowing them to cool 
before releasing the compression. If the wool had a regain not less than about 
20%, it was found that the heat and compression treatment caused a permanent 
set. Excess water was not necessary, as the saturation regain of wool at 100° c. 
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is in the region of 20%. Using this technique, densities up to 1-0 gm./cc. 
were obtained. For densities less than 0-1 gm./cc. the felt was carded dry and 
placed in laps between the plates, which were separated by spacing blocks of 
ebonite between the guard ring and the cold plate. The conductivity of each 
felt of density greater than 0-1 gm./cc. was measured at a low (0-7%) and a 
medium (10°7%) regain, the regains being based on the dry weight obtained by 
placing the felt in an oven maintained at 100°c. The values obtained are shown 
graphically in figure 5. 


(d) Variation of conductivity with regain 


The effect of increasing the regain is to increase the bulk density and also the 
conductivity, the magnitude of the effect depending on the bulk density. At 
low densities the regain effect can hardly be detected owing to the small part 
played by the fibres and also as at low temperatures the conductivity of air is 
negligibly altered by increasing the relative humidity. This has been verified 
in the present investigations, which are in line with the results of Gruss and 
Schmick (1928), who showed that introducing 5°%, of water vapour into the air 
caused the thermal conductivity of air to increase by 1%. 
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Figure 6. Effect of regain on conductivity of wool felt of density 0-335 gm./cc. 


The effect of regain on a felt of density 0-30 gm. per cc. when dry was investi- 
gated. ‘The density of the felt varied slightly with regain, and the values of the 
conductivity were corrected, using the curve of figure 5 to correspond in each 
case to a density of 0-335 gm./cc. The values are shown in figure 6, and the 
conductivity is seen to increase linearly with increasing regain. - The slope of the 
line is 7-7 x 10-% units of conductivity per unit change in regain. This value 
will of course depend on the bulk density of the material. The experimental 
points are all within 2°% of the straight line. 


§5. THE CONDUCTIVITY OF HORN 


Since both wool and horn are keratin, the estimated values for the conductivity 
of wool would be expected to be similar to the conductivity of horn. The thermal 
conductivity of horn, in a direction at right angles to the axis of the long chain 
molecules, was determined at various regains by means of the Lees disc apparatus. 
A piece of horn was flattened by boiling whilst clamped in a flat position, and then 
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turned on a lathe into a disc of the same diameter as the Lees discs, and of thick- 
ness 0-226 cm. when dry. The horn warps badly as the regain is altered, and the 
specimen was held clamped flat between the plates of the apparatus, paraffin wax 
being used to ensure good thermal contact between the specimen and the plates. 
The conductivity of the horn was investigated over a wide range of regains, and 
the values are given in figure 7. 

The thermal conductivity of mixtures can in many cases be expressed by a 
formula of the type 

Ry = k(1 ote aC), 


where &,, is the conductivity of the mixture, ky the conductivity,of the pure sub- 
stance, C,,, the concentration of the added substance, and @ a constant. Water 
molecules are absorbed by wool in two states, some being absorbed_on to localized 
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Figure 7. Effect of regain on the conductivity ot horn. 


sites and others existing in the liquid state. ‘The proportion of water in the 
two states at any regain has been worked out by Cassie (1945). Assuming these 
two types of water, the conductivity of horn may be expressed by the equation 


prea Koeln: Gy TSCA x) ate | Lae (8) 


where hy is the conductivity of the dry horn, &,, that of the horn when at regain m ; 
Cy is the concentration of water absorbed on the localized sites, and Ci4_ x) is 
the concentration of water in the liquid state, whilst « and 6 are constants. 

The values of Cy, Cy and C(4_x) for various regain values have been 
evaluated from table 3 of Cassie’s paper, and are given in table 4. 
' The curve in figure 7 is for equation (8), the constants being obtained by 
making the curve to pass through the points k,,=5-3 x 10-4, 4=10-8% and 
k,=6:6x10-*, A=25:2%. The values of the constants are «=49-4 and 
B = 133, indicating that the free water contributes more to the thermal conductivity 


than the water found on localized sites. 
8-2 
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Table 4 
Regain | Cy Cx EGSxz) 
(A) (gm./cc.) (gm./cc.) (gm./cc.) 
% 
3°6 0:045 0-041 0-004 
7:2 0-086 0-082 0-004 
10:8 0-125 Oils 0-012 
14-4 0-165 0-141 0-024 
18-0 0-200 0-159 0-041 
21:6 0-232 0-169 0-063 
25:2 0-260 0-174 0-086 
28°8 0-286 0-174 0:112 
32:5 0-314 0-174 0-140 


§6. INTERPRETATION OF RESULTS 


There are two possible equations which may be applied to the experimental 
points of figure 5 ; they are the Lees formula, equation (5), and Schuhmeister’s 
formula, equation (3). The orientation of the wool fibres in the compressed 
felts used do not agree with the conditions assumed in deducing either of the two 
formulae, and any equation which fits the data must be taken as a purely empirical 
equation without any real theoretical foundation. Curve (1) of figure 5 is for 
Lees’ equation (5), with k, = 0-63 x 10-4, the conductivity of air at the temperature 
of the experiments, k,=3-63 x 10-+, and curve (2) is for equation (5), with 
k,=0-63 x 10 and.k,=3:92x 10+. ‘The two curves, fit’ the experimental 
points quite well, except for densities less than 0-2 gm./cc., and predict the con- 
ductivity of wool at the lower regain (0-7%) as 3-63 x 10-4, and at the higher 
regain (10-7%) as 3:95 x 10-4. The values of the conductivity of horn at the 
same two regains are 4:62 x 10-4 and 5-32 x 10-4. 

Schuhmeister’s formula, equation (3), may be written 


lise = x(RyV} oe kyVp) Hee A hang A 82 (9) 
where x+y=1. 

If x is unity, we have &,, increasing linearly with density from k, to &,, and if 
y is unity we have a hyperbolic relation with k,,, increasing slowly with density 
at first and then more rapidly as the density increases. Equation (9) has been used 
with y=1 for water-vapour diffusion by Peirce, Rees and Ogden (1945), but this 
by no means fits the experimental points, as is seen by curve (5) of figure 5, which is 
equation (9) for y= 1, k, = oo and the density of the fibres 1:3 gm./cc. Asthere are 
no means of determining directly the values of x and y, the only method is to 
take the value of k, as determined for horn and then to fit the equation to the 
experimental points and so determine x and y. This method gives the values of 
« and y as 0-29 and 0-71 respectively. Curve 3 of figure 5 is for equation (6) 
with x=0-29, y=0-71, k,=0-63 x 10-4, ky =4-62 x 10-4, and curve 4 is with 
ky=5°32x10-*. The curves again fit the experimental points within experi- 
mental error. If we take x=} and y=, which are Schuhmeister’s values, and 
fit the equation to the low regain experimental points, k,, = 1-425 x 10-4 C.g.S. 
and p=0-6 gm./cc., we find that k, is 4-04 x 10-4 c.g.s. units, and with this value 
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We get curve (6). This curve is almost identical with curves (1) and (2) over the 
region explored experimentally. 

We have, therefore, to decide which equation, type (3) or type (5), should be 
used. Both equations give curves which fit the experimental points equally 
well over the range measured, but the Schuhmeister formula, with «=0-29 
and y=0-71, takes in the values found for the conductivity of horn and, therefore, 
may be considered to give a better empirical formula than the Lees equation or 
Schuhmeister’s equation with x=4, y=2. 

At densities less than 0-1 gm./cc. neither formula fits the experimental points. 
The apparent increase in conductivity as the density falls below 0-1 gm./cc. is 
probably due to convection. The temperature gradient along a fibre will be 
different from the temperature gradient in the air adjacent to it. This will give 
convection currents whose magnitude will depend on the fibre surface area and 
the average distance between the fibres. As the density increases, the surface area 
increases, but the mean distance between the fibres decreases and a maximum 
would be expected in the (conductivity v. density) curve. This is seen to be the 
case: the convection effect rises toa maximumat a density of 0:01 gm./cc. and then 


Table 
Suieianee oe 104 x thermal conductivity 
(gm./cm”*) (c.g.s. units) 
Polyvinyl chloride-acetate 1:34-1:37 3: 9D 
Polyvinyl chloride 1-2 -1-6 3-9-4-0 
Cellulose acetate (transparent) 1-27-1-32 5°4 
Cellulose acetate (pigmented) 1:37-1-56 5°3=8°/ 
Ethyl cellulose 1-10-1-25 5:6 


decreases as the density increases, and becomes negligible at densities greater than 
0-2 gm./cc. 

Equation (9) with x =0-21 and y=0-79, may be used to give estimates of the 
conductivities of the fibres for the three thin fabrics of table 3. Using the mean 
values of the density and fabric conductivity, the conductivity of the wool fibre 
is found to be 4-6 x 10-4, of the cotton fibres 11-0 x 10~*, and the viscose rayon 
fibres 6:9 x 10-4, the regains being approximately 10, 5 and 10%, respectively. 
Taking the conductivity of air, 0-63 x 10-4, as unity, the relative conductivity. of 
wool becomes 7-3, of cotton 17:5 and of viscose rayon 11:0. The value for wool 
lies between the estimates of Rubner and Schuhmeister, and the value for cotton 
is roughly half their values. ‘The value obtained for the viscose rayon compares 
favourably with those given for cellulose derivatives in table 5 below, although 
the conductivity of viscose—regenerated cellulose—is not known. The result 
for cotton, taken together with the value indicated by the experiments with liquids, 
suggests that Schuhmeister and Rubner’s values are too high. 

The thermal conductivity of synthetic fibres can be measured directly, as the 
materials are readily available in sheet form, and, knowing the fibre conductivity, 
the conductivity of a fabric can be estimated from its density and equation (9). 
The thermal conductivities of many of the synthetic materials used as textiles 
are not available, but a few are given in table 5 (cf. Simonds and Ellis, 1943). 
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No mention is made of the regain of the material at which the measurements 
were made. Wool is seen to have a conductivity comparable with that of the 
artificial materials, whilst cotton appears to have a considerably higher con- 


ductivity. 
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the obvious scepticism of a class of first-year students to a remark that 
water in a falling vessel would not issue from an orifice in the bottom 
of the vessel. 

To demonstrate the truth of the statement, a flask filled with coloured water 
was attached to a spiral spring about 3 feet in length, which permitted of an oscil- 
lation of the flask in a vertical direction (v. figure 1). To the bottom of the flask 
was sealed a short tube, the opening of which was drawn down to a capillary of 
some 2 mm. diameter. The period of oscillation of the flask was about 2 sec. 
and an amplitude of 6 inches could be given to it. Thus it was subjected to a 
vertical acceleration ranging sinusoidally from about 5 ft./sec./sec. upwards to 
the same value downwards, and the value of gravity thus effectively varied from 
(32 +5) =37 to (32 —5) =27 ft./sec./sec. 

This variation, producing, of course, a proportional variation in the pressure 


"Tw experimental demonstration here illustrated was suggested to me by _ 
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(above atmospheric) at the orifice, was sufficient to cause alternately cessation and 
occurrence of flow of the liquid from the orifice, as shown in figure 2. 


Figure 1. Flask containing inky water suspended on spiral spring and executing vertical oscillations 


of period about 3 seconds. Beneath flask is a rotating table covered with white paper. 


¢ Born o 


Figure 2. Shows intermittent character of flow of fluid from opening of fine tube. Due to 


capillary effect at the orifice, the interval during which flow is inhibited is longer than that 
during which it occurs. 
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ABSTRACT. A dynamical method is described for the observation of the advancing 
and receding angles of contact of mercury at very low velocities. It is suggested that the 
irreversible part of the work expended when the line of contact of the liquid traverses 
the solid surface is different for the two directions of motion. 


§4. INTRODUCTION 
TT" angle of contact of mercury differs widely according as the liquid 


is advancing or receding over the surface of a solid. It has already 

been shown (Yarnold, 1938, 1940) that the difference between the cosines 
of the advancing and receding angles of contact against walls of glass, Pyrex and 
quartz is independent of the velocity of the line of contact throughout the range 
from 3 to 30 cm./sec. The hysteresis of the angle of contact at very low velocities 
has now been studied by measurement of the force exerted by a rising or falling 
mercury surface on a suspended sphere which is partly immersed in the liquid. 
An expression for the force experienced by the sphere may be derived as follows. 


Figure 1. 


Consider a sphere of radius a whose centre is at height y above the plane 
surface of the mercury. Whether the mercury surface is elevated or depressed 
in the neighbourhood of the sphere depends upon the value of y and upon the 
motion of the sphere. In figure 1 the mercury surface is represented as depressed 
in the region of the sphere, with which it makes an angle of contact 6. The lines 
joining the circle of contact with the centre of the sphere make an angle ¢ with the 
vertical. ‘The upthrust of the liquid may be divided for convenience into two 
parts : that due to the displacement of the segment ABE and that due to the dis- 
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placement of the cylinder ABCD. The forces due to surface tension act in 
directions inclined at an angle (9+¢—7) to the horizontal. The total upward 
force on the sphere is therefore given by the equation 


U=7 apg (5 —cosd+ ; cos? 4) + 77a*pg sin® 6(a cos ¢—y) 
—2raT sing .sin(6+4), 
where p is the density and T is the surface tension of the liquid. 
If the sphere is immersed completely, the upthrust becomes 


U, =37a°pg. 
Hence, in general, 
OG BRAS lee Sie y 
U, =(5 -- pee q 008 s) + 7a #(coss—2) 
Piet 
see “SINOmSIn (GH d)a- facie (1) 


apg 
(c) (d) 


Figure 2. 
(a), (b). Mercury rising. (c), (d). Mercury falling. 


In this expression for U/U,, the first term is always positive whatever the value 
of ¢. ‘The second term is positive or negative according as cos¢ is greater or 
less than y/a, i.e. according as (6 —¢) is greater or less than 7. Clearly, also, the 
same condition determines whether the effect of the third term is positive or 


negative. 
Let it be supposed that the sphere is immersed so that the mercury 1s exactly 


horizontal right up to the lines of contact (figure 2 (b) or 2 (c)). It follows then 
that (@+¢)=7 and that cosé=y/a. If the upward force acting on the sphere 
under these conditions is denoted by Up, then 


Coes (9 \ oth fy\" 
m=35-3(2) +4(2) nha (2) 


This criterion makes it possible to recognize the point at which the liquid surface 

‘is horizontal right up to the surface of the sphere much more accurately than by 
direct visual observation. We have, therefore, a dynamical method of observing 
angles of contact which is entirely independent of the surface tension of the 
liquid. 
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§2. EXPERIMENTAL METHOD 


Calculation suggests that the most suitable diameter of sphere is about 1 cm. 
In some of the experiments, steel ball-bearings were used ; in others, small 
glass spheres suitably weighted with mercury before being sealed. Each sphere 
in turn was suspended from a Sucksmith ring balance (Sucksmith, 1929) of 
simplified design, by the displacement of which the force, upwards or downwards, 
exerted by the mercury could be measured. ‘The ring was of steel tape approxi- 
mately 0-008 cm. in thickness, 0-4 cm. in width and 5 cm. in diameter. Using 
a single mirror in conjunction with a telescope and a scale at a distance of 2 metres, 
a force of 1 gm. wt. gave a change in reading of approximately 4 cm., while the 
scale could be read to one-tenth of a millimetre. 

All grease was removed from the steel spheres by washing three times in 
ether (B.P.). In some of the experiments the glass spheres were prepared in the 
same way. In other cases they were cleaned in concentrated nitric acid, washed in 
many changes of distilled water, and finally allowed to dry under cover to prevent 
contamination by dust. Redistilled mercury, filtered to give a clean new surface 
at the beginning of each set of observations, was contained in a square glass 
trough of side 7 cm., sufficiently large, that is, to give a truly plane surface. By 
means of three similar calibrated screws, the mercury surface was brought up 
into contact with the lowest point of the sphere, then slowly raised in steps of 
about 0-01 cm. so as to immerse about half of the sphere, and subsequently 
lowered until contact was again broken. ‘The time taken for a total motion of a 
few millimetres was rather more than an hour. The height y of the centre of the 
sphere above the plane surface of the mercury throughout this procedure was 
determined from the rotation of the screws and the small change in the absolute 
position of the sphere as the ring balance was deflected. 

The variation of the quantity U/U, with the ratio y/a during the raising and 
lowering of the mercury is shown by the full curves in figure 3 for the case of a steel 
sphere of radius 0-635 cm. On the same diagram the thin curve represents the 
variation, with the ratio y/a, of the quantity U)/ Uj, as calculated from equation (2). 
The points of intersection of the experimental and theoretical curves give the 
values of y/a when the mercury surface is horizontal right up to the surface of the 
sphere. ‘The angles of contact are then given by the equation 


—cos#=cosd=y/a. 


Clearly the method is capable of considerable precision. 

It will be observed that, on first making contact (y/a=1), the mercury rises 
a short distance up the curved surface, as shown in figure 2 (a), and so exerts a 
downward force on the sphere, the second and third terms in equation (1) being 
negative and together greater than the first. This adhesion of the mercury to 
the solid is only observed with really clean surfaces. Adhesion takes place simi- 
larly when the mercury is falling, but now, the angle of contact being smaller, 
quite a pronounced neck develops (figure 2 (d)). The downward pull on the 
sphere passes through a maximum value, as in Ferguson’s experiments (1913), 
after which contact is broken, and a minute drop of mercury remains at the lowest 


point of the sphere, which may now be as much as a millimetre above the level 
surface of the mercury. 
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§3. RESULTS 


The table shows the measured angles of contact for a number of different 
spheres. Although in any given case the angles of contact may be measured 
quite accurately, the angles themselves are not exactly reproducible, even when 
scrupulous cleanliness is observed. It will be noticed, however, that the receding 


Figure 3. Force on steel sphere. a=0-635 cm. 
1. Mercury rising. 2. Mercury falling. 3. U,/U, (calculated). 


ingles all lie in a comparatively narrow range, while the advancing angles vary 
nuch more widely. The difference between the two angles is greater when the 
lass spheres are washed in ether than when they are washed in nitric acid and 
istilled water. No trace of residue, however, was visible on evaporating several 
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cubic centimetres of ether in a watch-glass. In no case, however, is the difference 
between the cosines of the angles of contact as great as the difference observed 
for a mercury index travelling in a glass tube at a velocity of 3 cm. per sec. or 
faster (Yarnold, 1938). The first group of steel spheres gives angles of contact 
hardly different from those observed with glass spheres cleaned in the same 
way. It must not be inferred, however, that this result is due to the presence 
of a film left by the ether, which forms the true surface in each case, since two 
other steel spheres, obtained from a different source, give much higher advancing 
angles. 

Following Adam and Jessop (1925), it is usual to ascribe the hysteresis of the 
angle of contact of a liquid to a frictional force F per unit length acting at right 
angles to the line of contact. The force F is tacitly assumed to be the same 
whether the liquid is advancing or receding over the surface of the solid, and, 
as a consequence, the cosine of the true equilibrium angle of contact is taken 
to be the mean of the cosines of the advancing and receding angles. ‘The intro- 
duction of the force F, however, is equivalent to the statement that part of the 


Radius —cos § @ (degrees) 
(cm.) Remarks 


Material 
Advancing |Receding | Advancing |Receding 


Glass 0:574 | Washed in nitric acid 0-723 0-600 136°5 127°0 
e 0-522 and distilled water 0-780 0-596 141-0 126°5 


Glass 0:574 | Stored in alcohol and 0-831 0-648 146-0 13025 
3 0:522 washed in ether 0-848 0-600" 148-0 127-0 


Steel 0-635 | Washed in ether 0-860 0-610 149°5 127-5 
5 0397 0-806 0575 144-0 12520 
0°857 5 


55 0-318 0-648 149-0 130- 
Steel 0-556 | Washed in ether 0-930 0-624 158°5 128°5 
Ss 0:437 0-967 0-640 165-0 130-0 


work expended in the formation of fresh surface is irreversible. The assumption | 
that the magnitude of F is independent of the direction of motion of the line of _ 
contact is therefore equivalent to the assumption that the irreversible work 
expended in the formation of unit area of the solid-air interface is the same as the 
irreversible work expended in the formation of unit area of the solid-liquid 
interface. 

Stated in this form, the assumption hardly appears to be justified, though 
the error may not be serious when the two angles of contact differ only slightly. 
Nor does it appear to be supported by the measurements of the angles of contact 
of mercury. ‘The approximate constancy of the receding angle, together with 
the remarkable adhesion of the mercury to the solid surface, suggests that a solid- 
liquid interface of more or less constant energy is formed in all cases, and that 
more or less the same quantity of irreversible work is expended per unit area in 
its destruction. On the other hand, the wide variation in the advancing angle 
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suggests that the irreversible work involved in the destruction of the solid-air 
interface and its replacement by the solid-liquid interface may be extremely 
variable. 
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organizing an exhibition of scientific instruments, thus giving an oppor- 
tunity to British scientific instrument makers to demonstrate once again 
the excellence and ingenuity of their workmanship. 

During the years of war we have of necessity become more self-reliant for our 
research and development in all branches of science and technology upon our 
own native manufactures. And judging by the outstanding results which we 
have achieved during the last six years there can be no doubt at all that our scientific 
instrument makers have served us well. Many of them have turned their efforts 
to the direct production of instruments for aircraft, tanks or ships of all kinds and 
have thus entered the field of semi-mass production or even of mass production 
itself, but this has not prevented them from continuing with the development and 
production of those many scientific instruments of extreme accuracy upon which 
so much of our present research and industrial technique is based. 

Indeed, instead of debasing their standards by introducing volume production 
they have carried into that mass production the skill and accuracy which they had 
long cultivated in their specialized and small quantity productions. 

In the result the war should have put us into a position to expand our scientific 
instrument production upon a large scale, preserving all the excellencies for which 
our leading firms have so justly been famed in the past. 

And the present moment is a most appropriate one in which to demonstrate 
this fact not merely to our own users of these instruments but also to the potential 
customers for them from all parts of the world. 

As in the first world war, so in this last one, we have learnt how backward we 
have been in our research efforts during times of peace. Research is easily and 
widely recognized as the sole basis for rapid progress and superiority in wartime ; 
so much so that every effort is made by the State to entrain every scientist available 
in one line of research or investigation or in another. ‘The scientist is the blue- 
eyed boy called in to help the solution of every problem and to force the pace of 
progress in a never-ending duel of wits and intelligence with the enemy. 


To Physical Society is much to be congratulated upon its enterprise in 
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And the scientist must be duly and properly equipped for his work. Great 
research establishments are improvised regardless of cost and filled with the most 
up-to-date instruments and machinery—which, incidentally, put to shame the 
accommodation and facilities of our peacetime universities and other research 
centres. 

Nothing is regarded as wasted that may contribute to the superiority over the 
enemy inany field. Tens of thousands of research man-hours are gladly sacrificed 
to the merest chance of some new development which will give us advantage over 
the foe. 

But as we experienced between the two wars, the atmosphere is apt to change 
when the crucial danger of war is past. We then weigh up the£ s. d. and evaluate 
whether this or that expenditure is worth our while. And in that calculation we 
tend to become shorter and shorter sighted. We are liable to forget that progress 
and superiority in peace are as much based upon research as they are in war. 
The tempo of our drive forward on the scientific front slows down and we fail to 
educate and train a sufficiency of highly skilled staff to maintain the intensity of 
our research. 

That at least was the experience after the last war, and it is a folly that we must 
at all cost avoid repeating after this war. 

Research is the life-blood of our industrial progress, and industrial progress is 
the basis for our continuing prosperity as a great exporting nation. 

To attain the degree of research activity which is essential to our future as a 
great inventive industrial nation four things are required. 

First we must have the men and women adequately trained and in sufficient 
numbers. ‘Thatisa matter of more and better educational facilities of all kinds and 
a free access to those facilities for the intelligent youth wherever they may be and 
whatever kind of homes they come from. 

Second we must raise the status of our scientists. ‘Too long they have been 
regarded as rather superior craftsmen definitely inferior to the leaders of industry 
or the principal administrators. ‘The Classical Tradition in education still claims 
a superiority which is out of date, and the scientist is still looked upon by many as 
of a lower grade. Whether we look at the salaries they are offered or the position 
to which they are relegated in the hierarchy of industry or of the State Service. we 
see that they are kept down in a position of inferiority. 

I am not one of those who claim that the scientists should rule in all spheres, 
nor am I a believer in technocracy, but I am firmly convinced that both scientists 
and technicians are worthy, and have amply proved themselves worthy, of a 
position equal to that of any other section of the community. 

But until we gain the recognition of this fact by the country as a whole we shall 
not get the best out of our scientists, nor shall we make the headway that we should 
in applying the results of our research to our every-day needs. 

‘The scientist must not be regarded as an individual to be called in to assist in an 
emergency when matters have got into a hopeless position—like the plumber with 
a burst water-pipe—rather he should be looked upon as a valuable and equal co- 
operator in every phase of the solution of our problems of living. 

Third, we require proper accommodation and facilities for our research. It 
is not necessary that every research institution should be a palace with chromium- 
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plated ware, nor do I believe that string and sealing wax provide the essential 
setting for all that is of value in research. There is a happy mean in which men 
and women can do their best work, provided with the most up-to-date instruments 
that they require. 

A good many of our research establishments today fall far below that mean in 
standard both as regards space—and adequate space is vital to good research—and 
as regards equipment. 

I hope that full advantage will be taken of wartime surpluses to see that all our 
universities and other research establishments are fully equipped with adequate 
machinery and instruments. 

Finally, there is the all-important matter with which this exhibition is con- 
cerned. We cannot afford to be dependent upon others for the instruments 
necessary for our education of scientists or for the research work which they do, 
nor, indeed, for all those many testing instruments which are today part and 
parcel of every up-to-date industrial concern. It is essential that we should 
have and maintain in this country a pre-eminent scientific instrument manu- 
facture both for our own use and to help to supply those many other countries 
who cannot themselves afford to set up such manufactures. 

I believe that if we could insist upon these four points we could lay the 
permanent and solid foundation of a live, active and vigorous research in our 
country which would enable us to keep in the van of the great industrial countries 
of the world. 

It is because I believe that this exhibition will help along the lines that I have 
sketched that I am so delighted to have the honour today of declaring this 
exhibition open. 


REVIEWS OF BOOKS 


Elementary Wave Mechanics, by W. HEITLER. Pp. viiit+136. (Oxford and 
London: Sir Humphrey Milford, at the Oxford University Press, 1945.) 
7s. 6d. . 


In this little book, Professor Heitler gives the substance of a course of lectures which 
he delivered in Dublin. It is truly elementary, both in the demands which it makes on 
the previous knowledge and mathematical ability of the reader, and also in that it deals 
only with the elements of wave mechanics. ‘There is here no metaphysics, but a plain 
introduction to the physical arguments which led de Broglie and Schrodinger to the wave 
equation. The uncertainty principle is introduced as a demonstrable theorem, and not 
as a vaguely supported principle, and is indeed at first named the uncertainty relation, 
a change in nomenclature which might be advantageous if generally adopted. 

The wave equation is used to investigate the hydrogen atom, the Zeeman effect and 
the two-electron problem in a manner which is not now unfamiliar in elementary texts. 
The Pauli exclusion principle is introduced in the last-named chapter. 

The best part of the book, however, is in the last three of the nine chapters, where 
the periodic system, the theory of the homopolar chemical bond and the general problem 
of valency are taken up. - In conformity with the plan of the book, these matters are 
treated in part in qualitative fashion, with guidance from the theory of simple cases, but 
the treatment is so lucid that these chapters must be of value to all who are concerned 
with the problems of chemical combination. The explanation of how it comes about 
that valencies can become saturated after a limited number of atoms have combined, 
or of how the carbon atom lends itself to the multitudinous compounds in which it is 
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found, are of more value for a clear picture than would be many pages of detailed algebra 
working out a more complete theory. 

In reviewing the same author’s Quantum Theory of Radiation, the view was expressed 
that the author would do a service if he wrote a more elementary introduction to the 
subject. Can he now put all students in his debt by writing a book intermediate in 
standard and difficulty between this very elementary treatment of the simplest wave 
mechanics and that very difficult work in which the whole range of quantum mechanics 
is called on ? Such an intermediate treatment would, presumably, deal fairly exhaustively 
with the mechanics of particles and would have some introduction to the subject (untouched 
in the elementary book) of the interaction between matter and radiation, but without 
attempting an exhaustive treatment of this more difficult topic. jens Baye 


Questions Actuelles de Géophysique Théoretique et Appliquée, par EDMOND 
Rorue. (Paris : Gauthier-Villars, 1943.) 275 Fr. 


Written and published during the German occupation of France, this book has only 
recently reached Great Britain. It covers a very broad field, and from the data that have 
been obtained in the realms of gravity observations, radioactivity of the earth’s crust, 
volcanology, seismology and terrestrial magnetism an effort is made to present a unified 
picture of the modern conception of the earth’s constitution. Isostatic compensation is 
discussed in the light of gravity observations, and is supported by some detail of the work 
of the Swedish geologist de Geer on varved clays and the retreat of the ice-sheet, leading 
to the subsequent elevation of the land masses due to the removal of the ice-load. The 
two following chapters on radioactivity and volcanic action are broadly discussed in relation 
to isostasy, so that these branches are integrated into a self-consistent whole. 

In view of the unity up to this point, it is rather surprising to find the-section on 
seismology so short, and in particular the importance of near earthquakes, and the informa- 
tion revealed by them on the sequence and thickness of layers in the crust, so closely 
allied to the previous aspects, is scarcely mentioned. ‘The chapter on terrestrial magnetism 
is also disappointing and does not appear to fit the main theme established in the early 
part. After a brief description of the earth’s field, and ending with the diurnal and annual 
variations, the main part is devoted to a speculative theory, due to Haalek, on the cause 
of the earth’s magnetism. ‘The work on isopors and its possible implications is ignored. 
At the end of each chapter a short section describes some of the practical applications to 
geophysical prospecting, and here the contents supplement a previous publication, devoted 
to these aspects, by the same author. 

A number of numerical errors occur in the text ; but in the majority of cases they 
are obvious, and will, no doubt, be corrected in later editions. Some of the illustrations 
leave much to be desired, the reproduction is poor and a few are not easy to follow. In 
figure 67, for example, the axes are not labelled, but the fortunate inclusion of a table, 
giving the same data, allows them to be disentangled. In spite of these defects the book 
is a good review of the present ideas on the subjects treated. J. M. B. 


The Mathematics of Ultra-high Frequencies in Radio, by L. N. BriLLourn. 
Pp. vi+210. (Brown University, Providence, R.I., 1943.) Price not. 
stated. 


The theoretical background to micro-wave radio has become of increasing importance 
during recent years, and the present book, which is a record of a course of lectures given 
at Brown University in 1943, is a welcome addition to the few existing works dealing 
with the theory. It provides a unified account both of the physical concepts associated 
with wave propagation and the mathematical methods used in analysing the electro- 
magnetic field. It is readable by the physicist and engineer as well as the mathematician - 
and by seriously studying the book, the student and the research worker will avoid the 
necessity of consulting a large number of the original papers which lie scattered throughout 
the scientific literature of the past decade. 

The book begins at an elementary level. Tuning circuits with lumped or with dis- 
tributed constants are discussed, and the usual introduction to phase and group velocity 
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is given. ‘Transmission lines and co-axial cables are studied and there is a section on 
the various types and relationships of units. Acoustic waves are studied preparatory to 
electromagnetic waves in pipes. 

Then the mathematics proper is introduced. The author assumes slight initial 
knowledge on the part of the reader and begins with a discussion of vector algebra, scalar 
and vector fields and the expression in curvilinear co-ordinates of various vector functions. 
Maxwell’s equations are set up and solved by the method of the retarded potential. The 
expression of these equations in tensor form (fundamental from the relativistic point of 
view) is given, and the Hertz potential and Poynting vector are discussed. A short chapter 
is devoted to the skin effect. 

The propagation of waves along a cylinder of arbitrary cross-section is studied 
mathematically in some detail ; the treatment is general at first, but particular cases 
discussed include tubes of rectangular, circular or elliptical cross-section and co-axial 
cables. All these cases lend themselves readily to a mathematical analysis, and the various 
types of possible waves are classified. 

The last chapter is given to cavity resonators. It is a pity that the treatment is not 
more comprehensive, and while attenuation, damping and the distribution of the higher 
modes are discussed, no mention is made of some of the practical problems in resonatoi 
theory, such as the calculation of resonant frequencies and losses in a compound resonator 
by the (numerical) relaxation method, or, in special cases, by the Hahn method of 
matching solutions. Similarly, in the chapter on wave guides no mention is made of 
discontinuities in the cross-section, although this has been the subject of discussion in 
several papers appearing in recent literature. It is disappointing to find no mention of 
antenna theory (except the reciprocity theorem). But here again recent work by the 
author on this subject may be found in the literature. 

The lecture notes are well presented and the student should not find it difficult to 
read them. It is gratifying to note the enterprise of Brown University in making available, 
in mimeographed form, the results of the Advanced Instruction and Research in Mathe- 
matics. L. S. GODDARD. 


mcience i Progress, Fourth Series.’ Pp. 331..-(Yale: University Press ; 
London : Sir Humphrey Milford, at the Oxford University Press, 1945.) 
$3.00. P 


This delightful set of essays is made up of contributions from eleven very distinguished 
specialists in various branches of knowledge, edited by George A. Baitsell and with a 
foreword by L. L. Woodruff. In his thoughtful foreword Mr. Woodruff says :—‘‘ The 
accumulation of knowledge during the millennia since man’s emergence from a condition 
of complete dependence on his environment to the. relatively masterful position that he _ 
now holds is, except in its broad sweep, beyond the compass of one mind, however talented. 
The polyhistor has been succeeded by the specialist.”” This is indeed true, and for this 
reason the reviewer feels quite unable to do much more than give a broad picture in bold 
outline of the various subjects presented. 

The contributions are mainly the substance of national lectures delivered under the 
auspices of the organization known as ‘“‘ The Society of the Sigma XI”’, which is devoted 
to the encouragement of research in science, and the book represents the fourth series of 
such lectures. 

The collection is intended to provide scientific men with general outlines, although 
accurate ones, of fields of research in which, while not actually working, they may be 
interested. From this point of view the purpose in view has been excellently achieved. 

Each author has provided a most excellent series of references attached to his chapter, 
and it is flattering to national vanity to see that so many Englishmen are mentioned as 
having pioneered or followed the exciting hunts for the new knowledge which goes to 
the making of the various chapters. : 

A striking war-time chapter is the first, by Professor Miles of Yale. The importance 
of the psychological approach to certain problems in military aviation has been proved 
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to be of paramount importance, and the interested reader will learn how an aviation cadet 
is examined from all points of view so as to ensure not only that he possesses the personal 
qualities so essential for an airman’s success, but also that he receives the training to develop 
fully these qualities. 

In the second chapter Professor Bronk of the University of Pennsylvania speaks also 
of war-time interests in the combined fields of Physics and Biology. He gives attention 
to the subject of sensory impulses from a pressure-sensitive nerve ending, and there is 
much interesting matter belonging to past history and the live present in regard to those 
rhythmic oscillations of potential in nerve fibres which are so important in impulse pro- 
pagation along nerves. 

Professor Hecht of Harvard follows with an extraordinarily interesting account of 
the amount of energy required in the process of vision. He explains the experiments 
which have led to the remarkable result that for a normal individual to perceive a faint 
flash of light on the threshold of visibility the minimum energy lies between 5 and 8 quanta. 

Professor Loewi of New York University, formerly of Gratz, gives a fascinating 
account of the discovery of the chemical transmission of nerve impulses. It is, indeed, 
remarkable that the nervous system produces so many of its effects through the liberation 
of chemical substances. ie 

From the pen of the late Professor Birkhoff comes a masterly historical account of 
some aspects of mathematical physics garnished with personal recollection and luminous 
with individuality—it makes excellent reading. 

Professor Debye of Cornell gives, as is to be expected, a first-rate résumé of modern 
methods of obtaining very low temperatures and some of the problems to which these 
methods have provided a key. The subject of his choice is clearly not a closed book, 
for he points out the possibility of extending the experimentally accessible range towards 
the absolute zero by removing the disorder which may exist in the nucleus itself and 
connected with nuclear magnetism, just as the methods he described have nearly removed 
all the disorder associated with atoms and their motion. 

Professor Eyring of Princeton contributes a valuable article on physical chemistry, 
including a consideration of the distribution of energy in a molecular system, rates of 
reaction, and luminescence inhibition. 

Professor Rabi of Columbia discourses on molecular beams and ‘“‘ radio-frequency 
spectroscopy ’’. 

Dr. Hickman, who is Director of Research of the Distillation Products Company, 
and is an Englishman, produces a valuable chapter on vacuum chemistry—in particular 
molecular stills and their applications, for example in the vitamin A industry. The 
new substance kitol, which is the mirror image in properties of B-carotin, itself the normal 
precursor of vitamin A, was discovered with the use of the vacuum still. 

Professor Elvehjem of Wisconsin gives a detailed and convincing account of the 
vitamin B complex. 'The formulae of the various substances, about a dozen, are discussed 
where possible and the dietary factors explained—altogether a most illuminating article. 

Professor Cohn of Harvard University, in this last chapter, gives a most interesting 
account of work on which he and his colleagues have been engaged during the war on the 
chemistry of proteins, particularly in the chemical interactions and the separation of the 
protein components in protein solutions—researches which involved the proteins of blood 
plasma, so important in the alleviation of human suffering. A great deal of information 
appears to be concentrated in a few pages, from the relative dimensions of various proteins 
through the chemical fractionation of the components of plasma to details in connection 
with the regulation of bodily function by hormones and enzymes. 

The impression gathered from reading through this remarkable collection of masterly 
essays is one of immense satisfaction at being privileged in one small volume to read 
such clearly presented accounts of current research in so many widely separated fields. 
The references at the end of the book will be valuable to anyone whose curiosity is aroused 
sufficiently to follow up particular points ; the excellent diagrams and photographs are 
‘fully illustrative and give the impression of great accuracy. If one may intrude a word 
of friendly envy, it is to say what a rare pleasure it is to handle a book produced so lavishly 


and SO perfectly on such beautiful paper. If it can be bought in this country, it is 
magnificent value for three dollars. R. WHIDDINGTON. 
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Chemical Crystallography, an introduction to optical and X-ray methods, by 
C. W. BUNN. Pp. xii+422 and 13 plates. (Oxford : Clarendon Press, 
19452) 025s. 


The book is addressed to “‘ students of chemistry who wish to acquire some knowledge 
of crystallographic methods and research workers who wish to make practical use of such 
methods”’. It is divided into two sections, the first dealing with the identification of solid 
substances by physical methods and the second with the determination of crystal structures. 
After an introductory survey, the author begins with a chapter entitled The Shapes of Crystals, 
whichgivesan introduction to crystal morphology and the relation between shape and struc- 
ture. Heapproaches the subject of crystal symmetry through the phenomena of crystal growth 
and achieves a simplified treatment of elementary crystallography, e.g. indices of crystal 
planes by means of the unit cell rather than in terms of axial ratios. This is followed 
by an account of crystal optics and the properties of crystals observed in the polarizing 
microscope, leading up to the procedure which is then given for the identification of trans- 
parent crystals by microscopic methods. ‘The last chapter of section one introduces the 
elementary theory of diffraction of x rays and is concerned with the practical details of 
taking x-ray powder photographs for identification. 

Over half the book is devoted mainly to the determination of the arrangements of 
atoms in crystals. This second section opens with an account of the single-crystal 
techniques, how unit cell dimensions are measured and how this knowledge is applied, 
amongst other things, to identification and determination of molecular weight. Succeeding 
chapters give a masterly treatment of the principles involved in the determination of 
crystal structures by trial and error, and numerous examples of both organic and inorganic 
compounds are worked out in detail. A chapter on electron density maps and vector 
maps includes an account of Fourier series methods, of the use of optical methods in place 
of calculations and the scope of vector methods. The author also discusses the evidence 
on crystal structure given by optical, magnetic and other physical properties, and concludes 
with a brief account of broadened x-ray reflections and diffuse spots. 

The well-balanced review of progress in method and technique since the publication 
of The Crystalline State, by W. L. Bragg, makes this an indispensable companion io that 
volume. Formal physical and mathematical treatment are avoided, but the book has a 
discipline of its own and demands careful reading. The author assumes more than a 
passing interest in both the organic and inorganic compounds, thus revealing the range 
of his own contributions to the subject. He has not hesitated to adopt an unorthodox 
development in order to reflect the course of his own experience in applying crystallo- 
graphic methods to chemical problems. 

The text is profusely illustrated and provided with useful appendices and excellent 
indexes. Both author and publisher are to be congratulated on producing so fine a book 


under war conditions at a reasonable price. F. A.B. 


Polarographic and Spectrographic Analysis of High Purity Zinc and Zinc Alloys 
for Die Casting, by various authors. Pp. 117. (H.M. Stationery Office, 


1945.) 5s. 


An account is presented of researches, carried out under a British Standard Institution 
Panel, to investigate the use of physico-chemical methods for the analysis of high-purity 
zinc and zinc alloys for die-casting. 

The work concerned has proved the basis for a British Standard Method on the same 
subject, and the present publication is in the nature of a detailed discussion of alternative 
techniques, rather than a working manual for use in the laboratory. Such being the case, 
its primary appeal is to the specialist and to the spectrographer in particular, but the body 
of evidence presented is sufficient to warrant consideration by all interested in analysis 
in the wider sense. In view of the emphasis in modern spectrographic research on the 
investigation of different excitation methods, two papers giving considerable detail of varia- 
tions on simple condensed sparks and an intermittent A.C. arc technique are of con- 
siderable interest. In a further section an arc technique is described in which the sample 
is placed on a suitable electrode in the form of oxide. This procedure allows of the 
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preparation of standards of known composition, enabling the spectrographer to be 
independent of chemical analysis. In this paper methods of plate calibration are dis- 
cussed in some detail, but at the same time it is agreeable to find the results presented in 
such a manner as to be readily intelligible to a non-specialist reader. 

The description of an alternative method of attack by the polarograph will be of interest 
to laboratories without spectrographic equipment, and also as providing an independent 
result by an entirely different technique. It is shown that results of considerable precision 
can be achieved by modifications in chemical pre-treatment, suitable to the particular 
problem in hand. 

In fact, the standard deviations of the results by all the methods discussed agree very 
closely and provide important additional evidence, not in vague generalizations but in a 
precise and concrete manner, of the value of physical methods to the critical analyst. 
Finally, it may not be out of place to mention that the amount of information presented 
for a modest sum compares exceedingly favourably with the usual scientific text-book. 

H. G. SHORT. 


Radiocristallographie, by ANDRE GUINIER. Pp. ix+294, with 15 plates. (Paris: 
Dunod, 1945.) Price not stated. 


This can be recommended as a really excellent book, the sort of book that has 
been badly needed for several years. It is written specifically for those who wish 
to use the study of crystals by x-ray methods as a means to some further end. It is not 
meant for specialists in crystal-structure determination. Yet all will benefit by reading it, 
especially those who lightly play with structure factor tables, Fourier series and calculating 
machines, without ever troubling to understand the foundations upon which the subject 
is built. True, it dismisses gas tubes in a contemptuous seven lines and gives 1 in 10,000 
as a desirable limit for precision work with powder cameras. But it introduces the reci- 
procal lattice where it ought to be introduced—right at the beginning ; it relegates all 
tiresome but necessary mathematical proofs to appendices ; it show the raison d’étre of 
some of those mysterious tables in.the International Tables for Crystal Structure Deter- 
mination that we ourselves never seem to need ; and it lays proper emphasis on the study 
of crystal texture andimperfection. Every theorem is ilustrated by practical and interesting 
examples ; and numerous useful hints about apparatus and methods are given from the 
author’s own wide experience. 

Yet the book gives a bad first impression. The pages are uncut and the title is not 
clear to an English technical reader. It sounds as if it might have something to do with 
wireless receiving sets, but the prospective buyer is left to guess. After all, x rays are not 
entitled to a corner in the expression “ radio”’. The cover is of thin paper of a bilious 
shade of yellow (or is the reviewer perhaps jaundiced through having read the book during 
an attack of gastric influenza ?). The paper is of extremely poor quality and bad colour ; 


> 


the type is crowded and difficult to read. There is a table of contents, but no subject or 


author index, and the references are not up to date (for instance, the reader is referred — 


to the 1932, instead of the 1940, edition of Clark’s Applied X-Rays’’). Finally, and most 
unkindest cut of all, the present writer’s name is misspelt. It was Robert Louis Stevenson 


who, when some American publisher (before the days of lease-lend) had pirated one of his — 


books, said that it was not the theft that angered him so much as the mis-spelling of his 
name. “I saw my book advertised as the work of R. L. Stephenson ”’, he says, “ and I 


own I boiled. It is so easy to know the name of a man whose book you have stolen, for — 


there it is, full length, on the title-page of your booty. But,no, damn him, not he! He 
calls me Stephenson ”’. 

But in this case the author and publisher must be forgiven. The book is a revelation, 
not only of the difficulties of writing and producing in occupied France, but of the 
indomitable spirit that would undertake such a work at such a time and under such 
conditions. And, with all its faults, it is certainly a book that will be referred to again and 
again. ‘The pity is that it will fall to pieces in the process ; but surely we may soon 
expect a more worthily produced second edition, with references brought up to date and a 
modernization of the text where necessary ?) Or cannot some British or American pub- 


lishing firm negotiate for the translation rights and produce an agreed English text in a ~ 


better setting ? The book deserves it. K. LONSDALE, 
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THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


(A NEW DEVELOPMENT IN 
HEARING AID TECHNIQUE.) 


We believe that we have made considerable advances 


in understanding the most complicated problem of 


fitting deaf people with hearing aids with a perfor- 
mance most suitable to their needs. In the past a great 
deal has been written and said of taking an audio- 
metric curve of a patient and adjusting the hearing aid 
to the curve with the same certainty as an optician deals 
with a prescription for lenses. 


Those familiar with the subject know that the problem 
is actually of great complexity, and it is possible for 
two people with identical audiometric curves to react 


_ quite differently to a particular type of hearing aid. We 


have given the name of ‘‘ Differential Compression ”’ 
to the type of amplification in which the relative 


‘amplification of the higher speech frequencies and the 


rest of the scale is a function of the total sound level 
handied by the amplifier. This it is believed will make 
it much easier to fit difficult inner ear deaf cases with 
hearing aids, and be of benefit to children in special 


schools for the deaf. We have incorporated this — 


feature in some of our wearable instruments as well as 
in new group hearing aid installations for schools. 
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